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I.  SUMMARY 

The  research  in  liquid  rocket  motor  combustion  instability 
at  Princeton  University  consists  of  a  number  of  closely  interrelated 
approaches  to  the  problem.  The  following  summary  covers  the 
theoretical  studies,  the  basic  experimental  investigations  and  the 
rocket  motor  parametric  studies  in  the  order  in  which  they  appear 
in  this  progress  report. 

Further  calculations  based  on  the  time-lag  model  of 
longitudinal  mode,  nonlinear  instability  presented  in  Technical 
Report  No.  677  have  been  performed.  These  calculations  are  presented 
in  Section  I1I-B  and  show  that  stability  criteria  vary  with  mean- 
flow  Mach  number  and  mode  number.  A  third  model,  in  addition  to  the 
two  previously  presented6,  has  been  developed  (see  Section  III-C). 
This  model  considers  both  the  presence  of  shock  waves  in  the  longi¬ 
tudinal  mode  and  small  time-lag  or  phase  between  pressure  and  energy 
addition  at  the  concentrated  combustion  zone.  Periodic  solutions 
of  finite  amplitude  are  found  only  in  operating  regimes  where  the 
steady-state  is  unstable  so  that  the  "triggering"  of  an  oscillation 
by  a  finite  disturbance  is  not  possible  when  the  phase  is  small. 

This  is  consistent  with  the  results  of  Technical  Report  No.  677, 
where  "triggering"  was  found  to  be  most  probable  for  a  critical 
phase  range,  and,  based  on  the  trend  of  numerical  calculations,  was 
impossible  in  the  limit  of  phase  going  to  zero. 

In  fact,  all  of  the  results  of  the  three  models  of  longi¬ 
tudinal  mode  nonlinear  instability  and  the  linear  model  are  shown 
to  be  perfectly  consistent.  The  wave  forms  vary  considerably 
depending  upon  the  various  steady-state  parameters.  Singularly 
important,  in  this  respect,  are  the  effects  of  above-resonant  versus 
below-resonant  oscillations.  (The  off-resonant  conditions  are 
results  of  the  time-lag). 

The  case  of  nonlinear  transverse  instability  is  being 
analyzed  presently.  An  important  part  of  this  theoretical  program 
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is  the  development  of  a  nonlinear  boundary  condition  at  the  nozzle 
entrance.  This  must  be  obtained  by  an  analysis  of  the  unsteady 
nozzle  flow.  This  analysis  has  been  carried  out  to  second  order 
and  is  reported  in  Section  I1I-D.  The  nonlinear  nozzle  analysis 
is  essentially  an  extension  of  the  linear  treatment  of  the  unsteady 

12 

nozzle  flow,  the  numerical  results  of  which  will  soon  be  available. 

Pulsed-limits  testing  to  determine  the  threshold  level  for 

initiation  of  nonlinear  transverse  combustion  instability  continues 

to  provide  important  data  on  trends  in  stability  for  the  parameters 

Involved.  In  Section  IV-A  LOX,  alcohol  rocket  tests  are  reported 

using  unbaffled  injector  designs  which  further  substantiate  the 

Important  role  played  by  inter-spud  spacing  in  reducing  the  tendency 

toward  nonlinear  instability,  especially  in  the  space  engine  range 

18 

of  150  psia  chamber  pressure.  Thrust  per  element  and  spacing 
concepts  are  viewed  as  being  closely  related. 

Elimination  of  propellant  injection  in  certain  locations 
along  the  injection  circle  were  indicated  to  be  effective  stability- 
wise:  if  1)  the  propellant  injection  gap  was  sufficiently  large 
resulting  in  behavior  similar  to  certain  baffle  effects,  or  2) 
where  the  gap  occurred  often,  resulting  in  effects  similar  to  those 
experienced  with  spud  orientation  changes.*-  Higher  pressure  (300 
psia)  operation  in  each  case  increased  the  incidence  of  linear  in¬ 
stability.  The  further  increase  of  chamber  pressures  to  600  psia 
and  above  produced  higher  mode  oscillations.  This  was  found  to  be 
a  very  familiar  stability  pattern. 

Testing  with  LOX,  RP-1  indicated  basically  similar  stability 
behavior  to  that  experienced  with  LOX,  alcohol.  However,  higher 
velocity  LOX  injection  and  lower  velocity  fuel  injection  tests  with 
the  1  x  12  injector  have  provided  increases  in  the  instability  regimes 
over  the  equal  velocity  injector  design.  Noticeable  stability  im¬ 
provement  has  again  been  provided  through  use  of  the  6x2  injector 
design.  Positional  and  directional  effects  noticed  in  the  pulsing 
studies  initiated  relative  mass  and  mixture  ratio  surveys  of  this 
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lnjector  design.  Marked  changes  in  the  deviations  of  relative  mass 
and  mixture  ratio  as  well  as  positional  and  directional  effects 
were  shown  to  occur  as  the  mixture  ratio  was  altered. 

Longitudinal  stability  testing  on  the  square-motor  which 
has  Included  changes  in  the  propellant  used,  injector  element 
spacing,  orifice  diameters  (and  hence  droplet  site  and  relative 
velocity),  numbers  of  elements,  and  spray  fan  interaction  have  shown 
certain  important  relationships  as  discussed  in  Section  IV-B. 

As  Injector  orifice  sizes  have  increased  (with  a  decreased 
number  of  Injection  elements)  the  incidence  of  instability  has  been 
altered.  Starting  from  a  linearly  stable  design  that  was  provided 
by  the  6x6  injector,  then  moving  to  a  range  in  which  small  pulses 
were  required  to  trigger  the  instability  (4x4  design),  a  distinct 
linearly  unstable  regime  was  finally  reached  for  the  largest  (.12(7') 
diameter  orifices.  The  4x4  nonlinear  instability  regime  was  the  largest. 

RP-1  was  also  used  with  LOX  in  the  later  case  showing 
similar  stability  behavior  as  the  alcohol. 

In  this  larger  orifice  injector  case,  pulsing  has  proven 
to  be  of  only  limited  effectiveness  in  extending  the  stability 
regime. 

Increasing  the  LOX  orifice  diameter  while  leaving  the 
RP-1  orifices  untouched  produces  more  rapid  combustion  with  peak  c* 
reached  at  a  shorter  chamber  length.  Instability  tendencies  were 
increased  via  this  change.  Larger  RP-1  orifice  tests  with  the  LOX 
unchanged  produced  slower  burning  conditions  and  eleminated  lonti- 
tudinal  mode  instability  for  the  chamber  lengths  tested  (to  38%  inches). 

The  study  of  the  displacement  effect  using  freon  shower- 
head  injection  into  a  "pseudo"  rocket  chamber  (see  Section  IV-C)  has 
yielded  the  following  data: 

1.  Tangential  mode  oscillations  of  less  than  50  psi 
peak-to-pcak  have  resulted  in  important  displacement  effects. 

2.  Not  only  have  large  quantities  of  vapor  been  cyclically 
displaced  but  stream  breakup  has  also  been  observed. 
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3.  Pressure  and  velocity  measurements  substantiate  the 
frequency  and  spinning  mode  character  of  the  oscillations. 

Another  basic  experiment  probing  the  effects  of  oscillating 
pressure  and  velocity  conditions  on  droplet  distributions  (Section 
IV-D)  has  indicated  that: 

1.  The  number  of  droplets  present  at  a  given  axial  and 
radial  location  in  the  resonating  chamber  changes  between  steady - 
state  and  oscillating  conditions. 

2.  Relatively  small  amplitudes  less  than  20  psi  peak-to- 
peak  are  capable  of  noticeably  disturbing  the  droplet  distributions. 

3.  The  photographic  techniques  evolved  allow  droplet 
sizes  down  to  10 diameter  to  be  measured.  This  direct  photo¬ 
micrography  approach  was  found  necessary  when  the  light  scattering 
technique  was  found  incapable  of  supplying  the  required  droplet 
characteristics  as  reported  in  Technical  Report  No.  648. 

Prediction  of  the  characteristic  shape  of  the  lower  length 
limits  of  longitudinal  combustion  is  found  in  Section  V.  This 
represents  one  application  of  theoretical  study  of  droplet  burning 
presented  in  Technical  Report  No.  671,  together  with  further  data 
based  on  a  model  of  unsteady  burning  in  the  droplet  wake.  These 
predictions  agree  well  with  previous  longitudinal  experimental  data 
indicating  trends  in  the  stability  limit.  Used  in  conjunction  with 
the  sensitive  time-lag  theory,  this  approach  holds  promise  in  allowing 
instability  regimes  to  be  predicted  from  calculated  droplet  sizes. 

Such  correlation  is  being  investigated. 

The  original  time-lag  theory  of  Crocco  and  Cheng  assumed 
that  the  mean-flow  Mach  number  in  the  chamber  was  small  and  its 
square  was  negligible.  That  assumption  does  not  apply  to  present-day 
rocket  chambers  of  low  contraction  ratio  design.  Therefore,  that 
theory  is  being  extended  to  include  higher  order  Mach  number  effects. 
Preliminary  results  for  the  corrected  stability  limits  in  the  case 
of  concentrated  combustion  are  presented  in  Section  VI. 
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II.  INTRODUCTION 

This  is  Che  fourth  in  a  series  of  yearly  progress  reports 
on  the  continuing  investigation  under  NASA  Contract  NASr-217 
(formerly  NASA  Grant  NsG-99-60)  into  the  general  problem  of  combustion 
instability  in  liquid  rocket  motors.  In  these  studies,  the  emphasis 
has  been  placed  on  the  nonlinear  aspects,  both  theoretical  and  ex¬ 
perimental.  The  research  involves:  theoretical  Investigations  of 
the  unsteady  combustion  process  and  the  mechanisms  involved,  rocket 
thrust  chamber  and  nozzle  acoustic  characteristics  and  the  inter¬ 
action  between  these  two  important  parameters;  basic  experiments  to 
test  the  validity  of  proposed  instability  mechanisms  and  provide 
vital  data  for  use  in  the  theoretical  models;  and  controlled  testing 
of  liquid  rocket  motors  with  parametric  variations  in  order  to 
isolate  and  explore  the  importance  of  the  controlling  factors  that 
influence  the  incidence  of  nonlinear  combustion  instability. 

In  the  study  of  nonlinear  combustion  instability  (i.e., 
high  frequency  instability  that  requires  an  energy  input  for 
initiation  which  may  be  internally  or  externally  supplied)  it  is 
important  to  remember  that  the  linear  mathematical  approaches  have 
often  proven  extremely  helpful  in  the  general  understanding  of  the 
phenomena  involved.  In  this  report  linear  treatments  will  be  used 
as  well  as  purely  nonlinear  theoretical  approaches  to  explore  the 
nature  of  the  instability  mechanisms.  A  most  important  result  of 
the  nonlinear  analyses  is  that  "continuity"  exists  between  linear 
theory  and  nonlinear  theory  if  the  feedback  mechanism  is  a  continuous 
function  of  the  amplitude  of  the  oscillation.  This  would  be  true 
for  combustion  processes  such  as  vaporization,  diffusion  and  chemical 
reaction,  but  would  not  apply  to  droplet  shattering.  These  processes 
are  classified  as  nonlinear  processes  under  the  proper  amplification 
conditions. 

In  this  research  one  of  the  prime  efforts  has  been  to 
closely  relate  the  theoretical  investigations  to  the  basic  experiments 
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and  rocket  motor  parametric  studies.  To  emphasize  this  aspect  of 
equal  importance,  this  report  has  been  written  with  a  somewhat  dif- 
ferent  format  than  in  the  past,  where  the  theoretical  and  basic 
experimental  studies  were  separately  covered  in  the  Appendices. 

As  is  the  case  in  any  status  report,  a  considerable  amount 
of  the  material  covered  is  in  the  state  of  active  investigation  with 
only  tentative  conclusions  or  observations  possible.  Every  effort 
has  been  made  to  include  some  mention  of  all  the  work  in  progress 
so  that  the  report  may  prove  of  maximum  benefit  to  those  readers 
engrossed  in  similar  research  or  confronted  with  development  problems 
involving  combustion  instability. 

Since  the  previous  yearly  progress  report  (June  1,  1963) 
a  number  of  technical  reports  have  been  distributed,  which  cover 
certain  aspects  of  the  research  in  far  greater  detail  than  is  possible 
in  this  report.  They  include  the  following,  listed  chronologically: 

"An  Optical  Method  for  Observing  Breakup  and  Vaporization 
of  Liquid  Jets,"  W.  R.  Seebaugh  and  D.  H.  Lee,  Aeronautical 
Engineering  Report  No.  647,  June  1963. 

"A  Theoretical  Study  of  Droplet  Burning:  Transients  and 
Periodic  Solutions,"  W.  C.  Strahle,  Aeronautical  Engineering  Report 
No.  671,  December  1963. 

"A  Theoretical  Study  of  Nonlinear  Combustion  Xnstabili  ty: 
Longitudinal  Mode,"  W.  A.  Sirignano,  Department  of  Aerospace  and 
Mechanical  Sciences  Report  No.  677,  March  1964. 

"Evaluation  of  a  Light  Scattering  Technique  For  Determining 
the  Spray  Characteristics  of  Impinging  Liquid  Jets,"  H.  R.  Bredfeldt, 
Department  of  Aerospace  and  Mechanical  Sciences  Report  No.  648, 

March  1964. 

Current  progress  in  each  of  these  general  subject  areas 
is  presented  in  this  report  together  with  many  other  aspects  of  the 
research  in  combustion  instability. 
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In  addition  to  the  Princeton  University  technical  reports, 
the  following  articles  have  appeared  or  have  been  accepted  for 
publication  during  the  past  year. 

’•Analytical  Investigation  of  Several  Mechanisms  of 
Combustion  Instability,"  W.  C.  Strahle  and  L. 

Crocco,  Bulletin  of  the  Fifth  Liquid  Propulsion 
Symposium,  13-15  November  1963,  Chemical  Propulsion 
Information  Agency. 

"A  Shock  Wave  Model  of  Unstable  Rocket  Combustors," 

W.  A.  Sirignano  and  L.  Crocco,  AIAA  Preprint  No. 

64-143  Solid  Propellant  Rocket  Conference,  Palo 
Alto,  California,  January  (accepted  for  AIAA 
Journal  publication  expected  this  summer). 

"Velocity  Effects  in  Transverse  Mode  Liquid  Propel¬ 
lant  Rocket  Combustion  Instability,"  F.  H.  Reardon, 

L.  Crocco  and  D.  T.  Harrje  (accepted  for  AIAA  Journal 
publication) . 

"Problems  in  Liquid  Propellant  Instability,"  D.  T. 

Harrje,  Proceedings  of  the  Fourth  Meeting  of  the 
Technical  Panel  on  Solid  Propellant  Combustion 
Instability,  Applied  Physics  Laboratory  TG  371-7, 

April  1964. 

"Periodic  Solutions  to  a  Convective  Droplet  Burning 
Problem:  The  Stagnation  Point,"  W.  C.  Strahle 
(Accepted  by  the  International  Combustion  Symposium 
to  be  held  in  Great  Britain,  August  1964). 

"Theoretical  Studies  on  Liquid  Propellant  Rocket 
Instability,"  L.  Crocco,  (Invited  paper  for  the 
International  Combustion  Symposium). 

Additional  background  and  history  of  the  research  on 
combustion  instability  at  Princeton  University  may  be  found  in 
References  1  through  5  and  in  previous  technical  reports  which  will 
be  referenced  at  appropriate  points  in  the  text.  Background  in¬ 
formation  on  the  specific  topic  of  nonlinear  combustion  instability 
may  be  found  in  References  1-3  and  6. 
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III.  NONLINEAR  THEORETICAL  APPROACHES- LONGITUDINAL  AND  TRANSVERSE  MODES 
A.  GENERAL 

Theoretical  studies  of  the  nonlinear  oscillations  which  com¬ 
monly  occur  in  rocket  combustion  chambers  and  exhaust  nozzles  have  been 
performed  at  Princeton  over  the  past  few  years.  The  concern  has  been 
with  the  ■'high- frequency"  type  of  oscillations  wherein  energy  is  fed- 
back  to  the  oscillation  by  the  combustion  process  and  the  nozzle  acts 
as  a  damping  device.  Shock  wave  dissipation,  if  present,  provides  an¬ 
other  means  of  damping  the  oscillation. 

Both  the  longitudinal  and  transverse  modes  have  been  analvzed. 
A  shock  wave  model  of  the  longitudinal  mode  with  no  time-lag  effects 
was  considered  and  found  to  have  application  to  premixed  gas  rockets. 
Preliminary  results  were  reported  in  Ref.  1,  while  final  results  are 
presented  in  Ref.  6  and  7.  A  second  model  of  the  longitudinal  mode 
with  time-lag  effects  but  no  shock  waves  has  been  investigated  and 
the  results  were  presented  in  Ref.  6.  Further  results  of  this  second 
model  are  presented  in  Part  B  of  this  section.  A  third  longitudinal¬ 
mode  model  has  been  analyzed  and  is  oresented  in  Part  C  of  this  sec¬ 
tion.  This  model  considers  the  presence  of  shock  waves  and  a  small 
time-lag  effect.  Analyses  of  the  transverse  mode  in  annular  chamber 
were  presented  in  Ref.  2  and  3.  The  analysis  of  transverse  oscillations 
with  time-lag  effects  in  a  full  circular  chamber  is  presently  being 
performed.  A  portion  of  that  work  is  presented  in  Part  D  of  this  sec¬ 
tion. 

Theoretical  work  on  nonlinear  combustion  instability  has  also 

g 

been  performed  by  Priem  and  Guentert  and  Chinitz,  Bursteln,  and 
9  10 

Agosta  ’  .  In  their  approaches  to  the  problem  ,an  initial  disturbance 

was  applied  and  the  resulting  flow  behavior  was  determined.  That  is, 
it  was  determined  (typically  by  numerical  integration  of  the  equations) 
whether  the  disturbance  grew  or  decayed  with  time.  Princeton's  work 
has  been  different  in  that  periodic  solutions  have  been  obtained,  chief¬ 
ly  by  analytical  means. 

In  all  three  longitudinal  models  which  were  investigated, 
the  assumptions  of  concentrated  combustion  zone  at  the  injector  end  and 
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of  short  nozzle  were  made.  The  Important  distinguishing  factor  of  the 
three  models  is  the  characteristic  time  of  the  combustion  process  as 
compared  tc  the  wave  travel  time  in  the  chamber.  The  characteristic 
combustion  time  in  the  first  model  is  negligible  compared  to  the  wave 
travel  time  while  in  the  third  model,  the  combustion  time  is  small  but 
not  negligible  compared  to  the  wave  travel  time.  Clearly,  one  may  con* 
sider  the  first  model  as  special  case  of  the  third  in  the  limit  as  the 
combustion  time  goes  to  zero.  The  small  combustion  time  implies  that 
the  phase  between  energy  addition  and  pressure  is  small  and  goes  to 
zero  as  the  combustion  time  goes  to  zero.  The  second  model  Involves 
a  combustion  time  which  is  of  the  same  order  of  magnitude  as  the  wave 
travel  time.  This  means  that  the  phase  between  energy  addition  and 
pressure  is  of  the  order  of  the  period  of  the  oscillation. 

Specifically,  the  characteristic  combustion  time  in  the  sec¬ 
ond  model  is  Introduced  by  means  of  the  Crocco  time-lag  postulate. 

The  validity  of  this  postulate  has  been  experimentally  demonstrated  for 
linear  oscillations  but,  as  of  this  time,  not  for  nonlinear  oscilla¬ 
tions.  The  combustion  time  in  the  third  model  is  introduced  by  means 
of  a  convenient  postulate  which  has  no  experimental  support;  the  gas 
velocity  at  the  time  t  at  the  combustion  zone  is  assumed  to  be  given 
as  a  power  series  in  both  the  thermodynamic  conditions  at  time  t  and 
the  thermodynamic  conditiorPat  the  time  t  Z*  has  a  small  value, 

and  may  be  either  positive  (lag)  or  negative  (lead).  This  postulate  is 
made  so  that  the  effects  of  a  small  phase  between  energy  addition  and 
pressure  upon  the  oscillation  could  be  determined.  For  the  purpose  of 
determining  these  effects  qualitatively,  the  exact  functional  form  of 
the  combustion  zone  energy  feedback  relation  is  not  important  provided 
that  the  phasing  properties  are  contained  in  this  relation.  For  this 
reason,  it  is  believed  that  the  use  of  this  relation  is  justified. 

With  the  first  model,  solutions  were  found  which  consisted 
of  shock  discontinuities  followed  by  exponential  decays  in  pressure  and 
gas  velocity  as  shown  in  Figure  III— 1.  Nonlinear  oscillations  were 
only  possible  in  the  region  of  "linearly"  unstable  operation  with  the 
amplitude  increasing  with  distance  from  the  linear  neutral  stability 


line.  The  amplitude  vent  to  zero  at  this  neutral  line  Implying  that 
the  nonlinear  and  linear  results  were  in  agreement.  The  results  of  the 
third  model  are  qualitatively  the  same  with  respect  to  wave  form,  regions 
of  instability,  amplitude  vs.  distance  from  neutral  line,  and  continuity 
between  linear  and  nonlinear  results. 

With  the  second  model,  periodic  solutions  without  shock  waves 
were  obtained.  These  finite- amplitude  solutions  could  be  either  in¬ 
side  or  outside  of  the  linearly-unstable  region.  Figure  III- 2  repre¬ 
sents  a  three-dimensional  plot  of  time-lag  (£*)  vs.  interaction  index 
(n)  vs.  amplitude  parameter  (C)  for  periodic  oscillations.  The  amp¬ 
litude  parameter  is  seen  to  be  proportional  to  the  square  root  of  the 
displacement  from  the  neutral  stability  line  in  the  Zf  ,  n  plane. 

The  periodic  solutions  without  shock  waves  outside  of  the  linearly  un¬ 
stable  region  (or  in  the  linearly  stable  region)  were  found  only  at  the 
lower  values  of  the  interaction  index,  n  (approximately  n<.2.0). 

Note  that  as  n  decreases  the  phase  between  energy  addition  and  pres¬ 
sure  decreases,  and  goes  to  zero,  as  n  goes  to  its  minimum  value  on 
the  neutral  stability  line.  This  type  of  periodic  solution  was  shown 
to  be  unstable  indicating  the  possibility  of  "triggering"  action  for 
lower  values  of  the  phase.  That  is,  a  periodic  disturbance  with  a  great¬ 
er  amplitude  parameter  than  some  critical  value  C  for  given  values  of  T 
and  n  would  grow  in  amplitude  resulting  in  a  fully-established  oscil¬ 
lation  (presumably  with  shock  waves)  while  a  disturbance  with  smaller 
than  the  critical  magnitude  would  decay  to  zero  amplitude.  (For  de¬ 
tails,  see  Ref.  6).  The  periodic  solutions  without  shock  waves  and  in 
the  linearly-unstable  region  were  found  only  at  larger  values  of  the 
interaction  index  (approximately  n>2.0).  At  these  larger  values  of 
n,  the  phase  between  energy  addition  and  pressure  was  larger.  These 
were  shown  to  be  stable  periodic  solutions  and  their  wave  forms  were 
meaningful  and  could  be  calculated.  The  numerical  results  for  the  funda¬ 
mental  mode  were  presented  in  Ref.  6  while  the  results  for  the  harmonics 
are  presented  in  Part  B  of  this  section. 

The  numerical  calculations  presented  in  Ref.  6  indicated  that 
the  amplitude  of  the  disturbance  necessary  to  trigger  an  oscillation 
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became  Infinite  as  phase  went  to  zero.  This  meant  "triggering"  action 
Is  impossible  in  that  limit  which  agrees  with  the  results  of  the  third 
model  with  small  phase  where  finite  amplitude  oscillations  were  found 
only  in  regions  which  were  linearly  unstable  and  never  in  linearly 
stable  regions. 

Further  agreement  is  obtained  between  the  second  and  third 

models  in  that  the  mean  pressure  level  under  oscillation  is  found  to 

be  higher  (lower)  than  the  steady-state  value  if  energy  addition  leads 

(lags)  the  pressure.  Finally,  in  both  models,  the  frequency  is  higher 

(lower)than  the  resonant  value  if  energy  addition  leads(lags)  pressure. 

In  the  analysis  of  the  transverse  oscillations  presently  being 

undertaken,  the  approach  is  similar  to  that  used  by  Maslen  and  Moore** 

in  that  a  perturbation  technique  is  employed  and  solutions  without  shock 

★ 

waves  are  sought.  However,  the  present  problem  is  considerably  more 
complicated  in  that  the  pressure  of  combustion  and  the  effect  of  a 
nozzle  are  considered  here  whereas  Maslen  and  Moore  considered  the  case 
of  no  mean  flow,  no  energy  or  mass  addition,  and  a  solid  wall  boundary 
condition. 

The  Crocco  time-lag  postulate  is  employed  in  order  to  repre¬ 
sent  the  feedback  of  energy  to  the  oscillation.  In  the  initial  attempt, 
the  assumption  of  concentrated  combustion  zone  at  the  injector  end  is 
being  made  as  it  was  in  the  studies  of  the  longitudinal  oscillations. 

The  short  nozzle  assumption  used  in  the  longitudinal  work  has  no  useful 
counterpart  in  the  transverse  case  since,  while  it  is  reasonable  to 
think  of  a  chamber  length  much  longer  than  a  nozzle  length,  the  chamber 
and  nozzle  diameters  are  of  the  same  order  of  magnitude.  For  this  rea¬ 
son,  quasi-steady  results  cannot  provide  a  proper  boundary  condition 

at  the  nozzle  entrance  and  it  is  necessary  to  analyze  the  unsteady  noz- 

12 

zle  flow.  Crocco  has  done  this  for  the  linear  case.  Following  similar 
lines,  the  work  is  being  extended  to  the  nonlinear  case  and  that  effort 
is  reported  in  Part  D  of  this  Section. 

*  Experimental  evidence  Indicates  that,  typically,  shock  waves  do 
not  occu^  in  the  transverse  mode  in  a  full  chamber. 


B.  FURTHER  CALCULATIONS  BASED  ON  THE  NONLINEAR  LONGITUDINAL  MODEL 
WITH  TIME-LAG  EFFECTS 

Calculations  based  on  the  second  nonlinear  longitudinal  model 
presented  in  Ref.  6  have  been  completed.  This  particular  model  consid¬ 
ered  the  case  of  one-dimensional,  unsteady  flow  in  a  combustion  chamber. 
The  assumptions  of  a  concentrated  combustion  zone  at  the  injector  end 
and  a  short  nozzle  length  were  made.  A  phase  between  energy  addition 
and  pressure  oscillation  was  introduced  by  employment  of  the  Crocco 
time-lag  postulate.  In  particular,  periodic  solutions  without  shock 
waves  were  sought.  Two  types  of  solutions  were  found:  stable  periodic 
solutions  and  unstable  periodic  solutions.  The  latter  type  indicates 
the  possibility  of  "triggering"  an  oscillation  by  the  application  of  a 
finite  amplitude  disturbance. 

Only  the  numerical  results  for  the  fundamental  mode  with  a 

mean  flow  Mach  number  equal  to  .2  were  presented  in  that  report.  The 

results  for  other  Mach  numbers,  .1  and  .3,  and  for  other  modes,  second 

* 

and  third  harmonics,  arc  presented  in  this  report.  Figure  III- 3  shows 

the  results  for  the  fundamental  mode  with  Mach  number  u  *  .2.  In 

o 

Figure  III-4  and  III-5,  the  results  for  the  cases  of  the  fundamental 

mode  with  u  *  .1,  .3,  respectively,  are  shown.  In  Figures  III-6  and 
o 

III- 7,  the  results  for  the  cases  of  the  second  and  third  harmonics,  re¬ 
spectively,  with  u^  *  .2  are  shown. 

In  the  analysis  of  Ref.  6,  the  solution  was  found  in  the  form 

of  a  perturbation  series  in  an  amplitude  parameter  £■  .  The  solution 

3 

was  approximated  to  an  accuracy  of  0  (  £  ).  To  first  order,  the  re¬ 
sults  of  the  nonlinear  analysis  were  identical  with  the  results  of  Crocco 
linear  analysis.  To  this  lowest  order,  the  wave  form  consisted  of  a 
single  Fourier  component  with  amplitude  £  .  This  component  could  be 
either  the  fundamental  or  any  of  the  overtones.  The  neutral  stability 
line  for  zero- amplitude  oscillations  was  identical  to  that  found  by 

The  calculations,  although  arithmetical,  were  quite  tedious  and 
were  performed  with  the  aid  of  a  7090  computer.  The  results 
were  plotted  on  a  cathode  ray  oscillioscope  and  photographs  of 
the  screen  were  taken.  These  photographs  are  presented  in  the 
figures . 
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Crocco  as  shown  in  (a),  (b),  and  (c)  of  Figures  III- 3  through  III-7. 

In  (a),  the  frequency  in  characteristic  coordinates  is  plotted 
versus  the  interaction  index  n  along  the  neutral  line,  (b)  shows  a 
plot  of  vs,  n  where  GJ  is  the  frequency  in  time,  space  co¬ 

ordinates  and  £  is  the  sensitive  time-lag.  ^  vs.  n  is  plotted  in 
(c). 

It  can  be  seen  from  (a)  that  there  is  a  band  of  frequencies 
possible  rather  than  discrete  values.  This  Is  due  to  the  time-lag 
which  introduces  another  characteristic  time  besides  the  wave  travel 
time.  The  natural  frequency  for  the  particular  mode  occurs  at  the  min¬ 
imum  value  of  n  for  neutral  oscillation.  Both  above  and  belov- 
resonant  frequencies  are  possible.  The  region  inside  the  "parabola¬ 
like"  curve  is  the  instability  region  while  the  outside  of  the  curve 
is  the  stability  region.  The  figures  show  that  the  larger  the  Mach 
number  the  wider  the  instability  regions.  (b)  is  identical  in  all 
five  figures  indicating  that  the  ratio  of  the  time-lag  to  the  period 
of  oscillation  is  independent  of  the  mode  of  oscillation  and  of 

Mach  number  •  *£  vs.  n  for  neutral  oscillations  is  seen  to  be  double¬ 

valued  from  (c).  The  minimum  value  of  n  corresponds  to  a  resonant 
oscillation.  The  branch  of  the  curve  with  high  (H)  values  of  £  cor¬ 
responds  to  be  low- resonant  oscillations  while  the  branch  with  low  (L) 
values  of  corresponds  to  above- resonant  oscillations.  The  region 

of  instability  is  inside  the  "parabola- like"  curve.  This  region  widens 
as  Mach  number  increases.  The  instability  region  moves  to  lower  values 
of  £  as  the  mode  of  oscillation  becomes  higher  in  such  a  manner  that 
is  unchanged. 

There  are  two  effects  upon  the  wave  form  which  appear  to 

★★★ 

second  order;  the  mean  pressure  level  changes  from  the  steady- state 

*  It  was  convenient  to  analyze  the  wave  phenomenon  by  trans¬ 
forming  the  equations  of  motion  to  a  coordinate  system  formed 
by  the  characteristics  of  the  hyperbolic  system. 

**  Actually,  if  Mach  number  squared  effects  are  considered,  there 
is  a  slight  dependence.  See  Section  VI. 

This  implies  a  change  in  the  mean  thrust  level  from  the  steady- 
state  value  dependent  upon  the  particular  wave  form  present  at 
the  nozzle  entrance.  The  significance  of  this  change  will 
be  determined  by  further  analysis. 
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value  and  a  second  harmonic  Is  added  with  a  phase  with  respect  to  the 

first.  ("Second  harmonic"  means  the  second  Fourier  component.  If  the 

first  Fourier  component  were  the  first  overtone  of  the  chamber,  the 

second  harmonic  would  be  the  third  overtone.)  The  amplitude  of  this 

★  2 

second  harmonic  for  one  of  the  Riemann  invariants  is  given  by  £  A 

with  the  phase  given  by  2  ©  .  A  vs  n  is  plotted  in  (d)  of  Figures 

2 

III-3  through  III-7  and  ©  vs.  n  is  plotted  in  (e).  £  is  the 
change  in  the  mean  value  of  the  Riemann  Invariant  and  is  plotted  in 
(f)  of  these  figures.  All  three  types  of  curves  are  seen  to  be  double¬ 
valued  with  an  above- re  sonant  branch  (L)  and  a  be  low*- re  sonant  branch 
(H). 

A  and  are  seen  to  tend  towards  infinity  as  the  resonant 

point  (minimum  value  of  n)  is  approached.  The  mean  values  of  A  and  C 
**  r 
are  higher  for  the  higher  modes.  Note  that  in  the  above- resonant  case 

(L) ,  the  sign  of  A  could  be  changed  (i.e.,  made  positive),  if  cor¬ 
respondingly,  the  phase  2  ©  was  translated  by  the  value  rr  .  In  this 
way  it  would  be  easier  to  see  that  the  effect  of  the  addition  of  the 
second  harmonic  is  not  very  different  for  above  and  be  low- resonant  cases. 
As  shown  by  (f),  the  mean  pressure  level  is  higher  than  steady- state  for 
the  above- resonant  case  and  lower  than  steady- state  for  the  below- resonant 
case.  There  is  no  simple  pattern  in  the  effect  of  Mach  number  upon  A 
and  ©  ,  however,  the  absolute  value  of  C  decreases  as  Mach  number 
increases . 

The  third  order  effect  upon  the  wave  form  is  that  a  third 

harmonic  (third  Fourier  component)  is  added  with  a  phase  with  respect 

3 

to  the  first.  The  amplitude  of  this  harmonic  is  C  B  and  the  phase  is 
3  -A  .  B  vs.  n  is  plotted  in  (g)  of  Figures  III-3  through  I II- 7  and 
A  vs.  n  is  plotted  in  (h). 

(g)  shows  that  B  tends  towards  infinity  as  the  resonant  point 
is  approached.  The  mean  value  of  the  absolute  magnitude  of  B  is  seen 

*  The  Riemann  invariants  are  certain  combinations  of  flow  properties 
which  are  constant  along  a  characteristic.  The  speed  of  sound 
perturbation  at  the  injector  end  is  directly  proportional  to  the 
Riemann  Invariant.  See  Ref.  6. 

**  This  means  higher  in  absolute  value. 
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to  be  higher  for  the  higher  modes  while  A  remains  essentially  the  same 
with  changes  In  mode  and  Mach  number.  No  simple  pattern  is  seen  In  the 
behavior  of  B  upon  Mach  number.  Note  that  the  amplitude  factor  B  In 
the  above- resonant  case  (L)  may  have  its  sign  changed  (made  positive)  If, 
correspondingly,  A  was  translated  by  ^/3.  After  this  change  the 
phase  would  be  smaller  over  a  wider  range  of  n  In  both  above  and  below- 
resonant  cases. 

On  the  basis  of  the  calculations  above,  one  may  calculate  the 
approximate  value  of  2*  and  n  where  an  oscillation  of  a  given  amplitude 
occurs.  For  values  of  ^  and  n  along  the  neutral  line  of  (c),  the  amp¬ 
litude  Is  zero.  The  amplitude  parameter  £  Increases  with  the  square 

root  of  the  normal  displacement  from  the  neutral  line  as  schematically 

2 

shown  in  Figure  III- 2.  Specifically,  at  a  displacement  of  DC  from 
the  neutral  line,  the  amplitude  of  the  first  Fourier  component  for  neutral 
finite-amplitude  oscillation  is  €  .  D  vs.  n  is  plotted  in  (1)  of 

* 

Figures  III-3  through  III- 7 .  D  may  be  positive,  indicating  outward  dis¬ 
placement,  or  negative,  indicating  inward  displacement.  It  was  shown  in 
Ref.  6  that  D^O  corresponds  to  unstable  periodic  solutions  and  0 
implies  stable  periodic  solutions,  (i)  shows  that  D<  0  occurs  only  at 
n  values  greater  than  2.0.  According  to  the  experimental  results  of 
Ref.  6,  n  <  2  is  the  range  of  physical  interest  so  the  theory  predicts 
that  stable  periodic  solutions  without  shock  waves  are  not  likely  to  oc¬ 
cur  in  practical  situations.  This  conclusion  agrees  with  experimental 
findings  since  shocks  are  usually  observed  to  be  present  in  the  wave 
form  of  rockets  unstable  in  the  longitudinal  mode.  The  fact  that  unstable 
periodic  solutions  are  found  in  the  range  of  n-values  of  practical  in¬ 
terest  indicates  that  "triggeringM  action  is  possible  there.  This  agrees 
at  least  qualitatively  with  experimental  findings.  (See  Section  IV-B). 

The  values  of  D  become  more  positive  as  the  mode  becomes  higher  implying 
that  the  triggering  of  an  overtone  mode  requires  a  lower  amplitude  dis¬ 
turbance  than  the  triggering  a  fundamental  mode.  The  lowest  Fourier 

*  "Outward  displacement"  means  into  the  region  of  linear  stability 
while  "inward  displacement"  means  into  the  region  of  linear  in¬ 
stability. 
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component  of  the  disturbance  must  be  Identical  to  the  mode  being  triggered, 
however,  so  this  does  not  imply  an  overtone  is  more  likely  to  be  triggered 
than  the  fundamental.  There  is  again  no  clear  pattern  in  the  relation¬ 
ship  between  D  and  Mach  number,  Indicating  that  extropolation  from 
results  for  one  value  of  Mach  number  to  another  value  could  be  gravely 
in  error.  This  points  to  the  need  of  the  work  of  Section  VI,  since  the 
nonlinear  parameters  presented  here  depend  upon  V  and  n  values  along 
the  neutral  lines  which  in  turn  depend  upon  the  Mach  number.  Because 
of  the  sensitivity  of  these  nonlinear  parameters  to  Mach  number,  a  more 
accurate  analysis  of  the  Mach  number  effects  would  be  very  useful. 

The  calculation  of  the  wave  forms  would  only  be  meaningful 
for  stable  periodic  solutions  which  are  seen  to  occur  only  outside  of 
the  range  of  practical  Interest.  However,  in  case  propellants  with 
higher  interaction  Indexes  are  developed  or  in  order  to  compensate  for 
possible  error  in  the  prediction  of  the  D  »  0  point  due  to  a,  proxima- 
tions  in  the  theory,  the  wave  forms  have  been  calculated  and  the  pre¬ 
dictions  for  n  ■  4.0  are  presented  in  Figure  III-8.  The  pressure  vs. 
time  wave  forms  were  calculated  for  both  injector  and  nozzle  chamber 
locations,  for  both  the  fundamental  and  second  harmonic  modes,  and  for 
both  above- resonant  and  be  low- resonant  oscillations.  The  value  of  £  is 
.1  for  the  results  of  Figure  I II- 8 . 

In  (a)  of  that  figure,  we  have  above- resonant  oscillation  in 
the  fundamental  mode  at  the  nozzle  entrance.  Three  approximations  are 
shown:  P3  vs.  T,  which  is  a  first  order  approximation,  P4  vs.  T2,  the 
second  order  approximation,  and  P6  vs.  T4 ,  the  third  order  approximation. 
The  P's  are  the  approximations  to  the  pressure  and  the  T's  are  the  ap¬ 
proximations  to  the  time  which  appear  by  means  of  the  characteristic 
coordinate  transformation.  These  T  corrections  are  related  to  the 
wave  distortion  and  are  larger  than  the  P  corrections.  The  second  and 
third  order  approximations  are  essentially  identical  implying  that  the 
third  order  correction  to  the  wave  form  is  negligible.  The  important 

nonlinear  effects  are  the  sharpening  of  the  positive  pressure  peak  and 

★ 

the  increase  in  the  mean  pressure  level. 


*  The  steady-stace  pressure  is  1. 
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In  Figure  III-8  (b),  the  sene  cese  of  above- resonant,  funda¬ 
mental  mode  Is  plotted,  except  that  now  the  location  Is  the  Injector- 
end  of  the  chamber.  This  wave  form  (b)  Is  different  from  the  wave  form 
at  the  nozzle-end  of  the  chamber  (a)  only  In  that  the  amplitude  Is  low¬ 
er  due  to  the  phase  In  reflection  at  the  Injector-end  which  does  not  oc¬ 
cur  at  the  nozzle-end.  This  similarity  In  wave  form  Is  expected  due  to 
the  nature  of  the  wave  propagation. 

Figure  III-8  (c)  shows  the  case  of  fundamental-mode  oscilla¬ 
tion  at  the  nozzle-end.  It  differs  from  (a)  In  that  the  oscillation  Is 
now  be  low*  re  sonant  causing  a  sharpening  of  the  negative  pressure  peak 
and  a  decrease  In  the  mean  pressure  level. 

Figure  III-8  (d)  displays  the  results  for  the  second-harmonlc- 
mode,  above- re sonant  oscillation  at  the  nozzle  end.  The  comparison 
with  (a)  shows  that  the  positive  peak  Is  somewhat  sharper  and  the  mean 
pressure  Is  higher  with  the  overtone  than  with  the  fundamental .  In 
the  case  of  be  low- resonant  oscillations,  the  negative  peak  would  be 
sharper  and  the  mean  pressure  lower  with  the  overtone  than  with  the  funda¬ 
mental  . 

For  both  modes,  the  peaks  tend  to  sharpen  and  the  deviation 
of  the  mean  pressure  from  the  steady- state  value  Increases  as  n  de¬ 


creases. 
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C.  A  SHOCK  WAVE  MODEL  OF  UNSTABLE  ROCKET  COMBUSTORS  WITH  SMALL 
TIME-LAG  EFFECTS 

The  case  of  longitudinal  oscillations  with  shock  waves 
in  a  rocket  combustor  is  investigated  in  this  section.  The  charac¬ 
teristic  time  of  the  combustion  process  is  assumed  to  be  small  com¬ 
pared  to  the  wave  travel  time.  A  time-lag  or  phase  between 
pressure  and  energy  addition  introduces  this  characteristic  time  of 
the  combustion  process.  The  relation  between  this  work  and  other 
nonlinear  work  performed  at  Princeton  is  given  in  Section  III  A. 

The  combustion  zone  length  is  considered  negligibly  small 
compared  to  the  chamber  length  so  that  the  limiting  case  of  con¬ 
centrated  combustion  at  the  injector-end  is  considered.  The  nature 
of  the  combustion  process  determines  the  boundary  condition  on  the 
oscillation  at  this  end.  The  perturbation  in  gas  velocity  is  related 
to  the  perturbations  in  pressure  and  temperature.  Typically,  an 
increase  in  one  or  more  of  the  thermodynamic  variables  causes  an 
increase  in  the  velocity  of  the  burned  gas  emitted  from  the  combustion 
zone  since  the  burning  rate  increases.  Other  assumptions  are: 

1.  The  flow  is  one-dimensional. 

2.  The  chamber  cross-sectional  area  is  constant. 

3.  The  lergth  of  the  convergent  portion  of  the  nozzle  is 
negligibly  small  compared  to  the  chamber  length  so  that 
the  limiting  case  of  the  short  nozzle  is  considered. 

4.  A  shock  wave  moves  back  and  forth  the  length  of  the 
chamber  with  a  constant  period,  reflecting  alternately 
from  the  nozzle  and  the  combustion  ends.  Here  only 
one  shock  wave  is  considered,  so  only  the  fundamental 
mode  is  allowed. 

5.  Flow  is  homentropic  outside  of  the  combustion  zone  up 
to  and  including  second  order  in  the  wave  amplitude. 

This  allows  shock  waves  to  occur  but  no  entropy  waves 
are  allowed. 
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6.  The  chamber  gas  is  calorically  perfect. 


Undei  these  assumptions,  the  steady-state  flow  properties 
are  constant  throughout  the  chamber.  All  thermodynamic  variables 
are  nondimensionalized  with  respect  to  their  steady-state  values. 

The  gas  velocity  is  nondimensionalized  with  respect  to  the  steady- 
state  speed  of  sound,  space  dimension  with  respect  to  chamber  length 
and  tinv.  dimension  with  respect  to  chamber  length  divided  by  speed 
of  sound. 


The  coordinate  system  is  taken  in  such  a  manner  that  x  *  ( 
at  the  nozzle  entrance  location  and  x  *  1  at  the  combustion  zone 
location.  Note  that  this  means  that  the  steady-state  velocity  is 
negative  in  sign. 

The  pressure  and  temperature  at  the  injector  end  may  be 
related  to  the  speed  of  sound  by  means  of  the  laws  of  thermodynamics 
and  an  equation  of  state.  Therefore,  a  relationship  between  the  gas 
velocity  perturbation  and  the  speed  of  sound  perturbation  may,  in 
principle,  be  found.  Here,  it  is  assumed  that  the  time-lag  effect 
appears  in  such  a  manner  that  this  relationship  may  be  written  as 
follows:  r  /  v  /  1 

U(t)-u*  fe  co  K.CaftJ-aj)*  Kg  La  (t-?)  -%  )  +• 

+  S  r  (a  ^  (a  j 


+  higher  order  terms 

where  u  is  the  gas  velocity,  a  is  the  speed  of  sound,  t  is 
the  time,  ^  is  the  time-lag,  and  O) ,  &,  ,  k^  ,  and 

k^  are  positive  parameters  which  describe  the  combustion  process. 
Zero  subscripts  denote  steady-state  values.  Note  that  the  relation¬ 
ship  has  been  written  in  the  form  of  a  Taylor  series  where  the  third 
and  higher  order  terms  are  neglected. 

With  no  loss  of  generality,  we  may  consider  k^  +  k^  ■  1  . 
The  time-lag  is  small  and  is  considered  aF  a  firsi  order  perturbation 
quantity.  This  means  that 

a  U-V)-a-Ct)=  \(eL(t-Y)-a9)-(a.a)-at 


T 
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must  be  a  second  order  perturbation  quantity  since  it  is  the  difference 
between  two  first-order  perturbation  quantities  whose  arguments 
differ  by  a  first  order  perturbation  quantity.  Now,  the  relation¬ 
ship  which  provides  the  boundary  condition  at  the  chamber  end  may 
be  rewritten  in  the  following  manner: 


=  aj(a<t;-ad)+  u>Kt 

+•  higher  order  terms 

(1) 

The  well  known  compatibility  relations  may  be  obtained 
from  the  equations  of  unsteady  one-dimensional  motion  for  a  fluid. 
Under  our  assumptions,  the  relations  (accurate  to  second  order)  are: 


l 

do¬ 

+■  du.  =  O 

along 

dx  s  ( U»+-Ol)  dt 

r-i 

-tr 

do. 

-  do.  »  O 

along 

Jx  *  C lx- cl)  dt 

Let  primes  imply  perturbations  so  that 

u  =  u  +  u 1  ;  a  =  1  +  a ' 

o 

Furthermore,  let 


11  =  -*-+(, 
Xt  st  1 


')  ihr 


~  st 


Then,  these  differential  equations  are  rewritten  as  follows: 


£  (  —  CL '  +  U. 

St  V.  ar-1 


-sr(-h  a'  ■|c’)  '°'r>-§T  (rfa 

The  left-hand  sides  of  Equations  (2a)  and  (2b)  arc  of  first  order 


(2a) 

(2b) 


while  the  right-hand  sides  are  of  second  order.  Consider  u^  and 


a^  as  the  fir-t  order  approximations  to 
Furthermore,  consider  u_  and  a. 


u'  and  a'  ,  respectively, 
as  second  order  corrections  so 


that  Uj^  +  u2  and 


a^  +  a^  are  second  order  approximations  to  u' 


and  a'.  Then,  Equation  (2a)  yields 

fr(  ft-  +u')=o 


(3a) 
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and  Equation  (2b)  yields 


The  general  solutions  to  Equations  (3a)  and 


to  be  the  following 


(4a) 


(3b) 

a.  -  u  \ 

(4b) 

(3b)  are 

found 

(5a) 

(5b) 

where  the  functions  F  and  G  will  be  determined  later  from 
conditions  imposed  upon  the  solution.  These  relations  (3a)  and  (5b) 
are  used  to  substitute  into  the  inhomogeneous  parts  of  Equations 
(4a)  and  (4b). 

The  particular  solutions  to  Equations  (4a)  and  (4b)  are 
found  to  be  the  following 


6,q  i 

r-i 


_  u  _  _  ^A±L  j  x  G-(t»T=Tr^frtft»rhrJ 


(6a) 


(6b) 


where  the  subscript  t  indicates  differentiation  with  respect  to 
time  or,  equivalently,  in  this  case,  differentiation  with  respect 
to  the  argument  of  the  function.  Note  that  the  homogeneous  solutions 
to  (4a)  and  (4b)  may  be  included  in  F  and  G  ,  respectively,  with 
an  error  of  third  order  only  in  (6)  which  is  negligible  here.  The 
integrals  may  be  assigned  zero  as  a  lower  limit  since  any  constant 
of  integration  could  be  combined  with  F  or  G  . 

Equations  (5)  and  (6)  may  be  combined  to  yield  the  following 


* 


J 
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relations  : 


■H  U' 


Ff*  -  j5=;)  +  W  *  F.dr^^t.-jk; 

rt. 


■(■¥)6^)F< (t_  *.)  | 


&CWJ5 


(7a) 


r-i 


-  u 


G*(t+  -i— )  -  /  it1-)  i&Lt*  rtsykitt+ifc 

‘“S  '  W  /  O-u.)1 


(7b) 


It  is  well  known  that  the  function 


n 


a  +  u  is 


continuous  to  second-order  accuracy  through  a  shock  discontinuity 

2 

travelling  in  the  negative  x-direction  and  the  function  y-y  a  -  u 

is  also  continuous  to  second  order  across  a  shock  moving  in  the 

positive  x-direction.  (See  Reference  12b).  If  the  second-order  terms 

in  Equation  (7a)  are  neglected,  it  is  seen  that  the  function  F  is 

continuous  through  a  shock  moving  towards  the  nozzle  entrance  to 

at  least  first-order  accuracy.  The  first-order  approximation  for 

F  may  be  substituted  into  the  nonlinear  terms  of  Equation  (7a).  If 

F  is  continuous  across  th'  shock,  so  is  F  .  Furthermore,  even 

though  G  is  discontinuous  through  this  type  of  shock,  its  integral 

is  continuous.  Therefore,  the  nonlinear  terms  in  Equation  (7a)  are 

2a ' 

continuous  across  this  shock.  Since  y  +  u'  is  also  continuous 
to  second-order  accuracy,  this  means  that  F  is  continuous  to 
second-order  accuracy  across  a  shock  moving  in  the  negative  x-direction. 
By  a  similar  argument,  Equation  (7b)  would  show  that  G  is  continuous 
to  second-order  accuracy  through  a  shock  travelling  in  the  positive 
x-direction. 

Multivaluedness  of  the  solution  occurs  due  to  the  presence 
of  the  shock  waves.  This  nonlinear  effect  is  caused,  for  example, 


by  the  intersection  of  some  of  the  lines  [t  -  -r—  ■ 

l+u0 

with  a  shock  moving  in  the  positive  x-direction.  F  (t 


constant ] 

_ x_ 

1+u, 


)  has 


different  values  depending  upon  from  which  side  the  shock  is  ap¬ 
proached.  This  is  solely  a  nonlinear  effect  because  the  deviation 
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of  the  shock  path  from  a  [t  -  777-  *  constant]  line  is  a  nonlinear 

i_Hj°  x 

effect.  Of  course,  the  same  phenomenon  occurs  with  G  (t  +  q -  ) 

near  shock  waves  travelling  in  the  negative  x-direction.  However, 
once  the  shock  path  is  known,  the  correct  value  for  each  x,t  point 
is  easily  discerned. 

The  short  nozzle  boundary  condition  is  used  to  obtain  a 
relationship  between  the  functions  F  and  G  .  The  Mach  number 
at  the  entrance  of  the  nozzle  (x  ■  0)  is  set  equal  to  a  constant 
with  time  since  the  wave  reflection  is  a  quasi -steady  phenomenon  in 
the  short  nozzle  case.  Therefore, 

u/(4tj  =U.  a'  (Ot) 

Equations  (7a)  and  (7b)  are  solved  for  u'  and  a'  in  terms  of 
F  and  G.  Then,  the  results  are  substituted  into  Equation  (8)  to 
yield  the  following  nonlinear  relation  between  F(t)  and  G(t): 


f  n)  ^Rt)*&it)J  -ijE  y  ” 

-  0  /.1  *  (9) 

where  the  definition  if  *  -  u  has  been  made.  (Note  that  1/  is 

1  o 

a  positive  number.) 

If  the  nonlinear  terms  in  Equation  (9)  are  neglected,  we 

have 

(ltVJP(t)+  G-(t)  =  0 


This  means  that  the  first  order  approximation  to  the  relationship 
between  F  and  G  is: 

Flt)  =  T77  e£t> 

If  this  relation  is  used  to  substitute  for  F(t)  in  the  nonlinear 
terms  of  Equation  (9),  the  error  would  be  of  third  order  which  is 
negligible  in  this  analysis.  This  permits  the  determination  of  F 
as  a  nonlinear  function  of  G  .  The  result  follows: 


Fit) 


Or  It)  ♦ 


_  i-n>  iuta 


'i&tOJt 


J 


(10) 


Now,  Equation  (10)  is  employed  to  substitute  for  F(t)  in 
Equation  (7a)  and  (7b)  with  the  results: 


4  * 


_l£.  Qr  )  1-  Jjd  (±=jrf  L&fLz 

*V-  l-HJL,  4  \l4tfy 


JJ_  i-v.  G-  (t  -  wr.  )  f  f 

t  ttr 

_  ±jr  jut,  GcE)Jfl 

l+V  l-u%  J  J 


s-ct;je 


(lla) 


_  LL 


-  G(t+”*)“J±i  (t  ■*»  i-u«) 

4  “ 

C*)J£ 


j- 1  i  +  uo  i-ir  Crt  (i+JL)  f 

-  i  “ 7  J 


(lib) 


Another  condition,  on  the  function  G  appears  from  the 
combustion  zone  boundary  condition  (1).  a'  and  u'  at  the  com¬ 
bustion  zone  are  given  by  setting  x  =  1  in  Equations  (11)  and 
solving  the  simultaneous  system.  Noting  that  ?  is  a  perturbation 
quantity,  we  find  the  following  first  order  approximation  from  this 
procedure : 


Or  (t  ♦  _!_  )  -  ItMf  Gr  Ct  -  -L-  ) 

j-tt#  1 4>  tr  i-wu  1 


(12) 
1 


It  is  readily  seen  from  Equation  (12)  that  if  co  >  1 ,  G  (t  +  j_u  ) 

1  ^ 

>  G  (t  -  )•  In  this  case,  a  first  order  analysis  shows  an 

exponential  growth  of  a  small  perturbation  (unstable).  If  CO  =  1 , 

then  G  (t  +  )  -  G  (t  -  )  and  the  system  is  neutrally 

1_uo  l+uo  1  1 

stable.  In  the  case,  uJ  <  1  ,  G  (t  +  )  <  C  (t  -  1+u^  ) 

and  a  small  perturbation  decays  exponentially  (stable). 

With  no  loss  of  accuracy,  Equation  (12)  is  used  to  substitute 

for  G  (t  +  — —  )  in  the  nonlinear  terms  which  result  wl  jn  Equations 

v  l-u0' 

(1)  and  (11)  are  combined.  This  operation  yields  the  following: 
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ZTT)1  J  <** 
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AdL  A 

G 
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•  4  U.# 

t 

J 

G (l)JZ 


I  -  o)V 


u  u. 
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where  the  definitions  have  been  made 


(13) 


A  m  1  +  1  ~  ^ 

i  4 


K 


i  -uro-  i  "  i+V 

Now,  we  will  seek  periodic  solutions  for  the  flow  properties 

The  period  of  the  oscillation  for  the  fundamental  mode  is  the  time 
required  for  a  shock  wave  to  travel  once  back  and  forth  in  the  chamber 
The  shock  follows  the  path  ABC  shown  in  Figure  1II-9.  The  flow 
properties  in  Region  III  are  identical  to  those  in  Region  I  due  to 
periodicity.  The  period  is  the  time  during  which  the  shock  travels 
the  path  ABC.  This  time  is  represented  by  T  *  T, 

If  V  is  the  shock  velocity,  we  have 


ABC 


tab  +  TBC* 


jr-r-v  •  ph 


Ae 


Vac  dt 


(14) 


The  shock  velocities  to  first  order  accuracy  are  given  by  the  well 
known  relations  (Reference  12b). 


=  Uj  f  1  + 


u. 

It 


4-  CL 


-L 

2 


(13) 


In  this  representation,  conditions  at  the  shock  AB  have  been  set 
equal  to  conditions  at  the  shock  CD.  The  flow  properties  in  (15) 
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are  to  be  evaluated  at  the  shock  with  the  particular  side  of  the 
shock  given  by  the  Roman  numeral  subscript.  (See  Figure  III-9). 

If  Equations  (5),  (14),  and  (15)  are  combined,  we  may  solve 
for  T  in  an  approximate  manner.  If  Tq  is  the  zero  order  approxi¬ 
mation  and  Tq  +  is  the  first  order  approximation,  we  obtain 


t  *  0  at 

the  point 

A  of 

o 

±zJL  f  K 

,1±-UI  4- 

l- 

»  L 

1-  u# 

1  *■  u.0 

.  r  k 

.i  r 

L  (i*Vg 

LJ  L 

_  r+i 


Since  the  flow  properties  are  T-periodic,  so  is 
2a 

j— -  -  u'  ,  and,  therefore,  as  shown  by  Equation  (lib),  G  is 
T-periodic  to  first  order.  On  the  contrary,  however,  the  periodic 
condition  implies  that  to  second  order  (as  is  shown  by  means  of 
Equation  (lib)) 


G  (t  ♦  T  )  «  Gr  If  +  *  ) 


¥ 


K  It  G  (t 
I  -  u.  t 


t  ♦  T  -  rt 


7*-) 


Introducing  no  greater  error  to  the  analysis,  it  is  consistent  to 
rewrite  this  as  follows  » 


G(/>-»  T)  *  G(/>)  +  3-f  K  j+*o  Gf(f)  J 

where  ^  replaces  the  fornu  rgument  [t  t  x / ( 1  u^)).  Since 

to  second  order  accuracy*  G  ( ^ +  TQ)^  G  (^  +  T)  +  (Tq  -  T)  X 
G  ^  +  T)  G  (^  +  T)  -  Gytf  (ytf+  T),  Equation  (16)  and  (17) 

yield 

*For  positive  T  ,  this  expansion  fails  very  close  to  the  end¬ 
point  0  of  the  range  0  <  T  due  to  the  shock  discontinuity. 
For  negative  ,  failure  occurs  near  T  for  the  same  reason. 
However,  the  approach  is  still  correct  since  the  relation  gives 
the  asymptotic  behavior  everywhere  in  the  range  as  0. 


*1 
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G  )  =  Gjp  +  _£dL  *  f  *  *  &(1)<U 

*  *#  ^  u#  “« 

-  T  G>  (/•) 

'  r  '  (18) 

Setting  y#  equal  to  the  argument  (t  -  1/(1  +  uq))  in  Equation 
(13),  that  equation  may  be  combined  with  (16)  and  (18)  to  yield 
the  following 

0  •  (A-.)  &(/>,*  JV j.  ^  ]  <?  </>  [  &v>- 

-i 

+  [&</-*>-  fry,]  (19) 

Note  that  (A-l)  appears  in  the  first  term  multiplied  by  a  first 

order  quantity.  Since  all  th-'  other  terms  in  (19)  are  of  second 

order,  (A-l)  must  be  of  firsc  order.  Therefore,  A  may  be  replaced 

by  unity  and  U>\f  by  if  in  the  nonlinear  terms  as  has  already  been 

done  in  (19).  This  also  means  that  G  can  only  be  determined  to 

first  order  accuracy  from  the  relationship  since  the  second  order 

correction  now  does  not  appear  until  third  order.  Furthermore, 

since  G  is  T  -periodic  to  first  order  accuracy,  the  integrals 
o 

which  appeared  in  the  nonlinear  terms  canceled  each  other. 

Since  is  small,  the  approximation*  may  be  made  by 

means  of  a  Taylor  expansion  that 

Gr(yO  -  V)  -  Cr(f)  =  -  V  Gfif) 

and  this  may  be  substituted  into  (19)  to  obtain  the  following  result 


i  .  r&»A6* 

■  Jf  C-& 

where  certain  definitions  have  been  made 


(20) 


*Note  that  near  one  of  the  endpoints  (0  or  T,  depending  on  the  sign 
of  *£  ) ,  this  expansion  is  not  correct  for  nonzero  due  to  the 
discontinuity  in  G(  f>)  ;  however,  the  correct  asymptotic  behavior 
as  Z*  - — »  0  is  given  even  near  these  points. 
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As  explained  in  References  6  and  7,  Equation  (20)  has  a  physically 
reasonable  solution  only  for  C  =  0  .  This  also  provides  a  condition 
on  the  solution  which  is  used  to  determine  the  constant  of  integration. 


-  r  -  A  G- 

)  +  2^p  =  0 


This  has  the  solution 


which  setting  =  0  can  be  shown  to  be  identical  to  that 

References  6  and  7,  in  the  case  where  no  time-lag  effect  was  considered. 
In  the  limit  as  A  — *0  (or,  in  other  words,  —♦0),  (21)  becomes 
a  linear  relation  indicating  a  sawtooth  waveform. 

Equations  (5),  (10),  and  (21)  may  be  combined  to  give  a 
first  order  approximation  to  the  wave  form  in  Region  I.  If  the 
isentropic  condition  is  used  to  relate  the  pressure  perturbation  to 


the  speed  of  sound  perturbation, 


(22a) 
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Similarly,  Equations  (5)  and  (21)  may  be  combined  to  give  the  wave- 

form  in  Region  II  if  note  is  taken  that  to  first  order  G  (t  +  rj —  ) 

2  x  i_uo 

=  G  (t  -  ^ -"—2  +  y-^  )•  Again,  using  the  isentropic  relationship, 

we  obtain  the  first  order  approximation 

p-Po~Hr [«*>  ^.K^.]^,^(22b) 


The  waveform  consists  of  a  shock  discontinuity  followed  by  an  ex¬ 
ponential  decay  (  is  a  small  number  so  that  [1  -  X  / r ] 

is  positive).  As  A  — —»0  ,  the  waveform  goes  to  a  sawtooth.  (21) 
and  (5)  can  be  shown  to  imply  that  the  pressure  jump  across  the  shock 
is  constant  as  the  shock  travels  the  length  of  the  chamber.  The 
strength  does  change  in  reflection,  however,  and  is  greater  for  the 
shock  travelling  in  the  negative  x-direction  than  for  the  shock 
travelling  in  the  positive  x-direction. 

When  ^  is  positive,  (1)  shows  that  energy  (or  mass) 
addition  lags  pressure  whereas  with  ^  negative,  energy  addition 
leads  pressure.*  (22)  shows  a  decrease  in  the  mean  pressure  from 
the  steady-state  value  when  £  is  positive  and  an  increase  in  the 
mean  pressure  when  'fc'  is  negative.  This  means  that  when  energy 
addition  leads  (lags)  pressure,  the  mean  pressure  is  higher  (lower) 
than  the  steady-state  value.  This  is  in  qualitative  agreement  with 
the  results  of  the  second  model  of  Reference  6  where  a  lead  or  lag 
was  introduced  by  means  of  the  Crocco  time-lag  postulate  but  no 
shock  waves  were  present.  Hcwever,  in  that  case,  the  absence  of 
the  shock  caused  the  "blowup"  of  the  solution  as  phase  went  to  zero. 
Note  that  the  deviation  of  the  mean  pressure  from  the  steady-state 
value  is  of  the  order  of  where  ^  is  a  perturbation  quantity 

of  the  order  of  the  amplitude  and  ^>  =  0  [i/  ) .  Since  is  a  very 

small  number,  this  deviation  is  much  smaller  than  the  amplitude.  In 
Reference  6,  it  was  found  to  be  a  second  order  quantity. 


*A  small,  negative  I'  or  a  small  time-lead  has  the  same  physical 
effect  as  a  time-lag  with  a  value  equal  to  the  period  of  oscillation 
minus  a  small  number. 
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It  is  seen  from  (22)  that  if  r/X  were  negative,  ex¬ 
pansion  shocks  would  occur.  From  its  definition  r  has  the  same 
sign  as  [  A- 1  ]  ,  so  if  X  >  0  and  A  <  1  (or  CJ  1)  ,  no 

physically  reasonable  solution  is  found.  This  can  be  shown  to  be 
true  even  if  X  were  negative,  by  means  of  a  topological  investi¬ 
gation  of  the  solutions  to  (20)  (See  Reference  6).  This  means  that 
in  the  "linearly"  stable  region  CJ  <  1  (and  ^  small)  there  can 
be  no  oscillation  which  precludes  "triggering"  action  as  the  phase 
between  energy  addition  and  pressure  asymptotically  goes  to  zero. 
This  result  is  in  agreement  with  one  of  the  conclusions  of  Reference 
6  where,  based  on  the  trend  of  the  calculations,  "triggering" 
seemed  impossible  in  the  limit  as  phase  went  to  zero. 

The  period  of  the  oscillation  may  be  determined  from  (16) 
and  (21)  to  be  the  following 
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The  coefficient  of  Z'  in  (23)  can  be  shown  to  be  positive  and  quite 
small  in  the  range  of  interest.  If  Z'  were  zero,  (23)  would  give 
the  nonlinear  resonant  frequency  of  the  chamber  with  a  mean  flow 
present.  The  time-lag  causes  a  modification  of  the  frequency  as 
shown  by  (23).  In  particular,  if  >  0 ,  the  frequency  is  below- 
resonant;  and,  if  0,  above- resonant .  This  means  that,  if  energy 

addition  leads  (lags)  pressure,  the  oscillation  has  an  above  (below) 

resonant  frequency,  which  agrees  with  both  the  linear  results  of 

4  6 

Crocco  and  the  nonlinear  results  of  Sirignano  .  These  frequency 

modifications  are  extremely  small,  however. 

Some  mention  should  be  made  of  the  mean  mass-burning  rate 

under  oscillation.  Since  there  is  a  deviation  in  the  mean  pressure 

and  the  mean  velocity  from  their  steady-state  values,  it  indicates 

that  the  mean  burning  rate  also  deviates  from  its  steady-state  value 

This  is  not  possible  for  fixed-injection  systems  since  the  amount 
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burned  must  equal  the  amount  injected.  It  can  be  demonstrated, 
however,  that  with  the  pressure  (or  density)  deviation  of  the  order 
of  rp  ,  the  burning  rate  deviation  is  of  the  order  of  u  , 

So,  if  uq  is  small,  the  burning  rate  deviation  from  the  injection 
rate  is  negligible  compared  to  the  pressure  deviation.  Furthermore, 
uq  fa  is  always  a  very  small  number.  Therefore,  while  this 
deviation  produces  a  mathematical  inconsistency,  the  numerical 
errors  are  quite  negligible. 

In  summary,  it  can  be  said  that  the  results  of  this 
analysis  fit  quite  well  into  the  "structure  of  understanding'  to¬ 
gether  with  the  results  of  the  other  two  models  discussed  in 
Reference  6.  Namely,  the  impossibility  of  "triggering"  action  as 
phase  goes  to  zero  is  shown,  the  mean-pressure  is  shown  to  be  higher 
(lower)  than  the  steady-state  value  and  the  frequency  higher  (lower) 
than  the  resonant  frequency  if  the  energy  addition  leads  (lags)  the 
pressure. 


k 


k  k  k  k 
*1  *  *2  y  3  *  4 


p 

r 


t 

T 


u 

V 


x 

Y 


NOMENCLATURE 


speed  of  sound 

parameter  defined  after  Equation 

parameter  defined  after  Equation 

homogeneous  solutions  to  partial 
Equations  (2) 

parameter  defined  after  Equation 

combustion  parameters 

pressure 

parameter  defined  after  Equation 
t  ime 

period  of  oscillation 
gas  velocity 
shock  velocity 
longitudinal  dimension 
ratio  of  specific  heats 


(13) 

(20) 

differential 

(13) 


(20) 
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£ 

X 

tX 

E 

t 

* 

P 

o> 


combustion  parameter 

parameter  defined  after  Equation  (20) 

parameter  defined  after  Equation  (9) 

dummy  variable  of  integration 

argument  of  function  G 

time-lag 

parameter  defined  after  Equation  (20) 
combustion  parameter 


Subscripts: 

o  zero  order  quantity 

1,2  order  of  correction  to  solution 


Superscript : 
Primes: 


perturbation  quantity 
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D.  TRANSVERSE  NOZZLE  ADMITTANCE  RELATIONS  FOR  FINITE- AMPLITUDE  WAVES 
Introduction 

There  are  several  fundamental  differences  between  acoustical 
oscillations  In  a  closed  cylindrical  chamber  and  the  oscillations  that 
may  be  present  in  the  combustion  chamber  of  a  rocket  engine.  In  the 
latter  case,  the  oscillations  are  superimposed  on  a  mean  flow  created 
by  the  combustion  process.  While  in  acoustics  the  amplitudes  of  the 
oscillations  are  small,  this  may  not  be  the  case  with  the  oscillations 
that  may  occur  in  a  combustion  chamber.  In  addition,  the  classical 
closed-end  boundary  condition  used  in  the  acoustic  problem  is  now 
replaced  by  a  boundary  condition  imposed  by  the  presence  of  a  con¬ 
verging-diverging  nozzle.  When  the  nozzle  is  in  supercritical  operation, 
the  flow  down-stream  of  the  nozzle's  throat  is  supersonic  and  no  finite 
disturbance  can  travel  in  the  upstream  direction.  This  condition  can 
be  expressed  by  requiring  that  the  solutions  for  the  flow  field  in  the 
nozzle  be  regular  at  its  throat.*  The  regularity  condition  can  be  re¬ 
placed  by  a  complex  relation  between  the  pressure,  entropy  and  velocity 
perturbations.  This  relation  is  expected  to  hold  in  any  location 
along  the  converging  portion  of  the  nozzle.  This  relation  for  three- 
dimensional  perturbations  has  been  termed  the  transverse  admittance  re¬ 
lation.  Evaluated  at  the  nozzle  entrance,  it  forms  the  proper  boundary 
condition  for  the  oscillatory  flov:  in  the  combustion  chamber. 

The  problem  of  supercritical  flow  with  oscillations  in  a 

13 

converging-diverging  nozzle  was  first  treated  by  Tsien  who  considered 

the  case  in  which  the  oscillation  of  the  incoming  flow  is  isothermal. 

The  solution  was  found  for  both  very  low  and  very  high  frequencies. 
Crocco^  extended  this  study  to  include  the  non  isothermal  case  and 
covered  the  entire  frequency  range.  Both  of  these  investigations  were 
limited  to  one- dimensional ,  or  longitudinal,  oscillations.  In  a  later 

*  No  continuous  disturbance  can  travel  upstream  from  the  nozzle  throat 
unless  the  amplitude  of  the  disturbance  at  the  throat  is  infinite. 

Shock  waves  are  not  considered  in  the  present  analysis. 
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work  Crocco  extends  his  one- dimensional  analysis  to  the  case  of 
transverse  oscillations  where  he  considers  the  behavior  of  three-dimen¬ 
sional  perturbations  on  the  one- dimensional  flow  in  a  deLaval  nozzle. 

It  is  the  purpose  here  to  extend  the  above  work  which  is  limited 
to  small  perturbations,  to  the  case  where  the  perturbations  about  the 
one- dimensional  mean  flow  have  amplitudes  of  finite  size.  The  trans¬ 
verse  admittance  relation  resulting  from  the  presence  of  three-dimen¬ 
sional  perturbations  about  the  one-dimensional  steady  mean  flow  will  be 
derived.  This  expression  will  provide  an  appropriate  boundary  condition 
for  the  case  when  the  flow  oscillations  in  the  combustion  chamber  have 
amplitudes  of  finite  size.  The  latter  can  result  from  either  the  growth 
of  the  originally  "small"  disturbances  or  the  "introduction"  of  dis¬ 
turbances  of  finite  amplitudes.  The  admittance  relations  corresponding 
to  isentropic  or  irrotational  oscillations  can  be  obtained  as  special 
cases  of  the  general  theory. 

The  work  done  in  this  paper  is  concerned  with  the  flow  in  an 
axi- symmetric  nozzle.  It  could,  however,  be  extended  to  the  special 
case  of  a  two-dimensional  nozzle.  The  latter  case  will  not  be  treated 
here.  The  admittance  relation  which  will  be  derived  could  be  used  as 
a  boundary  condition  in  the  determination  flow  behavior  in  any  pro¬ 
pulsive  device  followed  by  a  converging-diverging  nozzle  which  is  oper¬ 
ating  in  the  supercritical  range. 

Flow  Equations 

In  this  derivation,  the  flow  in  the  converging  part  of  an 
axi- symmetric  nozzle  is  assumed  to  be  adiabatic  inviscid  with  no  body 
forces  and  no  chemical  reaction.  The  fluid  is  assumed  to  be  a  perfect 
gas  with  constant  specific  heats.  The  conservation  equations  in  di¬ 
mensionless  form  are  then  as  follows: 

Conservation  of  mass 

+  V- (ff)=  o  (2.1) 

Conservation  of  momentum 


I 


*  T  V(|>(V,/)7  =  - 


VP 


(2.2) 
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Since  the  fluid  is  assumed  to  be  Inviscid  and  non*  heat  con¬ 
ducting,  the  energy  equation  in  its  simplest  form  expresses  the  con¬ 
stancy  of  entropy  of  a  fluid  particle  after  it  enters  the  nozzle. 


JT+;  -  vs 


(2.3) 


where 


s* 

CP  * 


constant 


(2. A) 


gives  an  expression  for  the  nondimens ional  entropy.  The  equation  of 
state  for  a  perfect  gas  was  used  in  the  derivation  of  the  last  equation. 

In  the  equations  above,  the  reference  quantities  were  chosen 
as  the  stagnation  quantitlves  of  the  unperturbed  gas  entering  the  noz¬ 
zle.  The  unperturbed  flow  in  the  nozzle  is  assumed  to  be  isoenerget ic , 
isentropic  and  irrotat ional .  Consequently,  its  stagnation  properties 
remain  constant  throughout  the  unperturbed  flow  and  therefore  are  suit¬ 
able  reference  quantities. 

Hence  we  define: 


f* 


-.5  ?* 


t* 


(2.5) 


where  *  denotes  a  dimensional  quantity,  o  stagnation  values  and 
superposed  bar  unperturbed  (steady)  values.  L*  is  a  suitably  defined 
characteristic  length  and  is  used  in  the  non- d imens iona 1 izat ion  of 
quantities  having  the  dimension  of  length. 

Linearization  of  the  Equations 

In  analyzing  the  flow  behavior  of  a  supercritical  nozzle  in 
which  the  amplitudes  of  the  unsteady  perturbations  about  the  steady- 
state  mean  flow  are  finite  in  size,  it  is  convenient  to  assume  that  the 
dependent  variables  appearing  in  the  problem  can  be  written  in  the  fol¬ 
lowing  form: 
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Where  €  is  some  amplitude  parameter  and  the  variables  written 
above  are  respectively  the  non-dimensional  velocity,  pressure,  density 
and  entropy. 

The  expressions  appearing  in  (3.1)  are  substituted  into  equa¬ 
tions  (2.1)  through  (2.4)  and  the  equations  are  separated  according  to 
different  powers  of  €  .  The  zero  order  equations  describing  the  mean 
flow  conditions  are: 

V * 

+  (V  X  j  )  1  p  =  -  yp  (3.2) 

%  •  75*0 

S  '*  Y  Jh.  P-  jt*.  J>  +  constant 

These  equations  can  be  replaced  by  simpler  ones,  obtained  from 

(3.2)  when  the  flow  is  irrotat tonal : 


(3.3) 


Using  the  definition  of  the  non  dimensional  sonic  velocity 


C  = 


together  with  (3.4)  we  get  from  (3.3): 


2 


C 


(3.4) 


(3.3) 
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which  specifies  the  mean  density  variation. 

When  the  nozzle  is  axi- synsnetrlc  the  first  of  equations  (3.3) 
can  be  used  to  define  a  stream  function: 


7l  “  *  V  V  0.6, 

where  V  is  the  nondimens  Iona 1  distance  from  the  axis  of  symmetry  and 
the  unit  vector  in  the  tangential  direction. 

Undei  the  present  assumption  of  irrotationality  of  the  unper¬ 
turbed  flow,  a  potential  function  can  also  be  defined  as 


(3.7) 


The  stream  and  potential  functions  introduced  are  nondimen- 

sional . 


In  the  following  the  unperturbed  flow  will  be  assumed  to  be 
in  the  meridional  plane. 

The  equations  of  order  ©  and  their  solutions  are  identical 
to  those  obtained  by  Crocco  on  his  derivation  of  the  transverse  admit¬ 
tance  relation  for  the  case  of  small  amplitude  perturbations.  For  the 
sake  of  completeness,  these  equations  are  repeated  here: 


ff’+  *lf-  + 
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=  0 


v) 


+  jr  (3. 


(3.8a) 
9a) 


(3.10a) 

(3.11a) 


2 

The  equations  of  order  0  are  similar  to  those  found  by  Crocco  ex¬ 
cept  for  their  inhomogeneous  parts.  They  are  written  as  follows: 


— 
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(3.8b) 


r*) 
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a  TP 


(3.10b) 


(3.11b) 


A  comparison  between  the  systems  of  equations  corresponding 

to  0(C)  and  0(£*)  shows  that  their  homogeneous  parts  are  identical. 

The  solution  of  the  second- order  system  of  equations,  along 
lines  paralleling  the  solution  of  the  first-order  system  of  equations 
will  be  found.  Before  this  is  done  an  appropriate  set  of  independent 
coordinates  will  be  chosen. 

Choice  of  the  Independent  Variables 

Abandoning  the  vectorial  representation,  it  is  useful  to  choose 
the  Independent  variables  in  a  way  appropriate  to  the  introduction  of 
the  boundary  conditions  at  the  nozzle  walls.  In  the  ad-synmetric  case,  following 
Crocco,  it  is  convenient  to  let  the  steady-state  potential  function  <P 
replace  the  axial  variable,  and  the  steady- state  stream  function  replace 
the  radial  variable.  Indicating  by  ^s  and  $ n  elementary  (non- 
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dimensional)  lengths  In  the  direction  of  the  unperturbed  stream  lines 
and  of  their  normals  on  the  meridional  plane,  see  Figure  III-10,  equations 
(3.6)  and  (3.7)  can  then  be  written  as  follows: 


- 

l  "  7^ 


rS  i-  J? 


(4.1) 


The  third  independent  variable,  0,  Indicates  the  variation 
in  the  tangential  direction. 

Separation  of  the  Variables 

In  order  to  be  able  to  solve  equations  (3.8b)  through  (3.11b) 
the  following  assumptions  are  made: 

(i)  The  unperturbed  flow  is  one- dimensional  and  the  cor¬ 
responding  variables  describing  it  depend  only  on 

(ii)  Assumption  (1)  also  Implies  that  the  angle  of  obliquity 
of  the  stream  lines  with  respect  to  the  axis  of  symmetry  is  sufficiently 
small  so  that  its  cosine  is  practically  1  and  the  element  of  normal 

cTn  along  the  surface  (>  -  constant  can  be  identified  with  dr. 

Hence  equation  (4.1)  can  be  integrated,  providing 


(5.1) 


(iii)  The  solution  will  be  restricted  to  the  case  where  the 
axial  component  of  vorticity  is  zero.  This  assumption  is  necessary  in 
order  to  achieve  separation  of  variables  in  the  continuity  equation. 
This  assumption  is  equivalent  to  saying  that  and  %:>  <*>• 

defined  in  Equation  (5.8),  are  equal  to  one  another.  This  assumption 
will  be  further  elaborated  on  other  parts  of  this  report. 

Rewriting  equations  (3.8b)  through  (3.11b)  in  a  ((f)  >  ©  ) 

coordinate  system  and  separating  the  momentum  equation  (3.9b)  into  its 
three  components  results  in  the  following  system  of  equations: 


A 
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Note  that  the  dependent  variables  appear  in  equations  (5.2) 
through  (5.7)  in  the  following  particular  grouping: 


Consequently  it  will  be  more  convenient  to  solve  for  these  quantities 
rather  than  the  explicit  unknowns. 

The  inhomogeneous  parts  of  equations  (5.2)  through  (5.7) 
are  composed  of  solutions  of  the  first-order  equations  and  are  thus  as¬ 
sumed  to  be  known.  It  is  also  important  to  note  that  the  homogeneous 
parts  of  these  equations  are  identical  to  the  system  of  equations  of 
order  €  . 

Equations  (5.2)  through  (5.7)  have  two  kinds  of  solutions 
corresponding  respectively  to  tangentially  spinning  and  standing  waves. 
Since  the  second-order  equations  are  linear  and  the  standing  or  spinning 
type  solutions  can  be  expressed  in  terms  of  combinations  of  exponentials, 
it  will  be  easier  to  solve  the  second-order  equations  once  by  assuming 
that  the  dependence  of  the  solutions  on  time  and  0  has  the  following 
general  form:  21  ^Z^*****"^  ihe  final  result  is  obtained  by  assigning 

m  and  n  their  correct  values.  These  will  be  determined  from  the  form 
of  the  expressions  appearing  in  the  inhomogeneous  parts  of  the  second- 
order  equations. 

Since  the  homogeneous  system  of  the  second-order  equations 
is  separable  but  the  inhomogeneous  system  is  not,  it  becomes  necessary 
to  assume  that  the  unknowing  appearing  in  the  second-order  eqjations 
can  be  expanded  in  terms  of  the  eigenfunctions  which  satisfy  the  homo- 
geneous  system  (i.e.,  e  g  /  Iff  )  is  such  an  eigenfunction) 

or  equations,  the  boundary  conditions  at  and  the  con¬ 

dition  of  periodicity  of  the  solution  in  the  0  direction.  The  secon- 
order  solutions  are  thus  assumed  to  have  the  following  form:^^ 


(1)  The  various  functions  appearing  in  (5.8)  may  be  complex  functions 
of  their  arguments. 
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v/here  fa)  represents  the  nondimensional  angular  frequency,  which  relates 
to  the  dimensional  fiequency  as  follows: 


U>  = 


(5.9) 


and  m  and  n  take  on  a  particular  set  of  values. 

It  is  also  assumed  that  the  known  expressions  appearing  in 
the  inhomogeneous  parts  of  Equations  (5.2),  (5.3),  (5.6)  and  (5.7)  could 
be  expanded  in  an  eigenfunction  expansion  of  the  following  general 
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form: 


F^G,  (Mr).  F4  Lj>)  Z  A^'  CsHI 


(A) 


i»JQ 


where  e"''rw'  and  T {\ JT  )  are  eigenfunctions  of  the 
homogeneous  system  of  equations.  In  order  to  illustrate  the  manner  in 
which  the  above  expansion  is  to  be  used,  consider  for  example  the  ex¬ 
pansion  of  the  first  term  in  the  inhomogeneous  part  of  equation  (5.2) 
for  the  case  of  a  spinning  wave: 


Vur 


*0  K$  f  '  °  **  ^  /^7 


=  c«.zc„tfv6)jT/(s^  Op  ) 


Tn\J)  «  _  .  rr?0 


where 


R<1^and 


are  assumed  to  be  known  from  the  solution  of 


the  first  order  equation.  The  above  expansion  shows  that  m  and  n  will 
take  on  the  following  values: 


m  *  n  *  0; 


m  *  n  *  2 ; 


m  *  n  =  -  2 


The  constants  A_  _  and  A  .  can  be  obtained  by  expanding 

^  >0  ^  >  Q 


the  function  J,  (3^  W  )  in  two  Dini-series  one  in  terms  of  JQ  (S#  )^kL  ) , 
and,  the  other,  in  terms  of  J9  i  (S  .  V  Uf  ). 

9  rvi  t 

To  achieve  separation  of  the  equations  corresponding  to  0(C  ) 
it  is  also  necessary  to  assume  that  the  terms  appearing  in  the  inhomo¬ 
geneous  part  of  equation  (5.4)  (^  -  momentum)  could  be  expanded  as  fol¬ 
lows  : 


geic"ut”*lajc  ? 
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It  is  also  necessary  to  assume  that  the  terms  appearing  In 
the  inhomogeneous  part  of  equation  (5.5)  (#  -  momentum)  could  be  ex¬ 
panded  as  follows: 


f,c0  2^*0 1 


i  L>nut* 


r 


(C) 


Expansions  (B)  and  (C)  could  be  obtained  by  integrating 
F2(4>)cyt,  V,  <9  )  and  F3(£)G3(t,f,0  )  with  respect  to  ^  and  $ 


respectively  and  then  expanding  the  integrated  functions  in  a  Dini- 
series  of  the  same  form  as  expansion  (A)  on  page  47.  Differentiation 
of  the  resulting  series  with  respect  to  ^  and  Q  results  in  ex¬ 
pansions  of  the  same  form  as  (B)  and  (C).  It  is,  of  course,  necessary 
to  assume  that  the  integrated  functions  and  their  expansions  are  of 
such  nature  that  term-by-term  differentiation  is  allowable.  The  va¬ 
lidity  of  the  above  statements  is  also  important  in  proving  that 


The  latter  is  a  necessary  condition  for 


the  separation  of  second- order  equations. 

It  is  convenient  at  this  point  to  introduce  the  function 


<*> 


defined  as  follows 

tt‘v>  =  'Qfl 


N|> 


(5.10) 


Substitution  of  Equations  (5.8)  and  (5.10)  into  Equations 
(5.2)  through  (5.6),  expanding  the  terms  appearing  in  the  inhomogeneous 
parts  of  these  equations  in  the  manner  described  in  the  last  paragraph 
and  use  of  assumption  iii  (i.e.,l£^  )  result  in  the  fol¬ 

lowing  set  of  equation*  corresponding  to  particular  values  of  the  sub- 
,  /  >v 

scripts  (  'h  )  • 

d  .  _  d2 


/For  brevity,  the  following  notations  will  be  used:  1  = 


dr 


*In  the  derivation  of  the  following  equations  the  dimensions  were  con¬ 
veniently  scaled  so  that  V  ..  =  1. 

wa  i  1 


* 
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MlU/  9*$^)'- s\P)V*'i0+ 9*  =  £l%>  «t>)  (5,11) 

r  *  o*,)  2  <* V  M/)  /  t»ip r 


xi  or 


C/^r  Qt>  *  (5-i2) 


(5.13) 


K^+ 1'  >^»  ♦  Ct&  - 


-<*>  /. 


r-rv  ■%<* 

Q*-*' %»<»•*' 


(5.14) 


(5.13) 


(’■.</)  '  ~  •  -*»ij) 

The  functions  ^Ufj^  ^4|/^ 

are  defined  in  Appendix  A. 

It  will  be  interesting  at  this  point  to  consider  the 


(5.16) 


expression  for  the  second-order  vorticity  and  examine  more  carefully 
the  statements  made  in  assumption  iii  VXf  ■  ?*f  t  Va^c  t 
where  "* 


‘Vt... 


t 


u> 


_5£ 

»» 

J 

Ur 

»<p 

)Y 

a>) 

.  ,fk) 

J4_ 

t 

r?l 

r£, 


r*/ 
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-jr 


**  a 


* »  g  <*'  *"""*’?  (.<■*  -  <,>)  T« 

♦n*-5 ?  ( «*  -  &)  * 

(5.17) 

From  Equation  (5.17)  it  can  be  seen  that  the  requirement 

.  (\) 

introduces  the  same  restriction  as  in  the  first-order  equa¬ 


tions;  i.e.,  the  second-order  vorticity  corapcnent  along  the  streamlines 
must  be  identically  zero.  Further  implications  of  above  restriction 
will  be  discussed  later  on. 


Reduction  of  the  System 

For  convenience,  the  superscript  (2)  and  subscripts  (»*7,  £  ) 
used  in  previous  sections  will  be  left  out  of  the  following  derivation. 
Substract  equation  (5.14)  from  (5.13)  to  get: 


g  *  (w-v)'  4  *** us  (w-v)  -  JO  -  C  =  O 


(6.1) 
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i 


(Because  3D  -  C  .  See  definition  of  these  quantities  in  appendix  A 
and  discussion  on  bottom  of  page  4g  . 

The  general  solution  of  (6.1)  is 


U/-V  =  c  f 

O  0 


<>n) 


(6.2) 


where 


On) 


i«t> 


)  J  ->niur  r 

(^)  =  e 


(6.3) 


The  constant  C  must  be  taken  equal  to  zero  in  order  that 

o 

V  *  W,  which  is  a  necessary  condition  for  the  separations  of  the  continuity 
equat ion . 

From  (5.15)  and  (5.16)  we  get: 


5- 

0 


-  b 


—  I 


Pop)-  C  R«p) 

where  cr  »  S(<f>  mo) 


f  1  &  (<p) 


(6.4) 


Integrate  (5.12)  to  get: 


9  1  $  ♦  P  •  i  f  ^  S6fJ.iL  *  f 

1  J*  i  i 


=T 


r 

Jf  • 


^6 


=.  <rf  +  K 


(6.5) 
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where 


The  pressure  is  eliminated  between  (6.5)  and  (5.13)  to  get: 

(<£-^=3^  -t  K  -  C  =  T  f(  ♦  L 

where 


L  -  K  -  C 


The  solution  of  the  above  equation  is: 


c 


[  rck 


(6.7) 


where 


F,  =  j  ‘1nd  71-f. 


(6.8) 


From  (6. A) 


ptw -  E?„t  C‘c^c‘  Jf*  ffV  *w-*‘*> 

Using  the  last  expression  for  P(^)  in  equation  (6.5)  yields  the  fol¬ 
lowing  equation: 

Jl  C‘  R«  V-f,  -  C  C*  ^  *  6  *  K 


-  tft  -  c  z-flr")  +  5  =  N 


M 
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where : 


and 


-  -  _  t>nj  rf> 

X  =  K  -  c  *  f  jA'f  € 

^  *  J  rv>  r 

O  0  r 


N-tr  (*-£*£">  1 


Combining  (5.11)  and  (6.7)  gives  the  following  relation: 

+  .11  ^  tc, V"!)  *1 s,^ 

-M-  -<5 

where  Q  -  ‘  )  ^  SA,  +  E‘  -M  ^  Si,; 

Eliminating  R  between  the  above  differential  equations  results  in 
the  following  second-order  differential  equation  for  $  : 

+  *iiw.('^c,)(‘^L)iZl-  j£is*  H-i't'btc*') 

J  0  “  **/  /  (ft 


(6.9) 


Defining 


r  =  — -s — 

and  u^ing  the  definitions  of  Q,  N  and  M  and  the  relations 
Jf  ("° 

—  t  ®Ia_  -t-  >***CJ'  -K 


(6.10) 


r  w* 


(6.11) 


and 


wr'^  » jW) 

[*s- +tr-f«)] 


~  >n»  w 


(6.12) 
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the  inhomogeneous  part  of  equation  (6.9)  becomes: 


-  z‘f~  T[r$-  ^  ^'^-n 

Now  equation  (6.9)  can  be  written  in  the  following  form 

C*f  (£)«  -  <r  F‘*0*  F‘* 

(6.13) 

where  <**(•$)  represents  the  left-hand  side  of  equation  (6.9)  and  the 
functions  F^  ^  J  *  1,  2,  3,  are  expressed  in  terms  of  known  functions 


F'"-  V 

r-  -  f  .  sju  jf 


1  ri 


(6.14) 


(E  H.  sitt 


*-(•<•  A  (-ft 


Admittance  Condition  at  an  Entrance  of  an  Axisymmetr ic  Nozzle 

The  general  solution  of  equation  (6.9)  is  a  combination  of 
the  general  solution  of  the  homogeneous  equation  and  the  particular 
solution  of  (6.13)  If  these  solutions  are  known  the  general  solution 
can  be  written  as  follows: 

§=  C.§°>  +  <r$ $>  +  c,  $1,  *  C,  (7.1) 

where  ^  and  are  two  independent  solutions  of  the  homogeneous 

part  of  equation  (6.9): 

^(§h)-o 

C„  and  arc  arbitrary  constants.  Now  note  that  above  equation  has 
the  following  singular  points: 


f.o.  f.e.cthrMt 

For  supercritical  nozzle  with  finite  nozzle  entrance  only  the  singularity 
at  the  throat  (where  (fa  *  O  )  is  of  interest. 


4 — )*  and 


•  co 


I 
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^ssuming  that  all  the  singularities  of  the  solution,  $  ,  ap¬ 
pear  in  then  the  condition  requiring  the  regularity  of  the  solution 

at  the  nozzle's  throat  can  be  expressed  by  requiring  =  0.  Consequently 
the  proper  solution  of  (6.9)  is: 

4*1  T  <»» 

*  C.  9, 

(7.2) 

From  (7.2)  together  (5.10),  (6.7),  (6.4),  (6.5),  and  assumption 
iii,  the  following  relations  are  obtained: 

-iff  «■*#* 

'  «  >.(<)  /  r 


p+  fl  U.  +  Mi  «*  £  K  *”C,  Miw  -  C£  wuur 


*(M) 


+  <r 


4-  O 


(7.3) 


If  these  are  considered  to  be  a  system  of  linear  equations 
for  the  determination  of  C^.O*  and  C 2  with  given  values  of  U, 

V  *  W  ,  P  ,  and  S  for  a  given  ,  we  see  that  only  three  of  these 

values  can  be  arbitrarily  prescribed.  In  other  words,  a  relation 
must  exist  between  any  such  set  of  four  values.  This  relation  is 
obtained  from  the  relation  of  compatibility  of  the  four  linear 
Equations  ;  i.e.,  from  the  vanishing  of  the  following  determinant; 


u-if. 


v-r’VM 


p*  f  'u 


X  »» 

$  -K 


-r'-  j* 

•  J 


(§,0-O  (*"H> 


fW 


■  M) 


Ji 

Jp 


—  >niur  $ 


—  O 


(7.4) 
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u  £*(C*  5“’-5)-^  ♦  e'tfVttt1  r'-j?) 


*»<•>  , 

5*  +  >»lur 


•4  >nit^ 


-  z  c* 


-  M  >nlur  c*  s*4  +  K  (£*  6‘4-  5; 

r-J'l  V*  ♦  ~«r  tl*  -  {  5] 


(7.6) 


Above  relation  holds  for  any  value  oi  .In  particular, 

(7.6)  it  holds  at  the  nozzle's  entrance  where  it  provides  the  admittance 
relation  for  the  combustion  chamber  flow.  For  the  second-order  equations, 
there  are  infinitely  many  admittance  relations  corresponding  to  differ¬ 
ent  values  of  the  subscripts  (nV,«)  appearing  in  the  eigenfunction  ex* 
pensions  used  for  their  solution.  The  function  $  which  appears 

in  the  Inhomogeneous  part  of  the  admittance  relation  is  the  particular 

(3) 

solution  of  equation  (6.13)  which  corresponds  to  F  (defined  in  equa¬ 
tion  (6.14))  and  it  is  assumed  to  be  known. 

The  transverse- admittance  relation  for  the  first-order  equa¬ 
tion  can  be  obtained  from  the  second- order  admittance  relation  by  putting 
its  Inhomogeneous  part  equal  to  rero  and  letting  m  ■  1. 


The  admittance  relation  simplifies  in  the  isentropic  case  by 
taking  S  ■  0  identically.  If  the  perturbed  flow  is  irrotational ,  it 
must  also  be  isentropic.  This  Implies  that  ■  ■  0  and  con¬ 

sequently  Equations  (7.3)  become: 

V'  ^‘*>m  Cj 

Pt  !/♦  $4  -  K*  "  Ct  $4  (7.7) 

The  above  relation  will  be  simplified  further  by  observing  that  now  the 

expressions  for  J  and  will  be  considerably  simpler. 
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Eliminating  C,  frcm  (7.7),  one  obtalnes: 


V-  5‘" 


<D 


i. 


flr) 

P+f'U+* ii<r 


->ncu)' 


(7.8«) 

(7.8b) 


which  represents  two  admittance  conditions  when  applied  to  the  entrance 
of  the  nozzle. 

For  purely- axial  oscillations,  Equations  (5.13)  and  (5.14) 
become  meaningless  and  one  has  to  disregard  Equation  (6.7)  which  Is  re¬ 
lated  to  Equation  (5.13).  The  corresponding  admittance  relation  can  be 
obtained  from  (7.3)  by  disregarding  the  second  equation  and  all  terms 
in  .  Using  this  procedure,  one  finds  the  following  admittance  re¬ 
lation: 


U  £  ♦  Xsar)  +  P  5  *■  5  C  * 


_  Mt  up- 


k  S  ♦  z  e* 


(7.9) 


For  one- dimensional  isentroplc  oscillations  above  admittance  relation 
is  replaced  by: 


(7.10) 


Note  that  this  equation  has  the  same  form  as  one  of  these  ap¬ 
pearing  in  (7.8b),  although  the  similarity  is  only  formal  since  the 


values  of  $  (or  iu  )  depend  on  entirely  different  equations. 
Note  further  that  these  nonlinear  admittance  relations  may  be  used 
only  for  an  oscillation  without  shock  waves.  It  is  very  likely 
that  shocks  would  be  present  in  a  purely  axial  mode. 


NOMENCLATURE 


r  - 

p  - 

s 

l  * 

u  - 

V  - 
w 

R(i)- 

P<1>- 

U(l)- 

V<1>- 

w(i)- 

c  - 

r 

t 

Y  - 

4  - 

e  - 

w*  - 

r  - 

c 

p 

$9  - 

$» l  - 


Dcnsity 
Pressure 
Entropy 
Velocity 

Axial  component  of  velocity. 

Radial  component  of  velocity. 

Tangential  component  of  velocity. 

Function  indicating  axial  variation  of  ( 

Function  indicating  axial  variation  of  (C4p 

Function  indicating  axial  variation  o£(Uyf  0  ) 

Function  indicating  axial  variation  of  (V/'rfY^ 

Function  indicating  axial  variation  of  (  Cu^V  ) 
Speed  of  sound . 

Distance  from  axis  of  symmetry. 

Time 

Stream  function. 

Potential  function. 

Angle,  indicating  tangential  variation. 


Frequency 

Ratio  of  specific  heats. 

Specific  heat  at  constant  pressure. 


Infinitesimal  variation  in  axial  direction  (along 
Infinitesimal  variation  in  radial  direction  (along 


Eigenvalue,  solution  of  the  equation 


P 


C  -  Amplitude  parameter. 

-  Unit  vector. 

C  -  Entropy  perturbation  at  the  nozzle  throat. 
«C(  )-  Second  order  differential  operator 


constant) 

constant) 
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*(#  -  Defined  as  u#  = 

C1  ,  C  ,  -  Constants  appearing  in  the  solution  for  $ 


NOTE: 

The  functions  /A,  B.  C ,  f),  C  and  C  are  defined  in  Appendix  A. 
The  functions  f^® ,  K,  L,  M,  N,  Q,  X  and 

(j"l,2,3)  are  defined  in  sub  section  D-6. 

Z  ,  £  and  (J-1,2)  are  defined  in  sub  section  D-7. 

SUPERSCRIPTS: 

*  -  Dimensional  quantity 

o  -  Stagnation  quantity 
—  -  Steady  8tate  quantity 

(i)-  i  ■  1,  2  .  .  .  .  indicates  order  of  the  equation. 


SUBSCRIPTS: 


f»^) 


w 

h 


-  n  -  0,  +  1,  +  2,  .  .  .  g  ■  1,  2,  .  Indicates 

mode  under  consideration. 

-  wall 

-  Homogeneous 

-  Throat 

-  Indicates  direction  of  unit  vectors. 


I 
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IV.  NONLINEAR  COMBUSTION  INSTABILITY  MECHANISMS:  EXPERIMENTAL 

A.  TRANSVERSE  MODE  NONLINEAR  ROCKET  MOTOR  STUDIES 

For  several  years  parametric  studies  of  transverse  mode 
nonlinear  combustion  instability  have  been  carried  out  at  Princeton 
University  using  9- inch  diameter  rocket  hardware.  These  pulsed- limits 
studies  have  been  Instrumental  in  defining  problem  areas  for  more  basic 
study  as  well  as  reproducing  in  model-size  rocket  chambers  certain  as¬ 
pects  of  the  combustion  conditions  found  in  high- thrust  engines.  Trends 
in  the  variations  of  the  Instability  regimes  encountered  in  this  test¬ 
ing  point  to  possible  means  of  avoiding  certain  linear  and  nonlinear 
high  frequency  instability  problems  in  full-size  hardware  designs. 

A  number  of  important  parameters  were  outlined  in  last  year's 
report ^and  will  be  briefly  summarized  here  before  proceeding  to  the 
current  studies.  Using  the  ethyl  alcohol,  liquid  oxygen  propellant 

combination,  together  with  a  pulsing  device  based  on  a  grain  charge 

1-3 

and  burst  disk  arrangement,  the  incidence  of  nonlinear  combustion  in¬ 
stability  was  influenced  by: 

1. )  Orientations  of  the  spuds,  where  tangential  injector 

orifice  alignment  (radial  spray  fans)  showed  major 
differences.  Extremes  of  stability  were  shown  by 
6x2  (six  groups  consisting  of  two  spuds  similarly 
oriented),  the  most  stable  arrangement,  and  12  x  1 
and  1  x  12  groupings  which  proved  to  be  the  least 
stable . 

2. )  Rotation  of  the  spuds  to  the  radial  orifice  align¬ 

ment  was  shown  to  markedly  improve  the  stability 
when  compared  to  the  1  x  12  tangential  arrangement. 

3. )  Inter- fan  spacing  (distance  between  fuel  and  oxidizer 

spray  fans  in  a  single  spud)  changes  influenced 
stability  with  the  zero  spacing  arrangement  clearly 
superior  (i.e.,  fuel  fan  end- Impinging  on  the  oxidizer 
fan) . 

A.)  Chamber  pressure  variations  indicated  that  increased 
pressures  favored  higher  modes  and/or  linear  insta¬ 
bility. 


and 


5.)  Inter- Spud  spacing  changes  showed  that  larger  spacing 
arrangements  were  more  difficult  (or  impossible)  to 
trigger  into  nonlinear  combustion  instability. 
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These  studies  have  been  continued  for  both  the  ethyl  alcohol, 
oxygen  and  RP- 1 ,  oxygen  propellant  combinations.  Certain  additional 
data  on  these  past  five  parameters  have  been  obtained  and  will  be  cov¬ 
ered  first. 

Before  proceeding  to  the  details  of  the  test  dats  It  would 
be  well  to  review  the  mechanism  believed  responsible  for  many  of  the 
effects  observed..  Associated  with  transverse  pressure  oscillations 
are  velocities  and  hence  particle  displacements.  Stability  of  the 
combustion  process  Is  Involved  with  these  displacements  since  the  ef¬ 
fect  Is  selective  l.e.,  the  more  readily  vaporised  propellant  Is  dis¬ 
placed  with  reference  to  the  droplets  of  the  other  thus  varying  the 
local  mixture  ratio  and  hence  the  local  combustion  rates.  The  theo¬ 
retical  treatment  for  this  phenomenon  (first  observed  by  the  research 
group  at  Rocketdyne)  was  covered  by  Reardon  In  Reference  16  .  This 
topic  Is  also  covered  In  a  soon  to  be  published  paper 

The  situation  existing  In  actual  rocket  engines  from  the 
standpoint  of  determining  the  mass  and  mixture  ratio  distributions 
of  the  propellants  Involved  and  then  charting  the  course  of  the  dis¬ 
placed  vapor  becomes  quite  complex.  In  the  studies  at  Princeton 
University  the  propellants  are  normally  Injected  at  a  single  Injection 
radius.  Distribution  data  for  this  relatively  simple  arrangement 
will  be  discussed  in  detail  in  this  section.  The  displacement  paths 
and  how  they  vary  with  distance  from  the  center  of  the  Injector 
(originating  as  a  circle,  evolving  Into  an  eclipse  and  finally  Into  a 
circular  arc  at  the  chamber  wall)  also  complicates  the  situation. 

A  discussion  of  these  paths  and  the  Inter  relationships  with  elllplti- 
cal  spray  fans  of  various  orientations  Is  to  be  found  In  Appendix  A 
of  Reference  1. 

Thus  with  regard  to  the  various  effects  that  are  discussed 
next,  it  must  be  remembered  that  due  to  the  Interwoven  considerations 
as  mentioned  above  trends  In  stability  characteristics  ss  Influenced 
by  the  parameter  variations  will  be  Investigated.  The  true  magnitude 
of  any  of  these  effects  could  only  be  ascertained  In  the  final  hard- 
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ware  design  with  numbers  of  radial  injection  locations. 

That  the  stability  threshold  of  a  rocket  chamber  can  be  in¬ 
fluenced  by  increasing  or  decreasing  the  inter-spud  spacing  was  shown 
in  Reference  1.  Perhaps  the  best  illustration  of  these  spacing  in¬ 
fluences  was  found  in  Figure  18  of  that  Reference  and  is  repeated  here 
for  comparison  as  Figure  IV- 1.  Here  the  standard  instability  regime 
can  be  seen  as  determined  for  the  1  x  12  grouping  (1.84  inch  inter- 
spud  spacing)  as  compared  to  the  completely  stable  mixture  ratio- 
perturbation  amplitude  field  for  the  1  x  8  (2.76  inch  spacing)  and 
completely  unstable  field  (of  Increased  amplitude)  for  the  1  x  16 
(1.38  inch  spacing). 

To  evaluate  whether  the  chamber  Mach  number  was  influenclal 
in  establishing  the  degree  of  stability  in  the  inter-spud  spacing  in¬ 
vestigation  was  the  purpose  of  a  series  of  tests  in  progress  at  the 
time  of  the  last  report.  One  group  of  tests  used  spuds  with  orifices 
of  twice  the  standard  area  (.084M  diameter  vs.  .059"  diameter).  These 
spuds  when  placed  in  the  8  spud  configuration  providing  the  same  total 
propellant  flow  as  16  standard  spuds.  Thus,  the  Mach  number  was  main¬ 
tained  constant  while  spacing  and  orifice  sizes  were  the  variables. 
Testing  at  the  150  psi  level  with  these  double- size  orifice  area 
spuds  yielded  the  same  stable  results  indicated  in  Figure  IV- 1  for 
the  1x8  s:andard  spud. 

Another  approach  taken  to  provide  equivalent  chamber  Mach 
number  conditions  involved  the  use  of  a  12  spud  injector  plate  to 
which  4  more  centrally  located  spuds  were  added  (referred  to  as 
1  x  (12  +  4),  Figure  IV-2).  Such  testing  preserved  the  12  spud  spac¬ 
ing  (1.84  Inches)  while  producing  the  Mach  number  associated  with  the 
16  spud  arrangement.  The  data  obtained  in  these  tests  are  shown  in 
Figure  IV-2  and  IV-3.  In  the  0— +F  (see  nomenclature)  tests  (Figure 
IV-2)  the  Instability  regime  is  actually  slightly  smaller  than  that 
found  in  Figure  IV-1  for  1  x  12.  Checking  with  the  instability  limits 

m 

data  in  Reference  1  for  0~*F  (also  in  Figure  IV-2)  and  for  com¬ 
parisons  for  the  higher  pressure  results  (Figure  IV-3)  the  following 
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can  be  said:  At  150  psia  the  0— *F  data  also  indicates  a  somewhat 
more  limited  range  of  instability  as  compared  to  that  previously  re¬ 
corded  for  1  x  12  (Figure  19,  Reference  1),  certainly  not  even  ap¬ 
proaching  the  1  x  16  case.  At  300  and  600  psi  the  instability  regimes 
and  modes  closely  compare  to  the  previously  obtained  1  x  12  data 
(Figure  26,  Reference  l)with  the  characteristic  standing  modes  and 
higher  modes  entering  the  scene  as  the  pressure  was  raised. 

The  conclusions  based  on  both  the  large  injector  orifice 
tests  and  the  added  spud  tests,  in  which  the  chamber  Mach  number  for 
the  16  spud  arrangement  was  duplicated  in  each  case  are  as  follows: 

1.  Inter-spud  spacing  has  been  shown  to  be  an  important 
stability  parameter  for  the  150  psia  chamber  pressure 
regime. 

2.  Minor  differences  indicated  that  Mach  number  influences 
were  present  but  were  not  a  controlling  characteristic. 

Reardon  has  discussed  the  influence  of  thrust  per  element 

18  19 

(F/E)  in  recent  Aerojet  reports  ’  .  There  exists  a  close  con¬ 

nection  between  (F/E)  and  inter-spud  (or  inter-element)  spacing.  In 
a  fixed  thrust  rocket  design  both  parameters  must  change  simultaneous¬ 
ly,  whereas  in  the  Princeton  tests  just  described  the  parameters  were 
independently  varied  over  a  limited  range,  i.e.,  F/E  varied  from  80 
to  160  lbs/element  while  the  spacing  (considering  all  the  tests) 
varied  from  1.36  to  2.76  inches. 

Other  spud  arrangements  to  test  non-uniform  inter-spud 
spacings  were  also  investigated.  Both  16  spud  and  12  spud  injector 
face  plates  were  used  in  these  studies.  The  approach  was  that  of 
leaving  certain  selected  spud  locations  blank.  For  example,  with  the 
16  spud  face  plate,  groups  of  3  spuds  followed  by  1  blank  space  were 
used  and  assigned  the  nomenclature  4  (3  +  1  B).  The  purpose  of  these 
tests  was  to  determine  whether  non-uniformities  in  the  propellant  in¬ 
jection  (spaced  appropriately)  could  alter  the  instability  regimes. 

Evaluating  the  4  (3  +  1  B)  tests  for  the  four  pulse  direction- 
propellant  placement  arrangements  (see  nomenclature  this  section)  the 
following  points  are  observed  from  Figures  IV-4  and  -5.  0— *F  tests 

produced  the  broad  instability  regime  characteristic  of  the  1  x  16  tests 
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(Figure  IV- 1)  with  decreased ‘amplitudes.  0  ■  »F  pulsing  (also  Fig¬ 
ure  IU-4)  showed  instability  to  2.3  mixture  ratio  similar  to  1  x  16 
data  .  F  — ♦O  studies  (Figure  IV-5)  failed  to  produce  the  re¬ 

versed  spin  Instability  shown  in  both  1  x  16  and  1  x  12  results^. 

And  finally  F  — *0  teats  indicated  the  exlstance  of  instability  in 
separate  regions  at  high  and  low  mixture  ratio  where  1  x  12  tests  were 
entirely  stable  and  1  x  16  tests  were  unstable. 

In  another  series  of  tests  the  12  position  injector  face 
plate  was  used  with  2  spud  and  1  blank  groupings.  Testing  at  the  150 
psla  level  Indicated  stable  operation  under  0  *^F  pulsed  conditions. 
This  would  correspond  to  1  x  8  results.  Thus  when  the  spacing  alter¬ 
nates  between  each  pair  of  spuds  (i.e.  1.84,  3.68,  1.84,  3.68  etc.) 
the  blank  space  appears  to  be  effective. 

Based  on  the  stable  behavior  encountered  at  150  psia,  tests 
were  also  run  at  300  psla.  These  higher  chamber  pressure  results  fes 
shown  in  Figure  IV-6)  Indicate  the  presence  of  spinning  or  standing 
first  tangential  mode  oscillations  for  0-*F  and  0  ~*F  pulse  gun  and 
propellant  placement.  For  comparison,  the  300  psla  tests  with  the 
1x8  double-size  orifice  area  spuds  are  summarized  in  Figure  IV-7. 
These  results  also  indicate  the  same  stability  behavior  although  of 
somewhat  higher  amplitude.  Spot  checks  with  the  1x8  arrangement 
and  the  standard  orifice  size  reduced  the  amplitudes  close  to  those 
found  in  the  testing  of  the  4  (2  +  1  B)  injector. 

In  a  final  series  of  tests  checking  the  effect  of  non- 
uniform  propellant  injection  on  combustion  instability  regimes  the 
spacing  was  Increased  to  a  two- spud  gap  at  180°  intervals.  Using  the 
16  position  injector  face,  the  arrangement  consisted  of  4  spuds  fol¬ 
lowed  by  two  blanks  or  2  (4  +  2  B).  The  test  data  are  presented  in 
Figures  IV-8  and  -9.  Although  possessing  the  same  number  of  spuds  as 
the  last  tests  described,  the  stability  behavior  was  entirely  differ¬ 
ent  at  150  psia  chamber  pressure.  First,  tangential  standing  mode 
oscillations  were  present  for  all  four  pulsings  conditlons(0  — ►F, 

0— -►F,  F— #0  and  F— *0)  except  for  the  lower  mixture  ratios  opera- 
*The  pressure  antinodal  line  passed  through  the  blank  spaces. 
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Cion.  Even  a  number  of  Che  lower  mixCure  raclo  Ceaca  damped  very 

alowly  (one  Caking  600  mllliaec.*  Figure  IV-8).  Thua  1C  would  ap» 

★  2,20 

pear  chac  Che  wider  gap  acCa  much  like  a  baffle  In  alcerlng  Che 

chamber  envlronmenC  Co  one  favoring  Che  aCandlng  type  oacllladona. 
Furthermore,  It  again  haa  been  ahown  Chat  four  apuda  In  a  group 
(under  Che  general  rule  of  chree  or  more)  behave  much  aa  Che  full 
compllmenC  (In  chla  caae  a  16  apud  Injector). 

Ralalng  the  chamber  preaaure  to  300  pala,  linear  flrat 
tangential  standing  mode  oaclllaClona  replaced  the  nonlinear  type. 

AC  600  pal  linear  lnatabllicy  peraiated  with  flrat  tangential,  flrat 
radial  and  mixed  IT  1R  modea  preaent.  Aa  haa  been  prevloualy  dla- 
cuaaed  both  In  thla  section  and  in  Reference  1  ,  thla  transition  first 
Co  standing  mode  oscillations  and  then  to  higher  modes  appears  to  be 
characteristic  of  the  Injection  elements  used  (data  in  Figure  IV-9). 

Briefly,  the  following  conclusions  can  be  drawn  about  this 
variety  of  spacing  tests. 

1.  Similar  to  the  effects  observed  with  orientation 

If  three  or  more  spuds  are  grouped, followed  by  a  spac¬ 
ing  change,  only  minor  Improvements  in  stability  are 
derived.  On  the  other  hand,  when  alternating  spacing 
was  investigated  (l.e.,  two  spuds  followed  by  one 
blank)  stability  was  improved. 

2.  Use  of  wider  gaps  in  the  propellant  distribution  (l.e. 
two  spuds  positions  blank  at  180  intervals)  resulted 
in  effects  similar  to  those  observed  with  the  diametral 
baffle,  l.e.,  previous  spinning  first  tangential  mode 
was  replaced  with  standing  IT  mode  oscillations. 

3.  Higher  pressures  (300  and  600  psla)  resulted  in  the 
same  trends  as  the  previous  studies  reported  in  Ref¬ 
erence  1  .  The  linear  instability  tendencies  increased 
with  pressure  as  did  the  incidence  of  higher  modes, 

2nd  tangential  or  1st  radial  and  combined  modes. 

From  a  practical  design  standpoint  certain  difficulties  would 

be  present  in  using  voids  in  the  propellant  injection  density  to  iro- 

prove  stability  characteristics.  One  would  be  the  problem  to  locally 

cool  the  injector  face  if  no  propellant  was  being  Injected.  This  region 

may  receive  equivalent  heat  flux  or  even  Increased  heat  transfer  (if 
♦Appendix  B  provides  a  now  approach  to  jimple  baffle  design. 
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a  strong  recirculation  pattern  is  established  in  that  location) 
and  would  demand  careful  design  attention.  Such  considerations 
would  also  apply  to  any  form  of  increased  spacing  and  must  be  allowed 
for  in  the  Injector  design  (i.e.,  no  stagnation  regions,  sufficient 
cross  velocities  in  the  ring  groves,  proper  material  choice  to  reduce 
hot  spots,  etc.). 

RP-1  Testing 

In  previous  rocket  testing  at  Princeton  jet  fuels  and 

20b 

especially  a  purer  hydrocarbon,  iso-octane  ,  were  used  with  liquid 
oxygen.  The  tests  referred  to  were  concerned  with  longitudinal  mode 
combustion  instability  and  the  results  were  compared  to  those  obtained 
testing  with  ethyl  alcohol.  The  conclusions  based  on  linear  combustion 
instability  data  were  that  the  instability  regimes  were  similar  but 
slightly  smaller  for  LOX,  iso-octane.  This,  incidently,  could  be 
predicted  by  the  droplet  and  wake  burning  theory  covered  in  Section 
V.  Therefore,  similar  values  for  the  combustion  parameters  and 
n  were  found. 

In  considering  RP-1  for  nonlinear  combustion  instability 
testing,  other  factors  must  also  be  taken  into  account. 

In  determining  any  vapor  displacement  effects,  how  slowly 
the  droplets  of  one  propellant  vaporize  accomplishes  similar  ob¬ 
jectives  as  the  rapid  vaporization  of  the  displaced  propellant. 

From  this  standpoint,  RP-1  would  appear  to  be  somewhat  more  prone 
to  vapor  displacement  effects  than  ethyl  alcohol  due  to  its  slower 
vaporization  rate. 

In  the  first  RP-1  tests  to  be  described,  the  identical 
injector  elements  were  used  as  in  the  alcohol  tests,  i.e.,  like-on- 
like  doublets  with  .059  inch  orifices.  Since  the  nominal  mixture 
ratio  for  LOX,  RP-1  is  2.3  (maximum  c*  point, versus  1.4  for  LOX, 
alcohol).  The  LOX  velocity  is  increased* producing  finer  droplets 
near  the  injector  face  (where  D^q  *  13^1 A  versus  16(^  for  L0X- 
alcohol  assuming  Av-Vj,  Equation  11,  Ingebo  )  and  hence 
providing  a  means  for  more  rapid  vaporization.  The  fuel  velocity, 
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on  the  other  hand,  la  decreased  thereby  increasing  droplet  sizes 

(°30  RP-l"  verau*  D3o  Ale"  155/<,  Both  chan8**  may  be 

looked  upon  as  favoring  displacement  effects  provided  that  the  com- 

* 

bustion  distribution  has  not  been  significantly  altered  in  the 
combustion  chamber.  Tests  are  currently  in  progress  using  the 
appropriate  sizes  of  injector  orifices  necessary  to  restore 
equivalent  stream  velocities  at  the  2.3  nominal  design  mixture  ratio. 
These  tests  will  be  discussed  later. 

The  injector  chosen  first  was  the  1  x  12  tangentially 
oriented  fan  arrangement  -  the  "standard"  for  the  pulsed-limits 
testing  with  LOX,  alcohol.  Equivalence  ratio  is  used  to  cross 
correlate  the  RP-l  and  alcohol  data.  It  would  appear  from  these 
tests  (for  all  four  pulse  conditions)  that  the  instability  regimes 
were  noticeably  widened.  These  data  are  shown  in  Figures  IV-10 

and  IV-11.  Only  the  0  *F  tests  indicated  a  small  instability 

regime  expansion  (<C  .2  mixture  ratio  at  the  low  mixture  ratio  end). 
The  amplitudes  were  characteristically  of  the  order  of  260  psl  peak- 
to-peak  approaching  the  equivalence  ratio  and  amplitude  conditions 
experienced  with  LOX,  alcohol  and  the  closer-spaced  1  x  16  arrange¬ 
ment.^  Higher  pressure  testing  with  this  arrangement  produced 
linear  instability  of  both  spinning  and  standing  first  tangential 
inodes  similar  to  so  many  of  the  higher  pressure  tests. 

Use  of  the  6x2  orientation  approach,  which  had  proven 
to  be  quite  successful  in  limiting  the  instability  regimes  with 
alcohol  and  LOX,  was  checked  next  with  RP-l.  Spinning  first  tan¬ 
gential  mode  instability  was  found  with  certain  propellant  and 
pulse  gun  arrangements.  Conditions  of  complete  stability  occurred 
a  number  of  times  in  mixture  ratio  regions  formerly  unstable 
depending  on  the  exact  location  of  the  pulse.  The  particular 
data  referred  to  is  shown  in  Figures  IV-12  to  -14  for  a  number 

*If  the  larger  fuel  droplets  produced  should  result  in  the  rate 
of  combustion  decreasing  near  the  injector  face  then  the  additional 
time  for  mixing  and  evening  out  of  fuel  and  oxidizer  concentrations 
that  would  take  place  axially  could  result  in  a  supression  of  such 
displacement  effects. 
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of  pulse  gun  and  propellant  arrangements  at  150  and  300  psl. 

Compared  to  the  1  x  12  LOX,  RP-1  tests  these  data  represented  im¬ 
proved  stability  but  fell  considerable  short  of  the  near  stable 
operation  at  150  -  600  psia  chamber  pressure  that  had  been  indicated 
with  LOX,  Alcohol  6x2  testing.  However,  F  »0  data  (not  shown) 
was  completely  stable  and  the  300  psia  chamber  pressure  results 
indicated  lower  amplitude  standing  first  tangential  mode  instability. 

To  determine  the  exact  distributions  of  fuel  and  oxidizer 
present  in  the  regions  near  the  points  of  injection  (trom  both  a 
mass  and  mixture  ratio  standpoint)  and  thus  hopefully  explain  the 

positional  effects  observed  with  the  pulses,  the  6x2  injector  was 

22 

placed  in  the  spray  booth.  There,  using  at  first  a  liquid 

23 

collector  probe  similar  to  that  developed  by  Rupe  and  later 

switching  to  9  x  9  matrix  collector  (81-V  tubes  at  V'  center 
22 

distances  ),  the  necessary  samples  were  taken  and  analyzed.  The 
local  mixture  ratio  and  relative  mass  were  determined  using  RP-1  as 
the  fuel  and  appropriate  density  sugar  water  to  simulate  the  LOX. 

To  speed  the  initial  data  processing,  the  local  mixture 
ratio  and  relative  mass  were  computed  on  IBM  1620  equipment  using 
the  height  data  from  the  collected  columns  of  RP-1  and  simulated 
oxidizer  as  an  input.  These  data  were  then  used  in  conjunction 
with  adjacent  samplings  to  provide  an  overall  mass  and  mixture  ratio 
map.  Simplified  versions  of  such  maps  are  shown  in  Figures  IV-13 
and  -17  for  1.4  and  2.4  mixture  ratios.  Studying  these  data  one 
can  see  where  three  relatively  large  oxidizer  rich  zones  exist 
between  spuds  12-1,  4-5  and  8-9  on  the  outer  edge  of  the  injector 
under  both  1.4  and  2.4  mixture  ratio  operation.  Similar  broad  fuel- 
rich  zones  are  found  at  2-3,  6-7  and  10-11  locations  toward  the 
center  of  the  injector  (especially  visible  in  Figure  IV-15).  Whether 
any  effect  of  sectionalizing  the  injector  via  propellant-rich  zones 
is  responsible  for  the  improved  stability  (similar  to  the  wide  gaps 
previously  mentioned  or  the  effect  due  to  baffles)  remains  an 
interesting  interpretation.  On  this  question  it  should  be  noted  that 
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when  an  8  x  2  arrangement  va«  used  with  LOX,  Alcohol  the  results 
were  not  nearly  as  good  from  a  stability  standpoint  as  6  x  2. 

Also  of  interest  is  the  fact  that  the  LOX  momentum  is  so 
predominant  in  the  1.4  data  shown  in  Figure  IV- 15  and  -16  that  the 
mixture  ratio  averages  greater  than  two  above  the  center  of  the  spud 
never  reaching  the  1.4  mixture  ratio  at  which  the  injector  was 
operated. 

To  take  all  these  data  into  account  in  the  prediction  of 
spin  direction  is  no  easy  task.  This  approach  for  the  moment  has 
been  side  stepped  in  favor  of  simply  seeing  if  the  injector  is  really 
"balanced"  in  the  circumferential  directions.  With  this  in  mind,  a 
one -dimensional  view  of  the  mixture  ratio  and  relative  mass,  as  seen 
along  the  Injection  diameter,  is  illustrated  in  Figure  IV-16  for 
1.4r  and  in  Figure  IV-18  for  2.4r.  The  broader  LOX  and  RP-1  regions 
as  previously  mentioned  are  clearly  visible  in  these  figures.  The 
mixture  ratio  for  the  1.4r  data  at  the  spud  center  is  more  clearly 
defined  in  this  circumferential  s1j,'„<;  and  it  would  be  difficult  to 
justify  the  overall  mixture  ratio  that  was  present.  Referring  back 
to  Figure  IV-15,  apparently  the  fuel  was  displaced  inwardly  (possibly 
due  to  oxidizer  momentum  considerations). 

The  data  obtained  for  the  2.4  case  points  out  that  the 
mixture  ratio  and  relative  mass  fluctuations  are  considerably  less 
intense  than  the  1.4r.  This  is  very  obvious  when  one  compares 
Figure  IV-18  to  IV-16.  Both  samplings  were  done  at  three  inches 
from  the  injector  face.  Viewing  the  local  mixture  ratio  as  a  measure 
of  the  combustion  that  has  occurred  up  to  that  point,  one  would 
conclude  that  in  the  2.4  case  a  higher  percentage  of  the  total  com¬ 
bustion  has  taken  place.  This  improved  mixing,  whether  due  to  the 
relative  velocities  between  fuel  and  oxidizer  and/or  a  more  rapid 
blossoming  of  the  LOX  due  to  the  higher^  p,  has  a  direct  influence 
on  any  displacement  effects  present.  How  to  weigh  each  of  these 
effects  is  a  subject  currently  being  pondered.  For  the  moment,  however, 
let  us  consider  the  following  attempt. 
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Uslng  the  one-dimensional  approach  some  terms  must  first 
be  defined.  In  considering  available  oxygen  for  displacement,  the 
availability  could  be  considered  to  begin  at  stoichiometric  mixture 
ratio;  approximately  3. A  if  the  injector  is  evaluated  for  LOX,  RP-1 
and  2.1  for  LOX,  Alcohol.  Then  the  simple  relationship  to  use  is: 

M  -  M  r  rstoich. 

o  o  ■■■■■■■  ■ 

avail  total  .  , 


where  Mo  ..is  the  relative  mass  available,  Mn  ,  is  the 
'’avail  '’total 

total  relative  mass,  r  is  the  mixture  ratio. 


Adding  up  these  contributions  at  each  spud  location  with 
direction  being  assigned  via  the  oxidizer  to  fuel  preference,  it 
was  found  in  the  1.4r  that  a  difierence  in  "balance"  of  197.  existed 


between  the  two  directions.  At  2.4r  there  was  an  847.  unbalance. 


The  direction  of  predicted  spin  at  2.4  mixture  ratio  agreed  with 

Figures  IV-12  through  -14.  These  can  only  be  considered  preliminary 

results  until  additional  mixture  ratio  and  relative  mass  surveys 

can  be  completed.  Such  surveys  are  being  considered  as  a  Senior 

Thesis  topic  beginning  this  summer. 

Returning  to  the  discussion  of  the  rocket  firings,  the 

first  series  of  data  on  the  injector  elements  designed  for  2.3 

mixture  ratio  have  now  been  completed.  The  results  with  0  —  >  F 

testing  place  the  spinning  mode  instability  regime  between  2.1  and 

3,15  mixture  ratio.  Compared  to  Figure  IV-10,  this  indicates  a 

considerable  contraction  especially  since  the  45-20  grain  charge 

burst  d'fsk  combination  was  used  in  these  tests.  Based  on  equivalence 

ratio,  these  limits  represent  a  contraction  from  the  LOX,  Alcohol 

* 

instability  regimes.  Figure  IV- 1. 

Zero  inter-fan  spacing  (item  3  on  the  first  page  of  this 

section)  was  shown  to  oe  extremely  effective  in  controlling  high 

frequency  combustion  instability  for  LOX,  RP-1  as  was  the  case  for 

alcohol.  The  tests  run  to  date  include  only  the  1.4  mixture  ratio 
*In  contrast,  however,  T— *0  pulses  have  produced  IT  Instability 
as  with  the  1.4r  RP-1  spuds,  a  situation  not  found  for  LOX,  alcohol. 
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design  spuds.  The  one  caution  in  interpreting  these  tests  is  that 
the  c*  performance  is  down  because  in  essence  only  one  ring  of  fuel 
and  one  ring  of  oxidizer  are  provided.^ 

To  provide  controls  on  inter-spud  spacing  beyond  those 
obtainable  with  the  orientation  of  spuds  non-evenly  spaced  locations 
have  been  provided  in  one  injector  face.  Preliminary  tests  indicated 
a  marked  influence  of  pulse  gun  positioning.  For  this  reason  the 
tests  have  been  interrupted  while  the  injector  is  analyzed  in  the 
spray  booth. 

It  is  hoped  that  the  basic  experiment  being  conducted  on 
the  vapor  displacement  mechanism  (see  part  C  of  this  Section)  will 
provide  the  vital  data  needed  on  the  amplitude  and  phase  of  the  vapor 
and  droplet  movement.  Only  with  a  greater  knowledge  of  these  inter¬ 
relationships  can  Intelligent  stability  predictions  be  made  with 
regard  to  this  instability  mechanism. 
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Injector  Code: 


m 

a  x  (b  +  cB),  d  L  on  L  e  f  g 

a  -  Number  of  groups 

b  «  Number  of  spuds  in  group 

c  -  Number  of  blank  spaces  (B)  in  group 

d  ■  Inter- fan  spacing  between  the  like- on- like  (L  on  L)  spray  fans 

e  ■  Design  mixture  ratio  based  on  equal  velocity  of  75  ft/sec 

f  »  Orientation  of  center  line  formed  by  llke-on*llke  orifices: 
tangential  (T)  or  radial  (R) 

g  ■  Propellant  spray  fan  away  from  chamber  wall 

■  Propellant  spray  fan  nearest  the  chamber  wall  either  0  or  F 


Pulse  Arrangements: 


0-»F  F  0  F-»0  0  — 


g— a»g,  where  -  ■—  Indicates  pulse  direction. 


lUs 
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B.  LONGITUDINAL  NONLINEAR  ROCKET  STUDIES  USING  THE  SQUARE-MOTOR 

The  square-motor  approach  to  study  the  effects  associated 
with  high  density,  closely  spaced  injection  elements  was  outlined  in 
Reference  1.  Before  covering  the  most  recent  results  determined  in 
the  study,  it  would  be  well  to  briefly  review  the  previous  data. 

The  square-motor  as  its  name  implies  is  a  thrust  chamber 

made  of  three  and  four  inch  long,  copper,  square  chamber  sections. 

The  design  allows  chamber  length  variations  to  forty  inches  using 

a  hydraulic  jack  to  hold  and  seal  the  various  elements  (injector, 

group  of  chamber  sections  and  nozzle,  see  Figure  28,  Reference  1). 

The  various  injectors  under  test  provide  a  propellant-injection 

2 

density  greater  than  .5  lb.  per  in  -  a  considerable  step  above 

other  Princeton  rocket  hardware  that  flow  less  than  207.  of  that 

amount.  The  injector  evaluations  are  based  on  the  incidence  of  linear 

type  longitudinal  instability  which  can  be  partially  controlled  by 

2 

use  of  a  multior ifice  nozzle  to  vary  the  nozzle  admittance  or  by 
checking  the  incidence  of  nonlinear  longitudinal  oscillations  using 
injector-end  or  nozzle-end  pulsing  (gun  powder  and  burst  disk  com¬ 
bination). 

The  tests  previously  reported*  with  the  larger  injection 
orifices  (and  associated  higher  flow  rates),  the  pulsing  effects  and 
nozzle  changes  indicated  that  large  amplitude  (200  to  300  psi  peak- 
to-peak)  longitudinal-mode  oscillations  could  be  produced.  That 
only  the  smallest  pulse  was  necessary  (15  grain  -  1  K  disk)  to 
initiate  these  high  amplitude  oscillations  was  also  shown.  Both 
pulse  locations  could  produce  the  instability  although  the  nozzle 
location  generated  longitudinal  mode  oscillations  from  the  start 
(versus  mixed  mode  initiation  from  the  injector-end  pulses).  These 
were~simiIaP‘"lu  lh#-natural  rapid  instability  transition  observed 
with  this  chamber-injector  combination  (Figures  31  and  33  of  Reference 
1). 

The  tests  described  above  used  an  injector  which  could  be 
fabricated  in  four  simple  patterns  using  a  4  x  4  like-on-like  design. 
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For  clarity,  Figure  29  from  Reference  1  is  repeated  here  as  Figure 
IV-19.  Type  I  Injector  uses  aligned,  like  propellant  spray  fans 
that  end  Impinge  on  themselves  in  alternate  rows  of  fuel  and  oxidizer. 
Fuel  fans  end  Impinge  on  oxidizer  fans  in  the  unlike  end  impingement 
Type  II.  Type  III  Injector  allows  for  end  impingement  of  one  propel- 
lant  on  itself  with  side  Impingement  from  the  other  propellant. 

Finally,  Type  IV  rotates  the  fans  approximately  26°  to  avoid  end 
impingement . 

In  Reference  1  the  stability  behavior  of  Type  I  was  reported 
for  LOX,  Alcohol  as  producing  linear  longitudinal  combustion  insta¬ 
bility  at  high  and  low  mixture  ratios  at  lengths  greater  than  14  inches. 
Multiorif ice  nozzle  tests  extended  the  instability  to  cover  the  entire 
operating  range  of  mixture  ratio. 

The  Type  IV  Injector,  on  the  other  hand,  Indicated  linear 
Instability  at  only  very  low  mixture  ratio  at  a  length  of  38  Inches. 
Using  the  15-1K  pulses,  however,  the  range  of  instability  was  extended 
to  ten  inches.  The  Type  II  injector  exhibited  similar  stability 
behavior. 

The  Type  III  injector  design  was  then  tested  and  the  results 
are  illustrated  in  Figure  IV-20.  There  are  two  versions  of  the  Type 
III  design,  the  one  shown  in  Figure  IV-19  (where  the  LOX  end  impinges 
on  the  side  of  the  fuels  which  are  end- impinging  on  themselves)  and 
the  reversed  design.  The  former  is  referred  to  as  Type  III  pulsed 
0  ■■■»  0  and  the  data  in  Figure  IV-20  depict  the  high  amplitude 
oscillations  over  the  entire  operational  mixture  .  These 

oscillation;  exist  for  lengths  greater  then  tweW.  -cries  for  this 
design.  All  of  the  tests  shown  have  been  pulses  at  either  the 
nozzle  or  the  injector  end. 

The  other  Type  III  injector  (pulsed  F  -  » F)  yielded  longi¬ 

tudinal  instability  amplitudes  of  approximately  one  half  the  magnitude 
as  those  shown  in  Figure  IV-20  at  the  shorter  lengths  (at  39  inches 
the  amplitudes  were  similar).  The  lower  length  instability  limit 
extended  to  eleven  inches  between  L.  1  and  2.1  mixture  ratio.  Thus, 


neither  design  exhibited  improved  stability* over  Type  II  and  Type  IV. 

To  investigate  even  larger  orifice  sizes  on  this  program 
(the  injector  orifice  diameters  for  the  4x4  type  injectors  was 
.059"),  an  injector  design  with  even  fewer  injection  elements  was 
conceived.  Using  four-spud  elements  placed  on  1.33  inch  centers, 
each  spud  Injects  only  a  single  propellant.  The  orifice  size  was 
.120  Inches.  Since  spuds  were  used,  both  the  unlike  end  impinging 
characteristics  of  Type  II  and  the  rotated  feature  of  Type  IV  can 
be  provided  by  merely  rotating  the  spuds. 

Typical  data  from  such  tests  are  shown  in  Figure  IV-21  for 
the  Type  IV  injector  design  using  LOX-Alcohol.  Here  it  is  found 
that  without  pulsing,  instability  extends  down  to  at  least  twenty 
inches.  At  that  length  the  amplitudes  are  quite  low  (38  to  151  psi, 
often  intermittent)  while  with  the  Type  II  design,  the  instability 
range  was  similar  but  the  amplitudes  at  twenty  inches  were  still  more 
than  200  psi  peak-to-peak. 

As  in  the  case  of  the  transverse  mode  tests  (reported  in 
the  previous  part  of  this  section)  RP-1  replaced  the  alcohol  in  the 
test  program.  As  was  previously  discussed  in  such  a  propellant 
switch-over  (using  the  identical  Injector  design),  changes  will  occur 
in  the  relative  velocity  between  the  two  injected  propellants  and 
in  the  droplet  sizes  produced. 

The  Type  II  pattern  shows  the  least  change  from  the  results 
described  for  LOX,  Alcohol.  Lower  length  limits  are  comparable 
extending  to  approximately  nineteen  inches  as  seen  in  Figure  IV-22). 
The  upper  limit  of  instability  is  found  to  increase  in  mixture  ratio 
as  the  length  increases  to  30  inches.  A  change  of  this  type  was 
also  indicated  in  Figure  IV-21  for  the  same  equivalence  ratio. 

The  Type  IV  injector  design  stability  data  shown  in  Figure 
IV-23  also  Indicates  a  changing  stability  limit  from  a  mixture  ratio 
standpoint  as  the  length  increases.  The  lower  length  limit  in  these 
data  is  shifted  to  over  26  inches.  In  both  Figure  IV-22  and  -23 
pulsing  at  these  upper  mixture  ratio  limits  showed  only  minor  effects. 

*  Based  on  injector-end  pulsing  to  emphasize  effects  of  the 
injector  pattern  where  Type  II  and  IV  were  stable  below 
15  Inches. 
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Several  size  pulses  were  used  in  these  small  shifts  of  the  stability 
limits.  It  is  hoped  that  these  or  similar  tests  will  provide  a 
means  for  checking  experimentally  some  of  the  theoretical  longi¬ 
tudinal  nonlinear  models  covered  in  Section  III. 

With  the  LOX,  RP-1  combination,  the  high  velocity  flow  on 
the  LOX  side  restricted  the  high  mixture  ratio  operation  to  3.2r. 

As  a  first  step  the  orifices  on  two  spuds  were  Increased  to  .170 
inch  diameter.  Testing  with  the  Type  II  arrangement  and  the  large 
LOX  orifices  increased  the  limits  of  the  instability  regime  to  the 
1.8  to  3.8r  range  for  lengths  greater  than  23  Inches.  Intermittent 
longitudinal  oscillations  were  present  above  2.8r  at  nineteen  inches. 
With  the  Type  IV  design  and  large  LOX  holes,  high  amplitude  instability 
was  found  over  the  operational  mixture  ratio  range  at  23  inches  with 
Intermittent  traces  down  to  nineteen  Inches.  Thus,  the  addition  of 
the  large  LOX  orifices  and  unchanged  fuel  injection  increased  the 
instability  regime  for  both  injector  types  but  especially  in  the 
case  of  Type  IV. 

Reversing  the  propellants  so  that  the  large  orifices  were 
on  the  fuel  side  produced  the  following  effect:  Type  II  (large  F) 
was  stable  for  all  mixture  ratios  at  29%  and  38%  inches  with  damping 
taking  less  than  30  milliseconds  after  the  Initial  pulse.  Type  IV 
(large  F)  was  stable  for  all  mixture  ratios  at  22'j,  25V  and  38% 
inches  with  even  more  rapid  damping  than  lype  II.  One  can  see  from 
such  tests  that  the  increasing  of  the  sizes  of  the  fuel  droplets 
and/or  slowing  down  the  fuel  velocity  provided  a  considerable  im¬ 
provement  in  the  longitudinal  mode  stability  characteristics. 

Changes  on  the  oxidizer  side  via  larger,  more  slowly  moving  droplets 
increased  the  unstable  regime.  Injection  characteristics  of  both 
propellants  are  shown  to  be  capable  of  markedly  Influencing  the 
incidence  of  instability. 

Looking  at  c*  versus  length,  data  for  unstable  operation 
of  the  Type  IV  injector  with  .12  inch  diameter  LOX  orifices  as 
compared  to  .17  inch  diameter,  some  idea  of  the  effect  on  combustion 


distribution  is  seen.  From  figure  IV-24,  the  .17  inch  diameter  LOX 
orifice  tests  are  found  to  peak  at  a  length  of  approximately  25 
Inches  for  1.8  to  3.4r  while  c*  peaks  do  not  occur  for  the  .12  inch 
LOX  orifice  until  34  inches  for  1.8r  and  over  38  inches  for  2.4r. 

Even  this  later  mixture  ratio  trend  Indicates  that  as  the  LOX 
velocity  is  further  increased  the  peak  c*  occurs  at  a  longer  length. 

With  the  Type  II  injector  both  the  .12  and  .17  inch  diameter 
orifice  reached  peak  C*  near  25  inches.  The  Type  II  instability 
regimes  also  changed  only  slightly  with  orifice  size  compared  to 
the  major  length  changes  in  the  instability  regime  for  the  Type  IV 
configuration.  This  is  further  proof  of  the  close  connection 
between  the  axial  combustion  distribution  and  the  longitudinal 
stability  characteristics  of  a  rocket  motor.  Local  mixture  ratio 
conditions  at  barious  distances  from  the  injector  face  will  be 
studied  to  see  how  well  cold-flow  tests  can  predict  the  final  com¬ 
bustion  distributions. 

The  experimentally  observed  effect  of  combustion  distribution 
on  the  stability  is  presently  being  compared  to  the  predictions  of 
the  time-lag  theory  for  the  distributed  case. 

In  addition,  this  distributed-combustion  theory  is  being 
extended  to  Include  higher  Mach  number  effects  following  similar 
lines  to  the  analysis  of  Section  VI  which  involved  the  concentrated 
combustion  case. 
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C.  VAPOR  DISPLACEMENT  MECHANISM  - 
BASIC  EXPERIMENT  IN  THE  "PSEUDO"  ROCKET 

The  significance  of  vapor  and  liquid  droplet  displacement 

as  a  factor  in  the  initiation  and/or  in  the  sustaining  mechanism  of  the 

tangential  mode  of  high  frequency  combustion  instability  in  liquid  pro- 
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pellant  rocket  engines  has  been  questioned  for  some  time  .  As  indicated 
in  the  preceding  Yearly  Report  (Reference  1)  the  results  of  analytical 
studies  (Appendix  A  of  Reference  1)  as  well  as  experimental  tests  carried 
out  on  liquid  propellant  rocket  motors  at  Princeton  University  (also 
covered  in  an  earlier  part  of  this  Section)  have  shown  that  the  dis¬ 
placement  of  vaporizing  propellants  by  the  oscillating  gas  flow  in  an 
unstable  rocket  motor  is  an  important  factor  in  the  combustion  process. 
The  experimental  program  described  in  this  portion  of  the  report  is  an 
outgrowth  of  the  desire  to  study  this  effect  in  a  systematic  manner. 

Activity  on  the  project  during  its  initial  year  encompassed: 
the  development  of  a  15KV,  spark  lighted,  schl ieren- shadowgraph  system 
(wherein  both  the  liquid  and  gas  phases  could  be  readily  distinguished); 
the  optimization  of  test  fluids  to  be  used  in  the  experiments;  and 
the  analysis  of  vaporizing  liquids  in  quiescent  surroundings.  Details 
of  these  results  can  be  found  in  Technical  Report  647  (Ref. 24).  Re¬ 
search  activities  during  the  past  year  have  brought  the  experiment  to 
a  fully-developed  state  and  the  ability  to  obtain  significant  data  de¬ 
scribing  the  phenomena  is  now  at  hand. 

A  new  graduate  student  was  brought  into  the  program  last 
summer  during  the  period  when  efforts  were  being  expanded  in  the  de¬ 
sign  and  construction  of  a  "pseudo"  rocket  motor.  The  intent  of  the 
experiment  utilizing  this  apparatus  is  to  simulate  as  nearly  as  possible 
certain  significant  rocket  engine  parameters  present  during  the  tan¬ 
gential  mode  of  instability.  The  inherent  limitations  of  this  "cold" 
motor  approach  were  well  recognized  but  the  value  of  gaining  at  least 
partial  insight  into  the  effect  was  felt  to  greatly  outweigh  the  short¬ 
comings.  A  schematic  diagram  of  the  "pseudo"  motor,  is  shown  in  Figure  IV-25 
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The  unit  was  modeled  after  the  1500  pound  thrust  rocket  engines  used 
at  Princeton.  In  this  way  It  was  hoped  to  utilize  data  from  each  pro¬ 
gram  to  the  benefit  of  the  other. 

The  "pseudo"  rocket  chamber  Is  9  inches  in  diameter  and  6 
Inches  long.  A  photograph  of  the  apparatus  in  its  present  form  is 
shown  in  Figure  IV- 26.  The  Injector  end  of  the  chamber  contains  a  porous 
metal  plate  through  which  nitrogen  gas  may  be  passed  to  simulate  the 
through-flow  of  combustion  products.  The  vaporizing  liquid  Is  Injected 
via  a  single  jet  of  large  length  to  diameter  ratio.  This  simple  type 
of  injection,  providing  well-established  turbulent  flow  conditions,  25 
represents  the  first  injection  element  to  be  tested  and  hopefully  will 
be  followed  by  impinging  types  if  the  experiment  will  permit  it.  The 
chamber  side  walls  contain  two  diametrically  opposed  one  inch  thick 
optical  glass  windows  two  inches  wide  by  four  inches  long.  The  nozzle 
end  of  the  motor  consists  of  a  flat  plate  containing  evenly- spaced 
multiple  exhaust  nozzles  to  assist  in  the  attainment  of  uniform  flow 
through  the  chamber.  A  nitrogen  shock  tube  is  located  in  the  chamber 
wall  such  that  the  delivered  pressure  perturbation  enters  tangentially 
into  the  chamber.  The  chamber  has  been  designed  to  provide  positions 
for  as  many  as  five  Klstler  type  crystal  pressure  transducers  in  the 
exhaust  plate,  two  Klstler  or  Dynisco  pressure  transducers  in  the  side 
wall  as  well  as  thermocouples  and  an  axial  ly- ad justable  hot-wire  ane¬ 
mometer  entering  through  the  exhaust  plate. 

One  item  of  great  importance  in  this  experiment  is  the 
acoustic  environment  provided  by  chamber  dimensions  employed  and  the 
use  of  nitrogen  gas  at  somewhat  less  than  ambient  temperature  (-'-528*R). 
Thus,  the  first  tangential  mode  of  this  cylindrical  chamber  is  approx¬ 
imately  880  cycles  per  second.  The  first  longitudinal  mode  would  ap¬ 
pear  at  1130  cps  while  the  1L,  IT  combined  mode  and  the  second  tangen- 
ital  modes  would  be  predicted  for  just  over  1400  cps. 

As  might  be  expected,  various  problems  arose  in  adapting 
the  original  optical  system  to  the  present  apparatus  and  indeed  in  per¬ 
fecting  the  apparatus  itself.  The  placement  and  focusing  of  the  optical 
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components  proved  to  be  very  critical  since  small  changes  in  the  position 
of  the  spark  source  and/or  mirrors  greatly  decreased  the  clarity  and  con¬ 
trast  of  the  photographs.  Alignment  of  the  windows  also  was  quite  criti¬ 
cal.  These  windows,  of  course  were  designed  to  be  perpendicular  to  the 
collimated  light  beam  within  close  tolerances  to  avoid  distortion. 

First,  the  optical  system  was  focused  without  the  windows. 

When  the  initially  procured  optically  ground  pyrex  glass  windows  were 
inserted  it  was  found  that,  due  to  internal  stresses  and  density  gradients 
in  the  glass,  the  liquid  droplets  and  vapor  were  partially  obscured. 

These  windows  were  replaced  by  optically  ground  aid  polished  plate  glass 
windows.  With  these  new  windows  the  fine  detail  of  the  jet  could  be 

observed  and  the  liquid  stream,  the  droplets  and  v^por  could  be  dis- 
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tinguished  as  in  the  previous  research.  Increased  contrast  attained 

by  varying  filter  and  film  combinations  was  also  partially  responsible 

for  the  improved  photographic  data. 

All  the  above  adjustments  were  made  using  either  helium  or 

F»*eon  21  as  the  simulated  propellant  with  the  chamber  at  atmospheric 

pressure.  It  had  previously  been  determined  that  the  normal  boiling 

point  of  Freon  21  was  such  that  it  vaporized  readily  thus  meeting  the 

specification  deemed  desirable  for  simulation  of  the  actual  propellant. 

During  the  de  elopment  of  the  equipment,  however,  as  the  chamber 

pressure  was  Increased  (by  blocking  off  some  of  the  small  nozzles  which 

make  up  the  exit  plate  as  well  as  increasing  the  flow  of  carrier  gas) 

it  was  found  that  the  Freon  21  no  longer  vaporized  sufficiently  due  to 

its  higher  boiling  temperature  at  these  pressures.  As  a  consequence  at 

increased  chamber  pressure  different  Freons  (Freon  14  for  Pc  75  pslg 

and  Freon  12  for  P  125  psig)  will  be  required  to  adequately  vaporize 

c 

and  allow  observation  of  the  phenomena. 

Concurrent  with  the  development  of  the  apparatus  a  compara¬ 
tive  analysis  of  gas  densities  in  actual  rocket  engines  and  those 
occurring  in  the  experimental  apparatus  was  made  which  indicated  that 
a  pressure  of  14.7  psia  in  the  apparatus  was  equivalent  to  230  psia  in 
a  typical  L0X/RP- 1  rocket  engine  or  360  psia  in  a  LOX/LH^  engine.  The 
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controlllng  relationship  is  simply: 
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where  P  is  pressure,  M  is  the  molecular  weight  and  T  the  temperature. 

The  subscript  p  is  for  the  pseudo  rocket  chamber  and  c  actual  com¬ 
bustion  chamber  being  modeled.  The  experimental  chamber  is  designed 
for  use  up  to  300  psig.  Thus  from  a  density  standpoint,  the  apparatus 
can  simulate  a  wide  range  of  chamber  pressures  which  are  currently  or 
practical  interest,  and  those  which  may  be  of  Interest  if  extremely 
high  chamber  pressures  (3000  psia  or  greater)  are  utilized. 

Although  the  problems  just  discussed  were  time  consuming, 
the  primary  problem  encountered  with  the  apparatus  was  involved  in  the 
operation  of  the  nitrogen  shock  tube.  In  its  original  form,  a  solenoid 
was  designed  to  drive  a  pointed  steel  dart  through  a  burst  disk  thereby 
initiating  disk  rupture  and  releasing  the  nitrogen  pulse.  In  actual 
practice  it  was  found  that  the  solenoid  did  not  have  sufficient  power 

to  break  the  disks  in  a  consistent  manner,  hence  the  release  of  the 

pulse  was  not  reproducible.  As  a  result,  the  shock  tube  was  redesigned 
so  that  the  dart  is  now  forced  through  the  burst  disk  by  energy  from 

an  air  cylinder.  The  latched  air  cylinder  shaft  is  released  by  the 

original  solenoid.  (See  Figure  IV- 26). 

Utilizing  previous  work  accomplished  under  Grant  NsG-99-60 
(Ref.20)the  gun  tube  was  designed  to  contain  a  volume  of  approximately 
1.4  cubic  inches  (3-3/4"  long  of  1/2"  tube).  Tests  were  run  using 
this  tube  with  1700  psig  in  the  gun  and  with  chamber  pressure  varying 
from  0-300  psig.  The  oscillating  pressure  in  the  chamber  was  measured 
by  Klstler  pressure  transducers  (Model  #601)  at  various  radial  positions 
in  the  exhaust  plate.  The  A.  C.  component  of  the  output  of  the  trans¬ 
ducers  was  recorded  on  tape  and  the  traces  were  then  printed  at  80 
lnches/sec  on  an  oscillograph.  Those  conditions  of  pressure  and  volume 
in  the  pulse  gun  produced  initial  amplitudes  of  the  pulse  in  the  order 
of  only  3  psl  peak- to-peak.  The  pressure  traces  Indicated  that  the 
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oscillating  portion  of  the  pulse  decayed  in  a  few  milliseconds.  In 
order  to  determine  relative  pulse  strengths,  to  check  the  instrumentation, 
as  well  as  compare  data  with  the  results  of  previous  work  (Ref. 20  )» 
several  tests  were  made  using  gunpowder  charges  instead  of  the  nitrogen 
gas.  In  this  way  peak-to-peak  pressure  of  approximately  100  psi  were 
attained.  These  tests  agreed  quite  well  with  the  earlier  investigation 
producing  frequencies  of  880  cps  as  expected.  It  was  known  that  the 
strength  of  the  oscillating  pressure  component  and  its  duration  were 
functions  of  pressure  level  of  the  initial  pulse  and  the  quantity  of 
gas  injected  into  the  chamber.  Therefore  the  tube  length  was  increased 
to  3  times,  and  finally  to  9  times  its  original  length  with  the  pres¬ 
sures  increasing  from  1000  to  3000  pslg.  Upon  studying  the  pulse  traces 
it  was  found  that  many  overtones  were  present  and  the  frequencies  did 
not  agree  too  well  with  the  theoretical  frequencies.  In  addition,  the 
peaks  of  the  pulse  were  much  sharper  than  expected. 

Upon  close  observation  of  the  apparatus  it  was  found  that  the 
porous  injector  plate  had  buckled,  apparently  as  a  result  of  the  powder 
charge  over  pressure.  When  the  tests  were  repeated  after  the  porous 
plate  had  been  repaired,  the  frequencies  agreed  well  with  the  theory 
and  the  peaks  of  the  wave  were  noticeably  smoother.  Utilizing  a  volume 
nine  times  the  original,  together  with  3000  psig  nitrogen,  a  peak-to- 
peak  amplitude  of  15  psi  was  obtained.  In  addition,  it  was  determined 
that  with  the  high  amplitude  pulses  the  overtones  died  out  more  quick¬ 
ly  than  the  primary  frequency  and  the  pulse  shape  became  more  like  a 
sine  wave  with  peak-to-peak  amplitudes  of  3-5  psi.  Band  pass  filters 
wore  installed  in  the  transducer  output  circuits  to  remove  extraneous 
noise  and  harmonics.  It  had  been  anticipated  that  peak-to-peak  oscil¬ 
lations  of  the  order  found  with  the  highest  pressure  pulse  would  be  re¬ 
quired  to  observe  significant  droplet  and  vapor  movement  in  the  jet. 

At  this  point  in  the  development  of  the  apparatus  a  variable 
delay  trigger  circuit  was  installed.  In  operation  this  circuit  is 
activated  by  the  initial  shock  wave  passing  over  one  of  the  pressure 
transducers  in  the  exhaust  plate.  The  transducer  output  then  initiates 
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the  15K  spark  source  at  desired  delay  times  between  0  mllllsec  and 
40  mllllsec  thus  producing  a  photograph  on  the  film  plate.  The  circuit 
proved  to  be  accurate  and  reproducible  to  within  3  tenths  of  a  milli¬ 
second  . 

Some  difficulties  were  also  encountered  with  the  apparatus 
In  attaining  reproducible  discharge  of  the  nitrogen  gas  with  the  shock 
tube  due  to  uneven  rupturing  of  the  burst  disks.  This  was  corrected 
by  redesigning  the  cutting  head  on  the  steel  dart  which  punctures  the 
d laphragm. 

At  this  point  the  entire  apparatus  was  considered  operation¬ 
al  and  specific  experiments  were  undertaken.  It  was  concluded  that 
the '^seudd'  motor  would  be  operated  initially  at  approximately  10  psig 
since  this  was  equivalent  to  400  psla  representing  a  reasonable  cham¬ 
ber  pressure  level.  This  chamber  pressure  was,  of  course,  somewhat 
less  than  that  required  to  assure  sonic  flow  across  the  exhaust  noz¬ 
zles,  but  It  was  felt  that  this  condition  would  not  effect  the  signif¬ 
icance  of  the  experiments.  Experiments  at  higher  pressure  will  be 
made  in  the  future  to  determine  the  validity  of  this  assumption. 

In  practice,  a  run  consists  of  setting  up  the  15KV  spark 
source,  darklng  the  room,  opening  the  film  shutter,  setting  the  de¬ 
lay  circuit  for  the  desired  time  period,  establishing  the  flow  of 
vaporizing  liquid,  establishing  the  nitrogen  thorough- flow,  energizing 
the  delay  circuit,  and  closing  the  shutter.  Each  test  sequence  takes 
approximately  30  seconds.  The  optimum  use  of  the  equipment  requires 
that  a  number  of  photographs  be  taken  before  they  can  be  conveniently 
developed  and  analyzed  as  a  series. 

Analysis  of  initial  test  data  indicated  that  with  1200  psig 
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In  the  shock  tube  and  with  the  smallest  volume  in  use  (1.4  inches  ), 
the  vaporizing  Freon  jet  was  blow  completely  out  of  the  viewing  win¬ 
dow  until  long  after  the  pressure  trace  showed  that  the  pulse  had  de¬ 
cayed.  This  result  was  somewhat  unexpected  since,  as  indicated  pre¬ 
viously,  it  was  felt  that  larger  perturbations  would  be  required  to 
disturb  the  jet  to  this  extent.  It  was  concluded  that  the  continued 


unidirectional  disturbance  of  the  Jet  was  caused  by  the  nitrogen  flow 
from  the  shock  gun  "choking"  at  exit  port  thus  creating  a  steady  ro¬ 
tational  flow  component  In  the  chamber.  It  was  also  evident  that  a 
severely  turbulent  zone  of  the  injected  nitrogen  gas  existed  in  the 
upper  portion  end  of  the  chamber.  The  specific  knowledge  that  this 
turbulence  was  caused  by  the  pulse  gun  was  confirmed  In  a  series  of 
auxiliary  experiments  where  pulses  were  Introduced  without  the  use  of 
liquid  injection  or  the  nitrogen  through  flow.  In  the  experimental 
set  up  the  pulse  gun  had  been  mounted  so  as  to  fire  directly  at  the 
Jet  approximately  one  inch  below  the  point  of  liquid  Injection.  The 
Interaction  of  the  rotational  flow  component  as  well  as  the  photographic 
Interference  of  the  turbulence  made  it  desirable  to  change  the  position 
of  the  gun  so  that  it  would  be  Introduced  on  the  opposite  side  of  the 
chamber  from  the  injector  near  the  exhaust  end. 

A  minature  probe  with  a  pressure  transducer  on  the  tip  was 
placed  Immediately  in  front  of  the  gun  port  to  determine  the  gun  empty¬ 
ing  time.  The  flow  from  the  gun  Indicated  that  a  "choked"  flow  con¬ 
dition  was  occurlng  and  that  the  time  to  decay  was  greater  than  de¬ 
sired  In  the  experiment.  It  was  concluded  that  while  this  effect  was 
undesirable,  that  by  reducing  the  pulse  strength  it  could  be  distin- 
quished  from  the  effects  of  the  oscillating  pressure  component  and 
thus  need  not  curtail  the  experimental  program.  An  Improved  technique 
using  a  piston  in  the  pulse  gun  so  as  to  discharge  all  of  the  flow  In¬ 
stantaneously  Is  being  investigated.  Several  series  of  tests  were 
made  on  the  system  in  the  foregoing  configuration  wherein  only  pressure 
transducers  and  thermocouples  were  employed.  These  tests  Indicated 
that  indeed  vapor  and  droplets  from  the  Jet  were  moved  in  some  inter¬ 
esting  patterns  but  the  movement  was  not  readily  capable  of  being  cor¬ 
related  with  the  pressure  traces  from  the  Kistler  pressure  transducers 
in  the  exhaust  plate.  It  was  felt  that  much  of  the  difficulty  arose 
from  the  general  turbulence  level  present  in  the  bottom  of  the  chamber 
as  a  result  of  the  shock  pulse.  Because  of  this  problem  as  well  as  a 
desire  to  obtain  velocity  data  a  hot  wire  anemometer  was  added  to  the 
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chamber  and  a  pressure  transducer  was  placed  in  the  chamber  side  wall. 
Typical  of  the  data  obtained  with  the  current  system  is  the  following: 

Chamber  Pressure  10  pslg 

Injector  Pressure  Drop  124.6-131.6  pslg 

Pulse  Gun  Pressure  1200  pslg 

Burst  Disks  1100  Aluminum  (0.016"  thick) 

Freon  Injector  Temp.  69°F 

Chamber  Temp.  61°-58°F 

Shock  Tube  1.4  In."* 

A  sequence  of  ten  photographs  were  taken  at  varying  Intervals  out  to 
approximately  5  milliseconds. 

The  oscillographic  data  obtained  from  the  test  sequence  Is 
shown  In  Figure  IV- 27.  The  figure  clearly  shows  the  poor  quality  of 
the  transient  pressure  record  when  the  transducer  Is  located  In  the 
exhaust  plate  and  the  Improvement  In  quality  obtained  by  moving  the 
pressure  transducer  closer  to  the  Injector  end  on  the  chamber  wall. 

The  traces  have  been  filtered  for  all  frequencies  below  400  cps  and 
above  2000  cps.  The  misalignment  of  the  phasing  of  the  traces  on  the 
record  is  due  to  the  fact  that  the  various  probes  are  located  at  dif¬ 
ferent  circumferential  positions  in  the  chamber.  The  timing  marks 
shown  In  the  figure  refer  to  times  from  the  initial  shock  wave  passing 
over  the  pressure  transducer  located  in  the  exhaust  plate.  The  upper 
trace  Is  a  velocity  measurements  using  hot-wire  anemometer  instmaen- 
tatlon.  Both  the  pressure  and  velocity  data  clearly  indicate  ap¬ 
proximately  850  cps  fluctuations  close  to  the  characteristic  frequency 
calculated  for  the  pulsing.  These  preliminary  data,  indicate  that  an 
inphase  relationship  exists.  This  phase  relationship  is  characteristic 
of  the  spinning  mode  and  Indicates  that  a  perturbed  state  consistent 
with  that  seen  In  an  actual  unstable  combustion  chamber  has  been  at¬ 
tained  . 

In  Figure  IV- 28  the  shock  wave  has  passed  the  far  viewing 
window  and  has  not  yet  passed  the  liquid  jet  Itself.  The  structure 
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of  the  Jet  is  typical  of  those  studied  earlier  under  quiescent  con¬ 
ditions.  The  opaque  positions  define  the  liquid  phase  and  the  cellular 
structure  defines  the  Freon  vapor.  Note  that  the  Jet  is  vaporizing  in 
a  uniform  manner  and  that  some  small  droplets  can  be  seen  moving  down 
and  away  from  the  core.  The  turbulence  from  the  nitrogen  pulse  is 
clearly  vlsable  in  the  lower  right  hand  corner  as  a  mottled  zone. 

Figure  IV-29  was  taken  approximately  0.8  of  a  millisecond 
later.  The  shock  wave  has  passed  the  jet  with  negligible  effect.  Note 
a  slight  indication  of  the  movement  of  vapor  to  the  right  side  of  the 
Jet  in  the  unobscured  upper  portion  of  the  photograph.  This  movement 
is  in  the  direction  of  the  pulse  (i.e.  from  left  to  right  on  the  near 
side  of  the  chamber).  The  liquid  portion  of  the  jet  appears  to  be  un¬ 
affected.  At  1.2  mllllsec  later  when  the  disturbance  should  have  made 
approximately  2  revolutions  in  the  chamber  (1^  in  relation  to  the  jet) 
the  jet  appears  as  shown  in  Figure  IV- 30.  Here  the  liquid  portion  of 
the  jet  has  been  broadened  and  many  droplets  appear  to  be  at  the  Jet  sur¬ 
face.  The  vapor  is  clearly  uniformally  displaced  to  the  right  side  of 
the  Jet  and  is  spread  out  over  a  larger  area  than  before. 

As  shown  in  Figure  IV-31  another  0.8  mllllsec  later,  the  jet 
was  spread  out  considerably.  Here  the  first  vestiges  of  a  sine  wave¬ 
like  vapor  sheath  begins  to  appear  along  the  axial  path.  In  Figure 
IV-32  (from  3.9  millisec)  the  jet  has  taken  on  even  more  obvious  dis¬ 
placement  variations  along  the  axial  dimension.  The  Jet  appears  to  be 
shedding  fine  droplets  and  vapor  at  the  peaks  of  the  waves  in  al¬ 
ternately  opposite  directions.  The  jet  also  appears  to  have  formed 
a  shape  of  approximately  three  cycles  as  it  flows  through  the  chamber. 
Theoretically  at  the  frequency  of  the  fundamental  tangential  mode 
there  should  be  2.6  cycles  in  this  length.  Here  the  dispersion  of 
the  jet  is  quite  noticeable  and  in  places  one  can  see  completely 
through  what  was  the  liquid  jet  core.  In  Figure  IV-33  the  jet  is 
beginning  to  drift  into  a  much  broader  pattern  and  seems  to  be  more 
of  a  unidirectional  behavior  (time  A. 8  ms.). 

The  results  of  the  experiment  to  date  have  established 
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more  clearly  a  picture  of  what  physically  happens  to  vaporizing 
liquid  streams  under  a  transverse  perturbation.  Specifically,  in 
resume,  the  result  to  data  show: 

1.  Transverse  pulses  can  completely  destroy  a  vaporizing 
jet . 

2.  Very  small  perturbations  can  readily  effect  the  position 
and  shape  of  the  jet.  This  may  lend  some  credence  to 
the  thought  that  linear  or  non-linear  instability  may 

be  one  continuous  phenomena  and  that  the  effect  can 
readily  grow  from  small  perturbations. 

3.  The  vapor  and  fine  droplets  surrounding  the  jet  while 
under  the  Influence  of  the  oscillating  field  take  on 
a  cyclic  axial  shape  at  a  frequency  approximating 
that  of  the  fundamental  mode  of  the  perturbation. 

4.  In  the  range  of  oscillation  amplitudes  tested,  the 
spray  and  principal  portions  of  the  vapor  produced  by 
the  oscillating  transverse  flow  appear  to  remain  in 

a  position  local  to  the  jet.  The  rotational  flow 
component  tends,  after  some  time,  to  displace 
significantly  the  droplets  and  vapor. 

During  the  course  of  the  experiment  it  is  intended  that  an  impinging 
jet  will  be  installed  in  place  of  the  single  jet  now  in  use  in  an 
attempt  to  observe  the  behavior  of  a  range  of  injection  elements. 

In  addition,  it  is  felt  that  a  much  better  insight  into  the  role 
that  baffles  play  in  inhibiting  instability  can  be  obtained  by  in¬ 
corporating  various  baffle  configuration  and  observing  the  behavior 
of  the  vaporizing  jet  both  in  the  baffle  cavities  and  beyond  the 
baffle  tips.  Such  studies  could  vastly  increase  the  knowledge  of 
how  a  baffle  accomplishes  it  stabilizing  functions,  Appendix  B 
discusses  one  approach  to  simplified  baffle  design. 


D.  BASIC  EXPERIMENTS  IN  DROPLET  RESONANT  RIFNOMENA 


The  research  conducted  on  this  particular  project  is 
undertaken  to  determine  what,  if  any,  mechanism  exists  between 
droplet  velocities,  droplet  size  distribution  and  pressure  oscil¬ 
lations  in  a  "pseudo1'  rocket  environment.  This  environment  is 
created  in  a  six  inch  diameter  cylindrical  chamber  of  variable 
length  (see  Fig.  IV- 34).  The  injector  being  studied  is  mounted  at 
one  end  of  the  chamber.  A  series  of  windows  positioned  parallel  to 
the  spray  axis  are  used  to  observe  the  first  4-1/4"  of  the  spray. 

If  desired,  the  windows  may  be  moved  to  greater  axial  distances 
from  the  injector  face.  The  chamber  is  pressurized  through  radial 
vents  around  the  windows.  Nitrogen  flowing  through  these  vents 
keeps  the  windows  free  of  fine  droplet  spray.  The  nitrogen  exhausts 
through  a  variable-speed  siren  at  the  opposite  end  <r  the  chamber, 
thus  providing  the  desired  pressure  oscillations.  These  oscillations 
simulate  certain  characteristic  fluctuating  pressure  conditions 
found  in  rocket  thrust  chambers  experiencing  combustion  instability. 
The  longitudinal  mode  of  pressure  variation  is  normally  produced  in 
the  chamber,  although  other  modes  can  be  studied. 

The  chamber  pressure  limit  is  450  psi  with  the  plexiglass 
section;  the  resonant  frequencies  of  current  interest  are  500  and 
1000  cps  (first  and  second  longitudinal  modes);  the  injector  pressure 
drop  has  a  range  of  0-400  psi;  peak-to-peak  pressure  changes  as  high 
as  20  psi  have  been  obtained  under  oscillating  conditions.  All  of 
these  values  are  in  the  range  of  the  rocket  engines  presently  being 
tested  under  this  contract.  Only  water  has  been  studied  to  date. 

One  phase  of  this  study  was  completed  last  fall  and  was 
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recently  published  as  a  technical  report.*”  That  work  evaluated  a 
diffractive  light  scattering  technique  which  predicts  the  Sauter 
Mean  Diameter  of  the  spray,  a  parameter  which  was  to  be  used  in 
observing  variations  under  oscillating  flow.  However,  the  inter¬ 
pretation  of  the  experimental  results  based  on  the  light  scattering 
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theory  was  found  to  deviate  considerably  from  theoretical  expectations. 
Results  from  theoretical  computations  indicate  that  these  deviations 
can  be  produced  by  a  particular  type  of  particle  distribution,  i.e., 
a  dominant  normal  distribution  of  particles  at  small  droplet  diameters 
followed  by  a  narrow  tail  of  droplets  extending  to  relatively  large 
droplet  sizes.  High  speed  photographic  spray  analysis  indicated 
that  this  form  of  droplet  distribution  does  exist  in  the  spray  fans 
of  the  impinging  jet  injectors  tested. 

It  was  shown  that,  in  addition,  for  impinging  jets  the 
correct  Sauter  Mean  Diameter  can  be  evaluated  within  the  theoretical 
limitations  if  specific  data  reduction  procedures  are  followed,  but 
that  the  Sauter  Mean  Diameter  in  itself  does  not  have  sufficient 
sensitivity  to  describe  the  spray  characteristics  adequately.  Had 
the  light  scattering  technique  fulfilled  all  the  experimental  re¬ 
quirements,  a  significant  simplification  would  have  been  possible 
in  the  analysis  of  the  droplet  variations  under  oscillation  as  well 
as  steady-state  operation.  Therefore,  a  somewhat  more  conventional 
approach  had  to  be  adopted:  photomicrography  and  individual  droplet 
analysis  which  would  provide  the  necessary  droplet  size  distribution 
data.  More  specifically,  the  method  used  is  required  to  stop  the 
motion  of  the  droplets  and  record  their  size  accurately,  i.e.,  to 
a  size  range  down  to  . 

Considering  the  camera  system  for  use  in  the  photomicro¬ 
graphic  approach,  a  number  of  requirements  were  taken  into  account. 
First,  it  is  well  known  that  any  lens  '■amera  has  a  depth  of  field 
which  is  dependant  upon  the  desired  circle  of  confusion.  In 
measuring  droplet  sizes,  however,  the  tolerable  circle  of  confusion 
is  a  function  of  the  droplet  size  (smaller  droplets  need  a  smaller 
absolute  error  in  depth  of  field  to  be  able  to  have  the  same 
relative  error  in  droplet  size  as  larger  droplets).  Thus,  a  drop 
size  distribution  obtained  from  a  camera  photograph  is  based  on  a 
count  of  droplets  in  a  volume  (the  area  seen  by  the  camera  times 
the  depth  of  field)  which  is  different  for  each  size  droplet.  In 
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comparing  the  data  obtained  in  this  research  to  earlier  data  from 
similar  tests,  it  was  necessary  to  take  this  correction  into  account. 

Since  the  primary  concern  of  this  experiment  is  in  dis¬ 
tinguishing  differences  in  spray  characteristics  between  steady-state 
and  oscillating  pressure,  consistent  errors  such  as  this  can  be 
ignored.  It  would,  of  course,  be  more  beneficial  to  be  able  to 
correct  for  these  errors  '0  that  the  data  taken  can  be  viewed  in  a 
proper  perspective  relative  to  other  similar  studies  and  be  compared 
to  other  works  in  the  future.  However,  since  the  evaluation  of 

steady-state  characteristics  of  sprays  (similar  to  those  here)  have 

21,26,  29 

been  undertaken  by  others,  no  attempt  will  be  made  to  repeat 

these  more  complete  studies. 

To  compensate  for  the  depth  of  field  effect,  the  following 
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relationship  from  Mack  and  Martin  is  useful: 
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where  the  depth  of  field  is  (As),  object  distance  is  s  ,  magnifi 
cation  is  M  ,  aperture  is  a  ,  and  tolerable  circle  of  confusion 
is  c  . 

For  values  of  c  encountered  in  these  tests,  the  second 
term  in  the  denominator  may  be  neglected  giving  a  linear  relation: 
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To  check  this  depth  of  field  Bredieldt  used  wires  of 
known  diameters  as  reference  targets.  By  controlled  camera  movement 
and  a  study  of  the  resulting  photographs  for  the  onset  of  out  of 
focus  condition,  the  depth  of  field  as  a  function  of  wire  size  was 
determined.  With  the  new  camera  arrangement  built  around  the  same 
telephoto  lens,  a  rechcck  of  the  equation  was  necessary.  It  was 
felt  that  targets  geometrically  and  optically  more  similar  to  the 
droplets  should  prove  even  more  useful.  A  microscope  slide  of  glass 
spheres  of  various  sizes  in  the  droplet  size  range  was  made.  Depth 
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of  field  tests  as  described  for  the  wires  were  conducted  using  this 
slide  at  the  precise  location  of  the  droplet  spray.  The  data  verified 
the  linear  relation  between  depth  of  field  and  circle  of  confusion. 

The  corrections  for  depth  of  field  which  are  mentioned  later  in  this 
report  are  obtained  by  dividing  the  population  of  each  droplet  size 
by  its  corresponding  depth  of  field.  The  resulting  drop  size 
distributions  represent  a  complete  distribution  for  all  the  droplets 
within  a  we  1 1 -determined  volume  (the  area  photographed  times  a  unit 
length  of  depth  of  field). 

Next  to  be  considered  was  the  choice  of  film  and  develop¬ 
ment  which  was  dependent  upon  photographically  satisfying  two 
conditions.  The  droplets  have  to  be  photographed  1)  under  lighting 
conditions  which  vary  with  different  optical  spray  depths,  and  2) 
with  sufficient  contrast  to  quickly  and  accurately  record  a  droplet 
image.  Since  films  with  high  contrast  have, in  general,  sma 1 1  ex¬ 
posure  latitudes,  these  two  conditions  cannot  be  satisfied  by  the 
choice  of  film  alone.  Therefore  tests  were  conducted  to  determine 
the  optium  film  and  developer  to  be  used.  The  final  procedure  uses 
Kodak  Plus-X  Portrait  Film  (because  of  it's  fine  grain  and  wide  ex¬ 
posure  latitude)  developed  in  Kodak  D-8  developer  providing  high 
contrast.  These  negatives  are  then  enlarged  to  a  final  magnifica¬ 
tion  of  about  30  on  "0"  grade  paper  (Kodak  Polycontrast  A,  without 
filter) . 

Also  to  be  considered  was  the  choice  of  the  final  value 
of  mangi f icat ion  of  the  camera  system  which  was  based  on  a  trade 
off  among  several  variables.  For  a  given  film  grain  size  as  the 
magnification  increases  smaller  and  smaller  objects  (in  this  case 
droplets)  can  be  recorded.  However,  there  are  three  limitations 
to  the  maximum  magnification  which  can  be  used. 

The  first  comes  from  the  depth  of  field  equation  discussed 
earlier.  That  equation  reveals  an  inverse  relationship  between 
depth  of  field  and  magnification.  Thus,  if  for  larger  magnifications 
the  depth  of  field  becomes  too  small,  the  number  of  photographs 


which  must  be  taken  to  obtain  a  large  enough  droplet  count  becomes 
excessively  large.  This  leads  to  cumbersome  and  time  consuming  data 
processing.  This  limitation  to  date  has  not  become  a  serious  problem. 

The  second  limitation  arises  from  the  need  to  maintain 
proper  exposure.  As  the  magnification  increases  so  does  the  distance 
between  the  lens  and  film,  and  therefore  the  exposure.  Since  the 
time  of  exposure  is  determined  by  the  speed  necessary  to  stop  the 
motion  of  the  droplets,  this  limitation  is  really  dependent  upon  the 
light  source  intensity.  With  sufficiently  fast  films  the  magnifi¬ 
cation  used  (five  times  size)  would  have  easily  fit  within  the 
exposure  limits.  However,  based  on  available  film  for  the  70  milli¬ 
meter  automatic  advance  camera  that  was  suitable  for  this  work, 
there  was  a  definite  magnification  limitation. 

Finally,  the  very  practical  limitation  arose  as  to  the 
overall  size  of  the  camera  system.  The  test  laboratory  which  houses 
the  entire  apparatus  is  of  relatively  limited  size.  This  laboratory 
geometry  was  certainly  one  factor  in  limiting  the  camera  to  five 
power  magnification  which  proved  to  be  adequate  in  resolving  the 
10  y*  droplets. 

With  these  considerations  in  mind  the  final  camera  system 
was  designed  and  checked  out  (see  Figure  IV-35).  The  lens  system 
consists  of  .1  200  mm  telephoto  unit  with  an  f  3.3  speed  rating.  The 
lens  system  is  mounted  at  one  end  of  a  four  foot  tube  with  the  film 
advance  apparatus  mounted  at  the  other  end.  Thus,  the  image  to  lens 
distance  's  fixed  resulting  in  a  constant  magnification  factor  at 
a  fixed  image  to  object  distance  (controlled  by  the  limited  depth 
of  field).  Movement  of  the  camera  system  is  remotely  controlled  via 
linear  actuators  to  provide  axial  and  transverse  positioning  of  the 
object  plane.  The  70  mm  film  is  remotely  advat ced  using  a  Beattie 
camera  back.  Thus,  a  series  of  photomicrographs  are  waken  in  rapid 
sequence  when  the  predetermined  operating  conditions  are  attained. 

A  typical  series  of  photographs  is  shown  in  Figure  IV-3fv  The  shutter 
remains  open  during  the  series  since  an  intermittent  light  source  is 
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used . 

To  satisfy  the  requirement*  for  the  light  source  portion  of 
the  overall  optical  system  it  is  necessary  to  provide  a  diffuse  light 
source  Reference  21.  Use  of  collimated  light  systems  would  lead  to  con¬ 
tributions  from  portions  of  the  spray  in  the  field  of  view  but  before 
or  after  the  depth  of  field.  The  diffuse  light  source  allows  for  even 
illumination  of  the  droplets  within  the  region  of  interest. 

Sufficient  light  is  another  requirement  for  properly  back¬ 
lighting  the  droplets  and  allowing  the  film  images  to  be  of  sufficient 
clarity.  Again  the  diffuse  light  source  can  more  easily  attain  this 
objective . 

Before  choosing  a  light  system,  a  method  of  photographing 

sprays  through  the  use  of  fluorescent  dyes  which  has  been  developed 

29 

at  the  University  of  Wisconsin  was  first  considered.  In  this 

technique  a  thin  sheet  of  high  intensity  light  is  focused  through  the 
spray  perpendicular  to  the  camera  axis.  The  light  causes  a  fluorescent 
dye  in  the  spray  droplets  to  radiate.  These  radiating  droplets  are 
then  the  only  ones  to  be  photographed.  Although  adoption  of  this 
technique  would  have  required  extensive  rebuilding  of  the  simulated 
combustion  chamber  used  in  the  tests,  it  did  have  one  desirable  ad¬ 
vantage.  By  narrowing  the  light  beam  to  the  camera  depth  of  field 
for  the  smallest  droplet,  every  droplet  in  the  light  beam  would  be 
photographed  in  focus.  The  correction  for  depth  of  field  would  then 
not  be  necessary. 

However,  although  Reference  29  did  use  this  fluorescent 
technique  to  obtain  spray  results  which  would  have  been  satisfactory 
for  our  purposes  it  was  felt  that  the  time  required  for  adoption  of  the 
chamber  would  not  be  worth  the  elimination  of  the  depth  of  field  vari- 
tion.  After  all,  as  mentioned  previously,  the  depth  of  field  correction 
can,  for  the  purposes  of  this  research,  be  considered  a  consistenc 
error  which  will  not  effect  the  differences  in  spray  characteristics 
which  are  sought. 

The  light  source  that  was  finally  decided  upon  uses  the  spark 
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from  a  15  kv  capacitor  discharge  of  about  V*  sec  duration.  The  light 
is  kept  in  a  diffused  beam  and  passes  through  the  spray  and  directly  into 
the  camera  lens.  This  back  lighting  of  the  droplets  combined  with  the 
photographic  processes  previously  discussed  provides  excellent  contrast 
as  can  be  seen  referring  again  to  Fig.  IV- 36.  The  duration  of  the  spark 
is  sufficient  short  to  stop  all  droplet  motion  in  the  data  taken  so 
far.  However,  at  higher  injector  pressure  differences  slight  blurring 
of  the  droplets  due  to  their  motion  may  be  encountered  especially  at 
locations  closer  to  the  injector  face. 

In  the  comparison  of  the  processed  data  to  earlier  work  the 
velocity  of  the  impinging  jets,  spray,  and  surrounding  gas  (N^)  were 
needed.  The  jet  velocities  were  obtained  by  calibrating  the  injector 
orifice  flows  for  a  given  pressure  drop  across  them.  The  gas  velocities 
are  determined  by  assuming  isentropic  choked  flow  through  the  outlets. 
This  assumes,  however ,  a  one-dimensional  flow  of  gas  which  in  the  region 
of  the  windows  is  not  quite  true.  As  mentioned  earlier,  all  of  the  N0 
is  injected  radially  into  the  chamber  through  vents  around  the  windows. 
The  velocity  of  the  spray  is  obtained  through  use  of  a  streak  camera. 

This  camera  will  also  be  used  in  studying  wave  patterns  of  droplets 
(see  Ref.  3  0)  close  to  the  injector  face  by  studying  the  frequency  of 
streaks  in  the  traces. 

The  first  series  of  tests  were  made  with  a  90°,  .040  inch  di¬ 
ameter  impinging  jet  injector.  It  is  similar  to  those  presently  used 
in  instability  tests  except  for  the  smaller  diameter  of  the  jets.  Fu¬ 
ture  tests  are  planned  using  significantly  larger  and  smaller  injector 

orifice  diameters.  The  present  injector  orifice  sire  was  chosen  based 

25 

on  previous  data  and  the  absence  of  window  fogging  problems. 

Before  any  attempt  was  made  to  determine  differences  in  dis¬ 
tributions  between  oscillating  pressure  and  steady-state  flow,  it  was 
necessary  to  obtain  steady- state  spray  characteristics  and  to  compare 
these  results  to  earlier  reported  evaluations  of  these  characteristics. 

(References  26  and  27.) 

The  wave  frequencies  (so  called  Christmas  tree  effect)  ob- 
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served  by  Bredfeldt  in  steady- state  tests  with  this  spray  injector  were 

27 

found  to  be  3100  cps.  These  agree  with  other  works  .  However, 

these  data  were  taken  for  sprays  in  quiescent  air  so  they  may  differ 
from  the  results  that  will  be  obtained  from  the  streak  camera  for  sprays 
in  the  simulated  combustion  chamber. 

The  processing  of  the  steady-state  data  for  number-size  dis¬ 
tributions  has  been  accomplished  by  comparing  the  individual  droplets 
on  the  30  times  magnification  photographs  to  a  plastic  template  of  cir¬ 
cles.  The  differences  in  successive  circle  diameters  correspond  to 
about  of  actual  droplet  size.  To  enable  these  data  to  be  cotm* 

pared  with  other  similar  data  for  impinging  jet  injectors  the  distribu¬ 
tions  had  to  be  grouped  to  represent  ranges  of  diameters  (Ad's)  in  which 
each  successive  range  was  in  a  ratio  of  \j  2  to  the  previous  range. 
Consideration  of  this  correction  along  with  the  depth  of  field  correc¬ 
tion  allows  the  distributions  obtained  to  be  plotted  four  different 
ways:  1)  uncorrected,  or  directly  as  obtained  from  the  photographs, 

2)  with  a  correction  for  the  depth  of  field,  3)  with  a  correction  for 
a  A  D  ratio  of  F  4)  with  both  the  depth  of  field  and  Ad  correc- 
t  ion . 


Data  obtained  3-3/4  inches  downstream  of  the  injector  for  a 
chamber  pressure  of  100  psig,  injector  pressure  drop  of  30  psi  (and 
a  corresponding  jet  velocity  of  52  ft/sec)  and  based  on  a  count  of  3032 
droplets  are  plotted  by  these  four  methods  of  analysis  in  Figure  IV-37. 
Bimodal  tendencies  as  reported  in  Reference  27  can  be  seen  here  in 
only  #1  and  #4.  This  again  is  in  agreement  with  the  literature,  Ref¬ 
erence  27  ,  The  data  taken  for  that  reference  were  from  a  collimated 
light  shadowgraph  technique.  Such  a  technique  has  no  depth  of  field 
effect.  The  data  were  processed  on  a  droplet  counting  machine  which 
had  a  Aw  D  ratio  of  F  The  conditions  correspond  to  the  #4 
plot  of  Fig.  IV-37.  The  data  from  an  earlier  work,  Reference  21,  which 
is  reevaluated  in  Reference  26  were  taken  using  with  a  photographic 
method  similar  to  that  used  here  and  processed  using  a  AD  ratio  of 
one.  No  mention  is  made  in  the  report  of  a  correction  for  depth  of 
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field.  These  data  correspond  to  the  #1  plot.  In  these  two  reference 

cases  bimodal  distributions  were  obtained  similar  to  those  obtained 

I  ere.  Furthermore,  the  variations  of  the  mean  values  of  the  modes 

27 

due  to  variations  in  the  injector  diameters  indicates  that  the 

mean  mode  values  obtained  in  this  experiment  are  ir  agreement  with  the 
tendency  indicated  by  these  earlier  works.  A  list  of  these  modes  ap¬ 
pears  in  Table  I,  Figure  IV-38. 

This  is  however  a  qualitative  comparison.  A  more  quantita¬ 
tive  comparison  with  previous  work  can  be  made  by  following  the  analysis 
of  Reference  21  .  In  that  work  the  drop  size  distribution  data  were 
analyzed  by  using  the  Nuk iyama-Tanasawa  expression: 
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where  R  is  the  volume  fraction  of  drops  having  diameters  <D  and  D 
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is  the  characteristic  volume  number  mean  droplet  diameter.  The 
values  obtained  from  tTIis  expression  were  then  used  to  formulate  the 
empirical  expression: 
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where  D  is  the  impinging  jet  diameter,  it's  velocity  and  AV 

the  velocity  differunce  between  the  background  gas  stream  and  the  liquid 


jet 


In  appl  ing  this  analysis  to  the  data  the  values  of  the  D 


30 

calculated  in  the  Nukiyama-Tanasawa  expression  do  not  agree  with  either 

3  .  2 

the  actual  values  (based  on  the  definition,  (D^q)  g  Sl  ND  ) 


IN 


or  the  values  needed  to  satisfy  the  empirical  expression.  Comparative 
values  are  listed  in  Table  II  ,  Figure  IV-38. 

Tliis  lack  of  agreement  in  the  data  is  not  understood.  It 
is  felt  that  when  a  more  complete  knowledge  of  the  stream  velocities 
ligament  breakup  region  and  final  droplet  size  distributions  is 
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obtained  from  photomicrograph  data  analysis  that  these  questions 
will  be  answered.  The  emphasis,  however,  will  be  on  determining  the 
changes  in  droplet  characteristics  between  a  steady-state  and  an 
oscillating  environment. 

Pressure  oscillation  data  was  taken  at  the  resonant  frequency 
of  500  cps  with  a  peak-to-peak  pressure  variation  of  approximately 
5  psi.  Preliminary  analysis  of  these  droplet  photomicrographs  in¬ 
dicate  a  tendency  for  definite  groupings  of  droplets  to  develop  under 
unsteady  conditions  (as  viewed  through  window  #2,  2-1/8  inch  point). 
These  groupings  are  more  pronounced  than  in  the  steady-state  photo¬ 
micrographs  and  possibly  a  phase  relationship  is  indicated  between 
the  droplets  and  the  oscillation  frequency.  This  preliminary  con¬ 
clusion  is  based  on  the  droplet  distributions  observed  under  oscil¬ 
lating  conditions  as  shown  in  Figure  IV-39.  Should  such  resonant 
effects  be  present,  then  coupling  between  the  combustion  process 
and  the  chamber  acoustics  could  result  in  combustion  instability. 

This  possibility  is  discussed  in  Reference  31  and  in  Section  V. 

Another  effect  also  appears  to  be  evident  when  the  droplet 
size  distribution  data  are  plotted  versus  similar  steady-state  infor¬ 
mation  (see  Figure  IV-40).  The  total  droplet  density  (number  of 
droplets  counted  per  picture)  is  considerably  lower  under  the  oscil¬ 
lating  conditions.  Assuming  that  the  same  mass  flow  is  being  provided 
by  the  injector  for  the  same  pressure  drop  across  it,  this  effect 
could  be  due  to  a  number  of  causes.  One,  the  water  droplets  may  be 
shattering  (or  even  evaporating)  into  particles  too  small  to  be 
photographed  (<10  JA  diameter)  prior  to  reaching  window  #2. 

Secondly,  the  opposite  case  may  occur,  i.e.,  the  ligaments  known  to 
form  in  such  sprays  may  be  extending  further  downstream  and  therefore 
would  contain  more  mass  than  they  do  in  steady-state.  However,  a 
visual  comparison  of  steady-state  and  unsteady  photomicrographs  in¬ 
dicates  a  reduced  number  of  ligaments  present  with  oscillating 
conditions.  Another  possibility  is  that  although  the  droplet  density 
may  be  lower  on  the  axis,  it  may  be  greater  out  in  the  radial  direction 
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indicating  a  transverse  broadening  of  the  spray  fan  in  terms  of  mass 
flux.  Preliminary  analysis  indicates  that  this  condition  is  present. 
Additional  photomicrographs  will  allow  these  possibilities  to  be 
fully  investigated.  It  is  hoped  that  the  results  from  the  streak 
camera  will  show  any  variations  in  droplet  wave  frequencies  and 
droplet  velocities  for  both  the  steady  and  oscillating  conditions. 

The  results  described  here  are  admittedly  from  a  very  small  sampling 
of  data.  The  visual  counting  of  droplets  is  a  relatively  slow 
process.  Many  photographs  have  been  taken  which  are  presently 
waiting  to  be  processed.  Now  that  the  techniques  of  obtaining  the 
data  have  been  established,  the  processing  and  analysis  of  the  data 
will  play  an  important  role  in  the  progress  of  the  research. 

It  is  planned  that  tests  covering  the  effect  of  variables 
such  as:  chamber  pressure,  frequency,  injector  pressure  drop, 
orifice  size,  transverse  and  radial  location,  etc.,  will  be  com¬ 
pleted  over  the  next  months. 


1 


V.  DROPLET  AND  DROPLET  WAKE  BURNING  -  MECHANISM 
FOR  COMBUSTION  INSTABILITY 


Introduction 

In  order  to  more  fully  understand  aid  control  the  occurrence 
of  combustion  instability  in  liquid  rocket  motors,  it  is  important 
that  the  responsible  mechanisms  be  well  understood.  Even  without  such 
information  the  lumped  parameter  approach  (utilizing  the  combustion 
parameters  rfr  and  n  from  the  pressure  sensitive  time  lag  approach 
of  Crocco ^  )  has  proven  extremely  useful  in  describing  instability 
regimes.  However,  without  a  detailed  knowledge  of  the  mechanism  in¬ 
volved  often  only  an  indirect  attack  on  the  problems  has  been  possible. 
For  example,  well-defined  limits  of  instability  could  be  used  to 

predict  other  stability  boundaries  or  instability  regimes  for  similar 

32 

rocket  motor  injector  designs.  However,  one  could  not  determine 
beforehand,  based  on  known  physical  or  chemical  characteristics, 
where  or  when  instability  would  occur  if  the  rocket  w.s  not  quite 
similar  to  those  previously  experimentally  investigated  (this  partic¬ 
ularly  pertains  to  the  injection  orifices). 

An  analytical  investigation  of  several  mechanisms*  has  been 
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carried  out  by  Strahle  and  Crocco.  The  mechanisms  under  consider¬ 
ation  were  (1)  the  oscillatory  burning  or  vaporization  of  liquid 
droplets  when  based  upon  quasi-steady  theories  of  vaporization,  (2) 
the  effect  of  droplet  displacement  from  mear  positions  caused  by 
chamber  gas  pressure  oscillations  at  the  injector  face. 

The  conclusions  reached  in  this  theoretical  study  may  be 
summarized  as  follows: 

There  were  strong  physical  and  mathematical  reasons  for 
rejecting  some  cotmonly  advanced  mechanisms  of  instability,  although 
these  mechanisms  could  be  a  contributing  factor.  The  evaluation  was 

*Other,  more  promising  mechanisms  were  also  analysed  and  are  summarized 
beginning  on  the  following  pages.  They  are  droplet  burning  and  droplet 
wake  burning. 
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based  on  whether  or  not  the  mechanism  could  account  for  the  initiation 
of  spontaneous  instability  resulting  fmm  small  perturbations  in  a 
rocket  engine. 

The  special  droplet  number  density  was  eliminated  as  a 
possible  mechanism  because  the  drag  was  estimated  to  be  too  small  to 
cause  the  necessary  relative  movement  during  the  cycle  time.  By 
this  same  reasoning,  positional  displacement  of  droplets  affecting 
the  positional  vaporization  rate  were  also  unable  to  contribute.  The 
vaporization  rate,  when  based  on  quasi-steady  theories,  was  judged 
inadequate.  This  was  based  on  response  to  a  pressure  wave,  response 
to  velocity  effects  when  the  droplets  were  close  to  a  velocity  node 
and  phase  problems  between  velocity  and  pressure.  Thus,  the  validity 
of  quasi-steady  vaporization  theory  was  questioned  when  used  in  con¬ 
junction  with  high  frequency  oscillations.  It  was  pointed  out  that 
a  valid  time-dependent  theory  of  the  vaporization  process  would  be 
necessary  to  explain  some  observed  instability  phenomena. 

Evaluation  of  the  pressure  oscillation  effects  on  the 
formation  of  droplets  near  the  injector  face  was  concluded  to  play  no 
part  unless  this  particular  mechanism  undergoes  a  strong  resonance 
in  the  frequency  range  of  interest  (see  Section  IV-D).  The  reason 
stated  was  that  since  combustion  occurs  in  a  region,  rather  than  at 


a  discrete  location  at  or  near  the  injector  face,  any  non-resonant 
effects  would  average  out  prior  to  combustion.  This  would  mean  that 
a  strong  mechanism  must  be  in  the  region  of  the  steady-state  combustion 

rather  than  only  at  the  injector  face. 
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In  References  the  analysis  of  instability  mechanisms, 

concentrating  on  droplet  and  droplet  wake  burning,  is  covered  in 
great  detail.  These  studies  pointed  out  that,  for  simple  droplet 
models,  frequencies  commensurate  with  the  diffusion  times  produced 
extremely  strong  response  in  vaporization  rate.  In  the  limit  of  very 
high  frequency  it  was  found  that  the  vaporization  rate  perturbation 
could  become  infinite  as  the  square  root  of  frequency,  but  that  the 


burning  rate  perturbations  (important  in  the  feedback  mechanism)  would 


remain  bounded. 

Using  a  flat  plate  model,  it  was  found  that  for  many  prac¬ 
tical  configurations  effects  of  periodic  liquid  heat-up  and  finite 
evaporation  kinetics  could  be  ignored. 

To  take  into  account  convective  effects  conditions  of  high 
free  stream  Reynolds  number  flow  were  considered  over  a  liquid  body. 
Only  for  sufficiently  low  Reynolds  numbers,  large  droplet  sizes,  high 
temperatures,  and  fast  reaction  kinetics  could  the  pressure  be  low 
enough  from  a  practical  point  of  view  and  still  retain  the  collapsed 
flame  assumption.  It  was  found  from  the  steady-state  theory  that, 
even  with  the  collapsed  flame,  a  significant  fraction  (order  of  1/2 
total)  of  the  vaporized  mass  would  be  convected  into  the  droplet 
wake  rather  than  burned  in  the  boundary  layer. 

The  burning  rate  perturbation  can  undergo  striking  changes 
with  frequency  for  the  assumed  wave  type  (longitudinal  mode)  if  one 
assumes  that  a  droplet  acts  partly  as  a  flat  plate  and  partly  as  a 
stagnation  point.  The  most  important  parameters  in  determining  this 
behavior  are  the  stoichiometric  heat  of  reaction,  q  ,  and  the  ratio 
of  the  free  stream  oxidizer  mass  fraction  to  the  stoichiometric  mass 
ratio  of  oxidizer  to  fuel,  Y  /j  •  The  velocity  effect  was  shown 
to  be  of  paramount  importance  to  this  behavior. 

The  theory,  in  many  ways,  only  roughly  approximated  the 

conditions  present  in  an  actual  combustor.  Nine  points  were  cited  in 

Technical  Report  671,  pages  151-152,  and  a  detailed  study  of  each 

would  be  required  to  determine  the  effects.  Cited  as  probably  the 

most  important  effect  was  that  a  great  deal  of  mass  would  be  carried 

into  the  wake  and  not  burned  in  the  leading  edge  boundary  layer; 

therefore,  it  was  concluded  that  the  problem  of  wake  flames  should  be 

treated  in  the  unsteady  state.  Even  with  such  limitations,  this 

droplet  burning  theory  explains  many  of  the  qualitative  features  of 
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longitudinal  rocket  engine  instability  if  the  burning  rate  pertur¬ 
bation  is  accepted  as  the  appropriate  feedback  function. 

In  a  paper  titled  "The  Unsteady  Laminar  Jet  Flame  at  Large 
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Frequencies  of  Oscillation  and  Application  to  the  Wake  of  a  Burning 

Droplet,"  Strahle  Introduces  the  important  case  of  wake  burning. 

37 

Below  Re  -  400  the  wake  should  be  laminar  and  stable  If  the 

flame  cannot  exist  from  the  forward  stagnation  point  and  the  droplet 

acts  as  a  flame  holder,  this  configuration  (i.e.,  droplet  burning  in 

a  free  convection  field  and  behaving  as  an  overventilated  laminar 
38 

diffusion  flame  )  can  still  be  expected  to  be  valid  but  with  the 
added  complication  of  an  unknown  ignition  point.  This  flameholder 
model  is  shown  in  Figure  V-l. 

Partial  premixing  of  the  fuel  and  oxidizer  components  would 
be  initiated  at  the  droplet  leading  edge.  The  leading  edge  would  be 
blown  by  vaporization;  therefore,  separation  would  occur  quite  early. 
The  wake  then  would  essentially  consist  of  three  parts,  the  re¬ 
circulation  region,  a  trailing  edge  vaporization  region  and  the 
burning  core  flow  region.  In  the  core,  conditions  would  be  fuel-rich; 
the  external  flow,  of  course,  would  be  oxidizer-rich.  Under  the 
proper  initial  conditions  near  the  droplet  (verified  from  steady- 
state  experiments)  modifications  to  the  controlling  relationships 
could  be  found  from  the  solutions  to  the  unsteady  laminar  jet  flame, 
and  should  hold  qualitatively  over  the  majority  of  the  burning  wake 
region.  It  is  this  droplet  wake  burning  theory  that  will  now  be 
applied  to  longitudinal  combustion  instability. 

Development 

One  result  of  the  unsteady  wake  theory  is  that  the  burning 
rate  perturbation  is  given  by 


-4-1  (1) 


^  ^  ^*S~'  fV  *e[dl 

Her®  \)/^  <P  is  the  local  burning  rate  perturbation  component 
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in  phase  with  the  local  pressure  perturbation,  x*  divided  by  the 
steady-state  local  burning  rate,  where  is  nondimensionalized  by  the 
stagnation  chamber  pressure.  This  ratio  is  the  important  quantity 
with  reference  to  engine  stability  and  must  be  approximately  one  or 
greater,  and  positive,  for  instability  to  occur.  %  is  the  ratio  of 
specific  heats.  q  is  the  stoichiometric  heat  of  reaction  divided  by 

'h  *k  — 

the  free  stream  enthalpy,  c^  T^.  is  the  wake  flame  radius  in  the 

steady-state  (a  function  of  axial  distance)  and  is  nondimensionalized 

by  the  droplet  radius.  It  therefore  has  a  range  of  approximately  one 

to  zero,  nip  is  the  steady-state  fuel  burning  rate  divided  by  2 

and  nondimensionalized  by  *(u£m  *  u*)  at  an  axial  location  for 

the  flame.  It  is  obtainable  only  by  a  solution  to  the  steady-state 

wake  problem  and  ranges  from  order  one  near  the  droplet  to  zero  at  the 

flame  tip.  T  is  the  temperature  distribution  in  the  wake  nondimens ion- 
★ 

alized  by  T  ^  and  is  a  function  of  axial  and  radial  distance.  It 
likewise  is  only  known  after  a  complete  steady-state  solution  to  the 
wake  is  carried  out.  is  the  flame  temperature  given  by 


(3) 


Y  /I  is  the  free  stream  oxidizer  mass  fraction  divided  by  the 
o 

stoichiometric  mass  ratio  of  oxidizer  to  fuel.  T  is  the  center-line 

c 

gas  temperature  near  the  ignition  point  in  the  wake.  Y^  is  the  fuel 

Q 

vapor  mass  fraction  at  the  same  axial  location.  r  is  radial  distance 
from  the  centerline.  Finally,  u)  is  the  dimensionless  frequency 


CJ  =  Co  *  .f  *  a.  *  *  (4) 

sQ 

This  high  frequency  solution  contains  no  terms  up  to  6  [  Zul 

of  the  "velocity  effect"  and  is  purely  a  "pressure  effect"  (see  Technical 
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Report  671  )  and  is  known  to  provide  good  qualitative  behavior  for 

an  analogous  stagnation  point  problem  down  to  U>  *  4  or  5 .  Such 
behavior  will  also  be  assumed  here  and  Equations  (1)  and  (2)  will  be 
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assumed  exact.  These  equations  hold  when  the  droplets  are  stationary 
with  respect  to  the  chamber.  It  is  claimed  here  from  unpublished  work 
that  to  0[l /CO]  it  is  also  valid  when  the  droplets  are  moving.  The 
reason  for  this  is  that  the  movement  causes  a  Doppler  effect  of  Q  [M] 
in  the  timewise  frequency  seen  by  the  droplets.  The  velocity  (or  in¬ 
compressible)  effect  results  being  ^[1]  are  therefore  not  affected 
and  still  die  out  as  1/to.  The  pressure  effect,  of  <5[M) ,  is  modified 
by  the  velocity  effect  times  the  extra  terms  introduced  by  this  Doppler 
effect  of  @[M]  .  This  creates  terms  of  &  ( M) ,  the  action  of  the 
Doppler  effect  on  a  pressure  effect  creating  terms  of  ©[M^]  which 
are  negligible.  Thus,  the  pressure  effect  is  only  altered  by  terms 
which  still  die  out  as  1/CO  which  are  neglected  in  Equations  (1)  and 
(2). 


Some  insight  into  the  lower  length  limit  of  longitudinal 

mode  combustion  instability  can  be  gained  from  Equations  (1)  -  (4). 

First,  from  Equation  (4)  we  can  calculate  the  frequency  for  such  a 

32 

lower  length  limit  based  on  Reference 

3000  cps  (approximately  8  inch 
cylindrical  length) 

For  the  .035  inch  diameter  injection  orifices  used  in  those  tests  we 

21 

cen  calculate  the  droplet  size  using  Ingebo's  data: 

a*  ^  3  x  10  ^  ft.  100^ 

The  chamber  density  at  600  psi  would  be: 

*  .2  lb/ft3 

From  Priem's  data  the  viscosity  would  be  approximately  equal  to 

■  5  x  103  lb/ft  -  sec 
Then  the  dimensionless  frequency  is 


CO  ~  6.8 


which  is  in  the  high  frequency  range.  Therefore,  the  equations  will 
be  assumed  valid  at  the  short  length  limit.  They  cannot  be  valid  at 
the  upper  length  since  the  functions  are  monotonic  in  frequency;  once 

*^F^  becomes  sufficient  to  cause  instability  (if  it  does  at  all) 
the  system  would  remain  unstable.  Supposedly,  the  neglected  terms 


I 


106- 


would  become  important  at  long  length.  It  now  remains  to  show  that 
Equations  (1)  and  (2)  do  indeed  predict  stability  below  a  certain 
length. 

Various  approximations  will  now  be  introduced.  Since  the 

2 

wake  is  axially  symmetric,  near  the  axis  T  goes  as  ar  +b  with  a 
and  b  constant.  If  this  functional  dependence  is  assumed  to  hold 
all  the  way  to 

T= 

where  5  is  axial  length  from  the  initial  point,  T  (0)  *  T  .  Now 

c  c 

the  integration  in  Equation  (2)  may  be  performed  with  the  result 


.  rf  f  ^ 

Now  assume  in  the  above,  that  T  is  linear  in 

c 

core  heat-up  with  axial  length 


Tt-TLOO  \(t)  j 

£  because  of  the 


T  (*)  *  (Tf  -  T  )€  +  T 

C  I  c  c 

Then  the  average  T  is  approximately  Tr/2  if  T  ^  1,  which  will 

c  t  c 

be  assumed.  Also,  assume  that  the  average  =  1/2  and 


Now  Equations  (1)  and  (2)  are  combined  to  yield 


%  - ^  % 

Concerning  chamber  instability  (  7nYtr^  <P  )  must  be 
a  certain  magnitude  to  overcome  the  nozzle  damping.^  Holding  the 
nozzle  type  and  chamber  pressure  fixed,  a  mixture  ratio  traverse 
produces  a  small  shifting  of  the  position  of  the  steady-state  flame 
front.  This  shifting,  however,  produces  only  a  very  small  change  in 


_ t  J 


-  t 
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the  required  feedback  unless  the  flame  front  is  very  near  a  pressure 

node.  Therefore,  this  magnitude  can  be  picked  for  all  mixture  ratios. 

It  is  analogous  to  twice  the  minimum  interaction  index  required  for 

instability.  Let  us  choose  (  K  /in  ^ 1  Then  there  results 

*f  *f 


Cj 


|  m'(  *-**-&>  1 
L  ^  +•  l Jn.  I  +r(6-0 j 


^  7 1  (5a) 

+r(6-0 

¥ 
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Now  it  is  known  from  Technical  Report  671  that  while  trends  are 
predicted  quite  low  in  the  frequency  range  by  Equation  (3),  absolute 
numbers  are  not.  Furthermore,  absolute  numbers  for  m^.  are  not 
available  without  a  complete  solution  to  the  wake  problem.  However, 
if  a  reference  point  is  known,  designated  by  a  subscript  zero, 
Equation  (5)  may  be  written 


CO  _ _ 

l  [  5.  +  2  JU+m-l)]/ 


and  it  may  be  possible  to  estimate  trends  for  m  /nu.  as  well  as 

*  f  t  f 

for  other  quantities.  Now  from  Equation  (4)  o 


where  L*  is  chamber  length  and  c* 
sound.  Then 


is  the  chamber  gas  speed  of 


This  states  the  lower  length  stability  limit  variation  if  one  point 
is  known.  It  remains  to  compute  these  quantities  for  an  actual 
system. 

★ 

In  a  rocket  chamber  the  ambient  temperature,  T^  ,  must 
increase  with  length  to  the  adiabatic  flame  temperature  for  the 


108- 


* 

particular  engine  mixture  ratio.  By  energy  considerations  must 

be  above  this  flame  temperature  initially  (since  burning  is  assumed 
at  stoichiometric  mixture  ratio)  and  decrease  to  it  in  time.  Thus 
on  the  average  T^>1.  Changes  in  the  adiabatic  chamber  temperature 
with  mixture  ratio  should  also  be  reflected  in  similar  changes  in 
T*^  .  It  will  be  assumed  that  T*^  /T*^  is  the  same  as  the  same 

ratio  for  the  adiabatic  flame  temperature.  This  quantity,  which 
appears  in  q/qQ  ,  can  therefore  be  calculated  from  equilibrium  cal¬ 
culations  concerning  the  chamber  gases  of  complete  combustion. 
Similarly,  c*/c£  an<l  j**/  9*  are  known,  is  assumed  pro¬ 
portional  to  Tjj^  so  that  is  known.  T^  affects  the 

results  only  slightly  if  it  is  low  enough  and  may  be  assumed  constant. 

Similarly,  T  is  primarily  determined  by  the  droplet  characteristics 
c 

and  may  be  assumed  constant.  This,  of  course,  is  open  to  question  if 
the  ignition  point  in  the  wake  seriously  changes  with  mixture  ratio. 
This  is  a  possibility.  It  will  be  assumed  that  Yq  is  proportional 
to  the  mixture  ratio  and  that  at,  say,  a  mixture  ratio  of  one 

(0/F)  Y  -  .2  on  the  average.  At  present  there  is  no  guide  to 

o*o 

i  ~ 

this  choice,  if  the  steady-state  distributions  of  oxidizer  are  not 
known.  For  the  alcohol,  oxygen  system  j  **  2 . 1  and  q^  &&  10.  Note 
q*  is  independent  of  T*  under  the  assumptions  of  the  wake  theory 

I 

so  that  only  T*^  in  the  denominator  affects  q  .  All  quantities 
are  now  available  for  the  computation  of  T^. 

The  computation  of  the  variation  of  m^.  is  extremely  dif¬ 
ficult.  If,  however,  the  flame  lies  on  the  outside  of  the  velocity 
inflection  point  in  the  mixing  region,  an  increase  in  Yq  will 
move  the  flame  toward  the  center  line  where  a  larger  velocity  gradient 
exists.  This  in  turn  will  increase  the  burning  rate.  It  will  be 
assumed  that  this  is  the  case  and  that  the  dependence  on  Yq  is 
the  same  as  that  for  the  droplet  theories.  , In  reality,  this  quantity 
depends  in  a  complex  manner  on  the  ignition  point  velocity  profile, 

Y  /J  ,  and  Y_  .  Exact  or  even  more  probable  behavior  of  this 
o  Fc 

function  can  only  be  determined  through  exoeriment.  At  any  rate  it 


\ 
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V*f  ■  )'  (7) 

o  o 

where  C  <v0.8  and  itself  varies  with  in  a  mild  manner. 

o  M 

The  final  quantity  for  determination  is  a*.  Now  strictly, 

a*  is  not  the  droplet  size  but  is  a  wake  thickness  and  will  vary 

depending  not  only  upon  droplet  size  but  upon  the  burning  and  convective 

environment.  Another  fundamental  assumption  is  that  a*  is  proportional 

to  droplet  size  alone.  This  is  to  be  interpreted  as  a  mean  size  of 

the  distribution  and  over  the  lifetime.  This  mean  will  be  assumed 

21 

proportional  to  the  Injected  mean.  Then  Ingebo's  correlation  suggests 


a* 


*  l-n*-n 
rj  r 


(8) 


where  n  is  between  1/2  and  1.  This  does  not  include  a  chamber 
density  or  pressure  effect  since  in  those  experiments  $  mo  *  or  p  * 
was  not  varied.  Such  an  effect  must  exist,  however,  and  the  most 
likely  candidate  is  a  P*  dependence.  Justification  for  what  as¬ 
sumption  should  be  made  can  be  found  in  unpublished  Princeton  data 
which  show  only  a  small  effect  on  the  short  length  limit  with  a 
pressure  Increase  as  long  as  the  pressure  is  sub-citical  for  the  fuel. 
Reardon's  Aerojet  compilations  show  a  very  small  increase  in  the 
length  with  a  pressure  Increase.  Equation  (6)  then  states  that 
a*  ,  where  p  m  .4  -  .5.  Thus,  p  will  be  assumed  equal  to 

.5.  Similar  justification  will  be  used  in  computing  n  .  A  knowledge 
of  an  actual  stability  limit  will  be  used  to  compute  n  ;  if  the 
numbers  are  reasonable,  i.e.,  n  >0  ,  it  will  be  considered  not  only 
fortunate  but  a  rather  important  confirmation  of  the  application  of 
the  theory.  It  should  be  noted  that  one  confirmation  is  already 
given  by  the  wake  theory  in  that  the  imaginary  counterpart  of  Equation 
(1)  is  negative,  indicating  a  lag  of  the  burning  rate  behind  pressure. 
This  was  confirmed  by  the  time  lag  theory  and  associated  experimen¬ 
tation  at  the  short  length  limit. 

The  final  quantity  for  consideration  is  the  fuel  injection 
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velocity  variation.  If  chamber  pressure  is  held  constant  and  mixture 
ratio  varied 

Vp  04  (0/F)"8  (9) 

where  s  A*  0.45  in  the  range  of  interest  and  is  actually  a  weak 
function  of  0/F  Itself. 

Now  combining  Equations  (5b)  and  (6) -(9),  there  results 
the  computational  equation  with  (0/F)  ■  1.0  and  (q/q  )  •  (T  /T) 


t  l(^> 


rv  0  )  10-*) 


(10) 


which  defines  .  The  enclosed  table  represents  computations  made 
for  the  old  3"  cylindrical  engine  fuel-on-oxldlzer  impinging  Injector 

if  if 

with  r  •  r  ■  0.035  in.  The  propellants  are  ethanol  and  liquid 

J  Jo 

oxygen  at  a  chamber  pressure  of  300  psia.  The  reference  point  is 
0/F  -  1.0,  Lq  -  8.0  Inches.  The  Injector  design  point  is  0/F  -  1.4. 


p  -  300 

psia  Yp 

c 

-  .9 

T  -  .2 
c 

q  -  10 
o 

0/F 

T 

chamber 

(°F) 

^chamber 

(ft/sec) 

T/T  - 
o 

Vq 

c/c 

0 

Y 

°/J 

Tf 

L/L 

o 

1.0 

4300 

3770 

1 

1 

0.10 

1.817 

1 

1.4 

5560 

3940 

1.292 

1.044 

0.14 

1.827 

.798 

1.65 

5840 

3725 

1.359 

1.019 

0.18 

1.968 

.833 

2.2 

5820 

3875 

1.352 

.988 

0.22 

2.150 

.937 

b  -  1 

■ 

to 

oo 

• 

> 

4-* 

*"t 

■ 

O 

fc 

-  2.: 

1 

211 

/ 

.rt/„..45n  -  .8  / 

(0/F)  -  V 

<l/lo)  (c 

o/c)(To/T) 

t 

Ill 


0/F 

t 

0/F,45n-*8 

n 

1.0 

1 

1 

1.4 

1.005 

.765 

.011 

1.65 

1.041 

.739 

.436 

1.8 

1.066 

.728 

.580 

2.2 

1.152 

.727 

.880 

It  is  noted  that  in  all  cases  n  >  0  which,  at  least  is 
reasonable  even  if  accurate  numbers  are  not  known.  It  should  be  borne 
in  mind  that  at  0/F  •  1.0  the  stability  map  has  almost  closed;  in 
fact,  all  stability  maps  close  rapidly  at  low  0/F.  It  is  felt  that 
such  application  is  stretching  the  theory  considerably  since  it  would 
then  be  applied  near  both  the  upper  and  lower  length  limit.  For  this 
application  it  is  preferable  to  only  consider  the  results  for  0/F>1.4. 
At  low  0/F  the  motor  performance  also  deteriorates  rapidly  making  the 
assumption  b  -  constant  invalid.  In  fact,  b  would  Increase,  tending 
to  raise  the  computed  n  to  more  realistic  values.  A  salient 
characteristic  of  the  low  length  limit  is  that  it  is  concave  toward 
long  lengths  on  a  plot  of  0/F  vs  length.  To  match  this  characteristic 
is  a  major  achievement  here.  It  can  be  done  with  a  variable  n  , 
although  for  the  injector  used,  there  is  no  information  to  tell  whether 
or  not  this  is  realistic.  It  can  also  be  done  with  a  fixed  n  .  For, 

say,  n  *  1/2  L/L  at  first  decreases  due  to  0/F  to  a  negative  power 

°  2 
in  Equation  (10).  At  sufficiently  large  0/F  ,  however,  becomes 

overpowering  and  tends  to  increase  L/ Lq.  This  accounts  for  another 
characteristic  of  stability  maps:  at  high  0/F  only  a  gentle  slope 
to  the  right  occurs.  This  is  due  to  the  compensating  effects  of  n 
and  The  fact  that  different  injection  types  should  produce  dif¬ 

ferent  types  of  instability  maps  clearly  appears  through  the  form  of 
variation  of  a*.  However,  exactly  what  these  dependences  are  is  not 
yet  clearly  defined.  For  Instance,  in  the  fuel-on-oxidizer  injection 
system,  the  resultant  momentum  axis  of  the  Impinging  jets  changes 
with  0/F.  This  is  not  true  f<?r  thq  like-on-like  system. 
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Concerning  changes  in  r  ,  a*  probably  Increases  with 
Tj  for  n  <  1.  This  in  turn  increases  the  lower  length  limit,  a 
fact  observed  experimentally  at  Princeton.  It  produces  an  effective 
increase  of  the  sensitive  time  lag.  Such  a  fact  is  also  borne  out 
on  the  larger  engines  built  today;  there  is  no  longitudinal  instability 
at  quite  long  lengths.  Injectors  of  the  same  orifice  size  but  of 
different  type  also  follow  the  predicted  trend.  For  a  showerhead 
which  produces  a  Larger  drop  size  than  impinging  jets  of  the  same 
size,  instability  does  not  occur  until  quite  long  lengths.  Concerning 
other  propellant  types,  RP-1  while  having  a  larger  heat  of  combustion 
than  ethanol  also  produces  a  larger  T*^  such  that  q  is  roughly 
independent  of  propellant  type.  Then  while  J/l  increases  c  also 
increases  and  there  is  a  compensating  effect  in  Equation  (6).  The 
question  of  propellant  change  once  again  revolves  about  the  drop  size. 
The  factors  involved  include  surface  tension  and  liquid  density  and 
viscosity  as  well  as  the  gas  phase  properties  and  injector  type.  It 
may  well  be  possible  to  determine  the  drop  size  change  from  instability 
data.  Therefore,  using  Reference  32  ,  it  is  concluded  from  the  lower 

length  limit  locations  that  a  switch  from  ethyl  alcohol  to  iso-octane 
as  a  fuel  produced  a  smaller  drop  size.  This  would  also  be  predicted 
from  the  injectors  used  since  r^  was  decreased  in  the  iso-octane 
injector  over  that  used  for  alcohol  (.031  vs.  .035  inch  diameter), 
lending  further  support  to  the  theory.  Unpublished  RP-1  data  suggests 
no  change,  all  other  factors  being  the  same. 


Comparison  with  n,  Theory 

Including  the  phase  effect  Equations  (1)  and  (2)  along  with 


our  approximations  for 


R^  may  be  written  from  the  wake  theory  as 


w,  <P  -  -  A  +  qfr” 

+  -y- A  *  “■ 


3  -4WO  Va  2s 


l 


(11) 


I 
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From  the  n-  2*  theory 

}ftr  / ,  -ccj  *  2*  \ 

-  s  *  Cl  -  c  c  )  (12) 

mff<P 

The  equivalence  of  these  two  relations  can  only  be  made  at  a  stability 
limit,  where  periodic  solutions  actually  exist.  Therefore,  variations 
in  n  and  £**  with  actual  droplet  parameters  are  under  the  con¬ 
straint  that  one  moves  on  a  stability  line.  Separating  into  real  and 
imaginary  parts,  Equations  (11)  and  (12)  become 


>1  * 


(13) 


(J*  S'*  -  -  V3A  <“> 

Abandoning  the  time  lag  theory  postulate  that  n  and  £**  are 
frequency  independent,  we  now  say  that  at  the  short  length  limit  they 
may  vary  with  U)  while  still  retaining  the  form  of  Equation  (12). 
Given  the  droplet  parameters  which  appear  in  A  and  B,  Equations  (13) 
and  (14)  may  be  used  to  determine  n  and  2** as  functions  of 
frequency,  CO  .  It  remains  to  be  ascertained  if  the  relationships 
between  n  ,  t*  ,  and  c *9  can  be  similar  to  those  obtained  from 
the  time-lag  theory,  and  if  so,  in  what  range  of  droplet  parameters 
does  the  similarity  occur. 
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NOMENCLATURE 


a 

c 

c 

P 

i 

J 

L 

m 

M 

7K 

n 

O/F 

P 


droplet  diameter 
speed  of  sound 

specific  heat  at  constant  pressure 
imaginary  number  <fT) 

stoichiometric  mass  ratio  of  fuel  to  oxidizer 
chamber  length  at  high-frequency  stability  limit 
total  mass  flux  in  radial  direction 
Mach  number 

space  dependent  part  of  m  perturbation 
interaction  index  or  exponent  in  Equation  (8) 
mixture  ratio 

pressure  or  exponent  in  droplet  diameter  and  density 
relation 

stoichiometric  heat  of  reaction 
radial  coordinate 

perturbation  magnitude  of  flame  movement 

Reynolds  number 

exponent  in  Equation  (9) 

exponent  in  Equation  (7) 

temperature 

mass  fraction 

defined  by  Equation  (10) 

ratio  of  specific  heats 

viscosity 

axial  variable 

density 

frequency 


Subscripts : 

c  core  centerline  value  (at  £  ■  0  unless  indicated  by 

£  in  argument) 
f  quantity  at  flame 


I 
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fuel  or  fuel  side  of  flame 

oxidizer,  oxidizer  side  of  flame,  or  initial  value 
free  stream  value 

Superscripts: 

-  steady-state  quantity 

*  dimensional  quantity 
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VI.  HIGH- SUBSONIC  MACH  NUMBER  EXTENSION  TO  THE  LINEARIZED  THEORY 

Using  a  simplified  model  for  the  rocket  engine,  whereby  all 

★ 

combustion  was  concentrated  in  a  narrow  zone,  Crocco  and  Cheng 
studied  the  stability  behavior  of  the  flow  system  in  the  rocket  chamber 
when  small  perturbations  on  the  steady-state  operation  were  intro¬ 
duced.  The  analysis  assumed  that  the  flow  Mach  numbers  in  the  chamber 
were  so  low  that  steady-state  thermodynamic  properties  were  constant 
on  either  side  of  the  flame  front.  In  effect  this  involved  ignoring 
terms  which  were  of  the  magnitude  of  the  Mach  number  squared. 

Since  that  time,  rocket  engine  configurations  have  been 
built  with  increasingly  lower  contraction  ratios  (M  ">  .5  is  not 
uncommon),  which  makes  the  terms  of  the  order  of  the  Mach  number 
squared  more  important. 

In  the  past  year  work  has  been  done  to  extend  the  analysis 
of  Crocco  and  Cheng  to  include  these  effects  of  higher  chamber  Mach 
numbers.  Most  of  the  assumptions  and  the  methods  of  analysis  of  the 
original  work  are  retained.  In  general,  a  steady-state  flow  system 
using  the  concentrated  combustion  front  is  used,  small  perturbations 
on  the  steady-state  are  applied,  and  the  stability  behavior  is 
examined.  The  chamber  oscillations  are  supported  by  means  of 
variations  in  the  burning  rate  which  are  described  by  means  of  the 
time-lag  concept. 

A  combustion  chamber  is  separated  into  two  regions;  Region 

0,  which  contains  stagnant  gaseous  combustion  products  with  liquid 

* 

propellant  droplet  flowing  through  it  with  constant  velocity  Ua  , 

^  o 

and  Region  1  which  contains  gaseous  combustion  products  of  the  same 

k 

composition  flowing  with  some  uniform  axial  velocity  U^  .  The  flame 
front  separates  the  two  regions.  The  propellant  droplets  are  con¬ 
sidered  to  be  instantaneously  and  completely  burned  at  the  flame 
front.  One-dimensional  flow  is  assumed.  The  following  equations 

describing  the  steady-state  conditions  in  the  chamber  are  written: 

★They  studied  the  stability  behavior  of  the  flow  system  for  the  more 
general  case  of  distributed  combustion,  as  well. 
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jf  *  r 

continuity 

(1) 

P.#*  fj  %'*-  !!*+?<' 

momentum 

(2) 

r  -f+k  u/‘ 

0 

energy 

(3) 

starred  quantities  are  dimensional. 

The  quantities  are 

non- 

dlmenslonallzed  with  respect  to  conditions  in  Region  0  in  the 
following  manner: 


The  equations  are  rewritten: 


ft  2i  *  ? 

h  *  5  “  7- +  j[  S-1 

T  ■ ¥  V 


The  perfect  gas  assumption  yields: 


Combining 

P-  '  + 


(4)  and  (5)  : 


(4) 

(5) 

(6) 

(7) 


(8) 
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or 


l- 


where 


i  LEll  ^"<) 

KWj 


*s  -y- 

u.. 


(9) 


Non-steady  state  conservation  equations  are  written  for  either 

Region  1  or  Region  0  as: 

continuity; 


Jt 

momentum: 


9  k 


(10) 


f&-  + 


Pu.  ilt-  +4 _ i 

J  )x  *  7  ) 


J£- 

**X 


(11) 


The  energy  equation  is  replaced  by  the  isentroplc  assumption 


as  follows: 


P*J*r 


(12) 


Let 


P9  P  *  P  f  *  f  *  ^  '  u’m  fL  +-  U 

where  p',  J3  '  ,  u'  are  perturbation  quantities  so  small  that 
their  products  and  squares  can  be  ignored. 

If  the  steady-state  equations  are  subtracted  from  (10) 
and  (11),  they  become: 


,  / 

”  **  ^  Ta*  ^  =  ^  continuity  (13) 

F f-  P  2  ■  m  -i—  P  momentum  (14) 

9k  m  t  9t 

These  equations  correspond  to  the  wave  equation  for  a  one-dimensional, 
uniform  flow  field.  The  general  solution  can  be  written: 


*4 


•  V  ll 


-Ilf- 


u  -  f,  [x-d.Dt 

f-  fe  «-f) 

p'-  '  -f  ({’--f,') 


t  f 

Note  that  the  functions  f^,  and  f2  are  not  the  same  in  Regions 
0  and  1. 

«  » 

Since  f.  and  f_  are  arbitrary  functions  of  their 
l  *  »  •  » 

arguments,  and  since  periodic  solutions  for  u  ,  V  ,  p  are 

i  i 


desired,  the  functions 

i  f^  and  f2  are  written: 

f  * c,  •«/» 

S  (t  -  ar  x ) 

(15) 

K  *c»  txP 

S(t-C^A) 

(16) 

where 

a--  “  ub 

*  a- c 

s  »  c 

n. 

Using  (15)  and  (16),  the  perturbations  become 


\JL  m  TAX)  St  (17) 

S'  *  £  ^X)  **P  St  (18) 

r? 

p  ""jr  <Tcx)  %xp  st  (i9) 

where  •*!«)  ■  Cr  «*p  (r^S  *)*  ^  9  **  (20) 

0“oo  =  Cr  exp(ar9*)-^  %xpC-o^Sx)  (21) 


(20)  and  (21)  represent  four  equations,  two  for  Region  0  and  two  for 
Region  1,  so  that  there  are  four  constants  to  be  determined. 
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Two  boundary  conditions  and  two  matching  conditions  are  available 
to  determine  the  four  constants. 

At  the  injector  face  ft-'  ■  0  and  x  ■  0  and,  therefore, 
Equations  (17)  and  (20)  yield 


C 


-  1 


(22) 


At  the  chamber  exit  plane  x  ■  1  and 


where  OC  is  the  nozzle  specific  admittance  coefficient  (Reference  4.) 
Then,  Equations  (17),  (18),  (20),  and  (21)  yield 


I* 

. _ 

I  - 

L  c. 

S  (CLr  -  43, ) 


where 


e*p  S  (a,.  -  cy  ) 


(23) 


Using  (22)  and  (23),  the  space-dependent  parts  of  the 
perturbation  become 


1^0 1)» 

<V. 

•xp  (•  S  O^x  )  £|  -  exp  [*  S  (flrp- 

V‘]] 

(24) 

*;<*>* 

c*. 

•XpC-SQ^x)  £l«-  exp[-s(»r#" 

(25) 

vjoo* 

c>, 

exp C~30^  *)£/-«  exp  [s  (a*. 

-<y)0-0]J 

(26) 

« 

•  ip  (-  3a  Xl)  ♦  6  exp  [s  ter; 

(27) 

At  the  flame  front,  the  steady-state  conservation 

equations  for  mass  and  momentum  are: 
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*i  %  - 1  u' 

+  ft  ai’  •  f,  v+  -f 


mass 


momentum 


The  assumption  has  been  made  earlier  that  J5  *i.  is  a 
constant,  independent  of  time.  Ibis  implies  that  the  same  amount 
of  mass  is  always  entering  the  flame  front  due  to  the  liquid  droplets. 
The  amount  of  mass  leaving  the  flame  front  is  dependent  upon  both 
the  amount  of  gas  and  liquid  entering  the  flame  front  and  the  amount 
of  gas  produced  in  the  flame  front.  Under  steady-state  conditions, 
no  gas  enters  the  flame  front,  and  gas  is  produced  in  the  flame 
front  at  the  same  rate  liquid  mass  enters.  Under  a  perturbation  of 
the  steady-state  pressure,  the  burning  rate  also  undergoes  a 
perturbation.  Under  such  a  perturbation  the  change  in  the  amount 
of  mass  produced  in  the  flame  front  is  equal  to  the  difference  between 
the  incoming  liquid  mass  flow  and  the  outgoing  mass  flow  less  the 
incoming  mass  flow.  That  is 


♦  J>,  -  fo  <28) 

Similarly,  the  amount  of  momentum  produced  under  steady- 
state  conditions  must  be  equal  to  the  incoming  liquid  momentum  plus 
the  pressure  difference  between  the  pressure  in  Region  0  and  the 
pressure  in  Region  1.  Under  a  perturbation,  the  amount  of  incoming 
liquid  momentum  i*-  constant.  The  change  in  the  amount  of  momentum 
which  must  be  produced  is  thus  equal  to  the  difference  in  pressures 
between  Regions  0  and  1  under  the  perturbation,  plus  the  outgoing 
gas  momentum  less  the  incoming  gas  momentum  less  the  incoming 
liquid  manen tun.  This  relationship  may  be  written: 

A  "V  =  ^  Vj‘  -f.  U.'*>  f,  (2’) 

The  assumption  is  now  made  that  ^mv,  the  amount  of  momentim  produced 
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in  the  flame  front  is  equal  to  the  product  of  Am,  the  change  in 
mass  produced  in  the  flame  front  under  perturbation  and  Ujj  .  That 
is,  it  1 8  assumed  that  all  mass  produced  in  the  flame  front°is 
released  with  constant  velocity  equal  to  the  velocity  of  the  incoming 
liquid  droplets.  This  relationship  is  written: 


Amr  »  (A-m)  .  oo> 

Using  (28)  as  substituting  (30)  into  (29),  (29)  may  be  rewritten 

'  2 

(neglecting  the  higher  ^rder  term  u  ): 


^“,VU<UV 


(31) 


Rewriting  the  quantities  used  above  as  the  sum  of  their  steady-state 
valines  and  a  small  perturbation  of  this  steady-state  value,  the 
following  equation  results  when  the  steady-state  momentum  equation 
is  subtracted: 


(32) 


This  equation  is  solved  to  yield  the  ratio  between  the  perturbation 


quantities  pQ  •  and  p^ 


?  ' 
Jj _ 

P‘ 


I  +  $  u, 


-  UJ 


(33) 


V 


Mass  conservation  at  the  flame  front  requires  that  A  m  *  m^  if 
m^  is  change  in  burning  rate  under  a  small  perturbation.  It 
follows  that 


Am-?,* 

*,  *• 


(34) 


where  A  b  is  the  fractional  increase  in  burning  rate. 

The  assumption  is  now  made  that  all  propellant  elements 
have  the  same  sensitive  time-lag.  In  this  case,  ^ib  may  be 
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expressed  as  (see  Crocco  and  Cheng,  Ch.  I) 

where  is  the  sensitive  time-lag. 

In  general: 


7-7-- 


n 


ret 


f 


p  it')  Jt' 


(35) 


f(t)  represents  the  processes  which  oust  accumulate  before 
combustion  can  occur.  That  is,  combustion  occurs  at  time  t 
according  to  the  criterion 

Cr  W>*'m  CA 

where  has  a  certain  value. 

n  is  an  interaction  index  dependent  upon  propellant 
properties  and  defined  as  71 ■  where 

Z  is  an  arbitrary  variable. 

Since  all  preparatory  processes  take  place  in  Region  0, 
f(t)  is  a  constant  and  p(t')  ■ 

Thus,  (35)  may  be  written 

r-r«-n.  f*  P'ct)H‘ 

Jtf 


JT 

it 


f  pxtt- 

/tk  *  [V ■  pt*  ■ 


The  time  used  in  the  last  expression  above  is  a  Lagrangian  Coordinate 
time;  in  order  to  be  consistent  this  is  rewritten 


P'iyu-  ?(*,*-*) 


(36) 


t  is  now  time  in  the  coordinate  system  stationary  with  respect  to 
the  combustion  chamber.  £  is  the  position  of  the  flame  front, 
is  the  axial  position  at  which  the  propellant  element  that  burns 


at  time  t  first  becomes  sensitive. 

Tha  assumption  is  now  made  that  la  ao  cloaa  to  ( 

•  • 

tha  tlaa  lata  undtr  conaldaratlon  that  p  (  fc  »t)  p  (K  »t). 
Thla  Implies  that  tha  production  of  mass  aay  vary  with  tlaa  but 
location  of  tha  flaa»  front  doaa  not. 

Undar  thla  assumption  (36)  aay  ba  written 


^  [exp  (st)- exp  met-W] 
•  ■nr  ^  «xp  $t  p-  (r &&)] 


This  axprasslon  may  ba  combined  with  (34)  to  ylald 


(37) 


(38) 


wr  OJ  exp  9-t  [i-*»p(-» '?’)]-  X (39) 

f  V 

J,  » 

This  may  ba  vrlttan  In  the  form 


it.  }  n  1  [i- *3"  r{'  ~  o 

’»  Ja 

From  (17)  and  (18) ,  It  la  found  that 

If  these  two  relationships  and  (33)  are  substituted  in  (40)  the 
following  equation  results: 


(40) 


H  f  n  r  D-  **p  ** 

-  u,  [i»  a  V]*?,  &*8< 


bw&ut)  (41) 


where 


D0*e««pV)-(»-it)u1c;  TD-8«pY] 

VmsXO-x),  DmT,-  3u5-»  %•  ■=£}%* 


X  represents  the  axial  position  of  the  flame  front.  This  equation 
may  be  reduced  to  the  form: 


P+M ««p  (-i*x)*[N+L*ipC-m)]  «»p  V 
*  [etE «xp(-«»x)JexpC- if')* [f*r**t t-f»i>J«xp(f-**)-0 


(42) 


where  E,  F,  1,  M,  N,  end  P  are  defined  later. 

In  order  to  find  the  neutral  stability  relationships, .A. 
is  set  ■  0,  S  ■  la!  ,  and  the  corresponding  ^  f  • 

If  (42)  is  separated  into  its  real  and  imaginary  parts, 
the  following  two  equations  are  obtained: 


CL-t-C  uUA-#-tjm.<A-P>0  (*3> 

b*J  +  c 

where 

si  (J^l  (f)  *  CJ  S 

a  ■  P+ M  ee*  C  «>x  j  b  «  -  M  -awe  t  o)X 

c  ■  W(6rc  cm  a«t  *  16M0  ♦*  (£«b,  **£•*>*«**)♦  ^•**  ft**  >A*f) 

d  ■  *)(0u,  CM  *w>X  -  *«*)♦  L.  [fo»  CftJ|*-cX)-  C*mX*M )] 

eiE»C  e«x  j  f  a  -  C  mwv  aux 

g  BP  Cs^ (tec  SMX U4.  tt*»x  ♦  -ax) ♦  ( 4wi  i&A  Uu>x  +**>*)] 

h  *  F[^  (uciui  ♦  tea  lw«  a-cx)-  (W  (*i*it*t  *4*  <****  ++*))] 

M SCO- B.  e. T C«- *))( cr,»r-iHv?Xi-f  5,) 

N«(o*u,r,  x<i-*>Xs.J  f  u.) 

i*(D»E.S,XU-»))(5,  jnrO-(fc,*2.Xi-  fs,) 

*■-«-.  *r(D- E, I.  Tit-*)) 

ps-t.j  >iro>»tt,e,T;  u-*0 

B  and  B.  are  the  real  and  imaginary  parts  of  B  • 
re  Tf" 

»r.- 

eCre  and  are  the  real  and  imaginary  parts  of  the 

noxsle  specific  admittance  ratio. 

After  some  manipulations,.  (43)  and  (44)  are  solved  to  give 


tha  following  relationships  for  a 


CdtM 


XF*  +C* 


H*  +  J* 


(45) 


(46) 


CL*  +  + 

50*  L  COlI  i-bJ  -*f-  H) 

Q  -  t  (bo-ad  +  «.)(  -£y) 

H  5  i  (afcf  If  *lfc  -  d  ^0 

j  5  2  (af -ie 

The  right-hand  side*  of  Equations  (45)  and  (46)  above  are 
complicated  function*  of  CO  and  n  .  Since  no  simple  solution  has 
been  found  for  the  general  case,  an  iteration  procedure  has  been 
adopted.  In  this  procedure  for  given  values  of  H^,  k  ,  and  Y  , 
an  6>  is  selected.  A  trial  value  of  n  is  then  chosen  and  the  right 
and  left-hand  sides  of  Equation  (45)  are  calculated  and  compared  to 
each  other.  Succeeding  values  of  n  are  chosen  until  the  left-hand 
side  is  equal  to  the  right  hand  side.  This  will  give  a  value  of  n 
which  corresponds  to  the  selected  value  of  CO  .  Using  these  values 
of  a)  and  n  0  -  uS  may  be  calculated.  Thus,  a  value  of  CO  S 
for  the  given  n  is  also  obtained.  Using  this  technique,  the  neutral 
stability  curves  for  u)  vs  n  and  vs  n  ,  may  be  plotted.  Since 
the  neutral  stability  curve  for  S  vs  n  is  also  obtained. 

In  the  simple  case  of  combustion  concentrated  at  the 
injector  face  and  short  nozzle,  X  ■  0,  ■  0,  and  the  equations 

are  greatly  simplified: 


Cos  s  -  - 


Cos  0 


(47) 


(48) 


where  A  ,  I  ,  I  ,  H  ere  no  longer  functions  of 
Iteretlon  procedure  Is  not  necessary. 

Results  end  Conclusions 

At  the  present  time,  numerical  results  ere  available  only 
for  the  simple  case  of  short  nozzle  and  combustion  concentrated  at 
Injector.  Numerical  analysis  Is  procedlng  for  the  more  Involved 
cases  of  arbitrary  axial  location  and  long  nossles* 

Equations  (47)  and  (48)  were  solved  using  Chamber  exit 
■  .417  end  Mg  ■  .207.  The  solutions  ere 
shown  In  Figs*  VI-  1,  *2,  *3.  Some  effects  of  increasing  Mach  number 
are  readily  seen.  First,  there  Is  a  shifting  of  the  minimum  n 
which  will  support  unstable  oscillations  toward  larger  values  of 
n  .  This  can  be  seen  by  comparing  the  minimum  n  for  Mg  «  .417 
and  Mg  ■  .207  in  Figure -1  or -2,  The  linear  (in  Mach  number) 
analysis  of  Crocco  and  Cheng  predicted  a  minimum  n  of  .46  regard¬ 
less  of  Mach  number.  As  can  be  seen  from  Figure  4  or -2,  agreement 

■  .207.  However,  for  Mg  •  .417 
a  noticeable  increase  occurs  (  n^  -  .54).  This  seems  to  indicate 
that  a  larger  forcing  function  Is  required  at  the  fundamental 
frequencies  of  the  chamber  in  order  to  produce  instability  as  the 
Mach  number  in  the  chamber  is  Increased.  The  effect  is  not  too 


with  this  value  is  close  for  Mg 


Mach  numbers  of  Mg 


large,  for  a  Mg  ■  .620  ,  nmin  Increases  only  to  .56. 

A  second  effect,  due  to  consideration  of  all  orders  of  Mach 
number,  may  also  be  seen  in  Figures  -1,  -2 ,  and  -3.  This  effect  is 
the  existence  of  a  maximum  n  above  which  no  stable  region  can  be 
found.  At  this  maximum  n,  w  ,  the  dimensionless  frequency,  is- 
exactly  half-way  between  the  fundamental  frequencies  of  the  chamber. 
The  fundamental  frequencies  of  the  chamber  are  the  frequencies  which 
would  be  present  in  the  chamber  if,  at  both  ends,  it  reflected 
compression  waves,  as  compression  waves,  and  expansion  waves  as 
expansion  waves  with  no  phase  or  delay  in  reflection.  In  Figure -1, 
these  fundamental  frequencies  are  found  as  the  frequencies 


The  frequencies  of  oscilletion  at  n  max  ara  furthaat  from  tha 


fundamental  vaiues  and  thus  ara  hardaat  to  support.  At  however 

thaaa  oscillations  can  ba  supportad.  This  implies  that  for  n 

graatar  than  n  ,  oscillations  at  any  fraquancy  can  ba  supportad 

and  tha  chamber  is  absolutely  unstable. 

Chamber  Mach  number  has  a  profound  affect  upon  the  value 

of  n  .As  can  be  seen  from  Figure  -1,  n  max  for  M-  -  .417  is 
max  £ 

equal  to  1*85.  For  M_  -  .207,  however,  n _ is  so  large  as  to  be 

h  max 

off  the  scale.  This  is  in  agreement  with  the  results  of  Crocco  for 


the  case  of  low  Mach  number.  He  predicted  no  n  ,  and  it  is  found, 

max 

when  all  orders  of  Mach  number  are  considered,  that  for  low  Mg, 

n  of  practical 


n  has  a  value  which  is  above  the  range  of 
max 


interest,  and,  in  the  limit  of  very  small  Mach  numbers,  n  tends 

max 

to  infinity. 

Another  consequence  of  considering  all  orders  of  Mach 
number  is  the  shifting  down  of  the  0>  vs.  n  curve  and  the  shifting 
up  of  the  f  vs  n  curve.  This  effect  can  be  seen  in  Figures  -1 
and  -3  if  one  looks  at  the  curves  for  the  same  mode  numbers  at  the 
two  different  Mg's,  .417  and  .207.  It  should  be  noted  here  that 
the  m  numbers  come  from  Figure  2.  From  the  solution  of  (47),  it 
is  found  that  has  the  closed  curve  shape  seen  in  Figure  2  and 

that  the  horizontal  line  of  symmetry  may  be  located  at 
where  m«0,  1,2,3,  ...  In  Figure  -3  results  using  m  «  0 
and  m  «  1  are  shown.  Since  an  m  and  an  /,  must 

be  specified  for  a  given  S  curve.  In  Figure  -3,  the  m  and 
on  the  various  curves  correspond  to  the  m  and  JL  of  the  corres¬ 
ponding  a >  and  curves  which  were  cross-plotted  to  yield  the 

5  vs  n  curve. 

Ihe  shifting  of  the  (J  vs  n  curve  to  lower  cJ  for  higher 
Mach  numbers  at  the  fundamental  frequencies  can  be  explained  as  an 
effect  due  to  the  change  of  wave  propagation  time  in  the  chamber. 
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The  time  required  for  a  vcvi  to  travel  from  on*  ond  of  tho  chamber 
to  tho  othor  ond  bock  again  in  steady-state  ia: 

L*/  ia  tho  factor  vo  uaod  oar  liar  to  nondim(ut.iic 
c0 

time.  Thus  tho  equation  above  becomes  in  non-dimensional  fora: 


Thia  oxprosaion  ia  precisely  equal  to  A.  Since ,  for  fundaeental 
frequenciea  CJ  *  HF  ,  then  uaing  (49): 


Thia  implies  that  ea  the  wave  travel  time  increaaes,  the  fundamental 

frequency  decreaaes  for  given  JL  .  Since  T  Increases  with  Mach 

number,  CJ  oust  decrease  with  Mach  number. 

The  values  of  at  nnln  are  (***•)  TT” 

therefore,  does  not  shift  at  n  .  for  Increasing  Mach  number 

min 

The  quantity  o  is  the  product  of  ^  and  **»  o  and  therefore  it 
shift  up  at  n^n  for  Increasing  Mach  number. 


Some  general  comments  should  be  made  about  the  curves  in 

Figures  -1,  -2,  and  -3.  First,  all  the  curves  plotted  are  continuous 

and  have  infinite  slopes  at  n  _  and  n  .  .  This  can  be  Shown 

r  max  min 


analytically  for  the  simple  case  of  short  nozzle  and  concentrated 

combustion  at  the  injector.  Secondly,  if  the  Mach  number  is  allowed 

to  become  very  small,  the  curves  exactly  duplicate  the  results  of 

the  linear  analysis  for  n  in  the  range  of  practical  Interest  (up 

to  about  S).  At  very  large  n  ,  however,  there  will  exist  an  n 

max 

even  for  very  small  Mach  numbers  and  the  curves  will  still  be 
continuous.  The  curve  of  Figure  -3  will  become  several  curves,  all 
continuous,  if  all  possible  combinations  of  m  and  JL  are  considered. 
These  curves  will  cover  practically  the  whole  area  between  nmin  and 
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with  unstable  rations,  The  curves  arc,  in  general,  slightly 
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mar*  spread  out  for  given  si  and  JL  for  larger  Mach  numbers,  thus 
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predicting  an  increase  in  the  alee  of  the  Instability  regions. 


•In  conclusion,  for  the  simple  case  of  short  nossle  and 


concentrated  coabuetlon,  two  effects  on  stability  are  to  be  recognised 
as  higher  Mach  nunbers  are  considered.  The  first  is  the  shifting 


■> 1  1 


of  n^  to  higher  values  of  n  indicating  Increasing  stability 
The  second  is  the  spreading  out  of  the  S  vs  n  curve  away  fron 
n  .  which  indicates  decreasing  stability* 


IHC.fijaQiLL 


e  .  L* 
6  Co* 


local  speed  of  sound 

speed  of  sound  in  stagnant  burned  gases 

chamber  length  from  injector  to  entrance  of  nozzle 

time  required  for  wave  to  travel  chamber  length  in  stagnant 


axial  length  from  injector  end  of  chamber 
mean  gas  velocity  in  axial  direction 


Mach  number  at  entrance  to  nozzle 


pressure  of  burnt  gas 

density  of  burnt  gas 

temperature  of  burnt  gas 

enthalpy  of  burnt  gas 

density  of  liquid  propellant  element 

velocity  of  liquid  propellant  element 

constant  relating  tT  and 


amplification  coefficient 


angular  frequency 
neutral  angular  frequency 


interaction  index 


time  lag 


critical  value  of  time  lag  corresponding  to  neutral 
oscillations 
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specific  admlttanca  ratio  of  daLaval  nossla 
V*  axial  location  of  tha  concantratad  combustion  front 

superscript*  indicataa  dimensional  quantity 


*  i 


4 


•uparacript'  indicataa  a  aaall  parturbation 

aubacrlpto  indicataa  quantity  avaluated  in  Ration 

subscript^  indicataa  quantity  avaluatad  in  Ragion  1 
subscript^  indicataa  quantity  partaining  to  liquid  propallant 

_  ovar  a  quantity  indicataa  steady-state 

subscripts.  and  _  indicate  tha  raal  and  imaginary  parts  of  a 
•  m  quantity,  raapactlvaly. 
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momentum: 
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Momentum: 

%. j(\  {*’ 4-'  if  jjf  !• -  *,  [  1 1 
*+(%  •  tf&jfj-  «* 

0  -  Momentum: 

Ji,,  fa  ff  r‘j] 

-f  £  iut  (•■**"•  ="•**'")  *,  S /T’']‘,/ 

Entropy  Equation; 
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Equation  of  State: 

E  f  < r  fT'-tb  ?‘',lhTjS.tW 

A  /  Tt  tij*  nv)S  /  Aicwt.MW  M>i 

+  ^\(iR  ^  p  4(t  *-6  H><V^ 

-  H  +(e“<^  •-“*%  (S^/F) 


II.  SOLUTION  FOR  STANDING  WAVE 

The  symbols  A,  IB ,  C,  1)  and  £  appearing  in  the  equations 
below  have  the  same  definitions  as  the  corresponding  symbols  appearing 
in  the  expressions  of  part  (I)  page  132  (corresponding  to  spinning 
wave  solution). 


Continuity  Equation: 


0  Mo»entu»: 


- liwfr 
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5  [»,,  [*.*<•“-  .**)]? 

-  Momentum ; 


£  {%{i +*"***)}& [*(*** W*9} 
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Entropy  Equation: 


Equation  of  State  ; 

£[«,[*•* *-']  J.  <5,^* 

*  G„,  [i  («".  l  <S„  /wj 

From  the  expansion  forms  of  the  inhomogeneous  parts  of  the 
second-order  equations  (for  standing  as  well  as  tangentially  spinning 
waves)  it  can  be  seen  that  m  and  n,  defined  in  Equation  (5.8)  take 
on  the  following  values:  0,  +  2. 

The  expansion  forms  of  the  functions,  depending  on  the  in¬ 
dependent  variable  ,  appearing  in  the  inhomogeneous  parts  of  the 
second-order  equations  will  follow: 


T/ (5* /?")-£  ^  Z* 

^  *,  (s^  m 

u.  j;^/f)-£c 

vfo  r» (s*  mf  ■  f  Jw  ^  /v) 

V  I,  (S^sr)^  q;  fsu<^)-  y  ^  J5'<V*> 

For  the  second  order  solution  n  tskes  on  the  vslues  0  end  +  2.  Also 
note  that  the  relation  (S2J^  ^  "  J -2})  ^2jJ  ^  (when  2 

Is  an  Integer  )  was  used  in  expanding  the  functions  appearing  in  the 
inhomogeneous  parts  of  the  second-order  equations. 
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APPENDIX  B 

CONCEPT  IN  ROCKET  BAFFLE  DESIGN 

Introduction 

It  is  well  known  thst  the  occurrence  of  transverse  spinning 
mode  high  frequency  combustion  instability  in  the  liquid  propellant 
rocket  engine  can  often  be  eliminated  by  the  use  of  baffles.  These 
baffles  are  solid  surfaces  which  protrude  from  the  injector  plate, 
are  perpendicular  to  it,  and  extend  some  determined  distance  into 
the  rocket  combustion  chamber.  The  pattern  or  arrangement  of  these 
baffles  required  to  prevent  combustion  Instability  cannot  as  yet  be 
predicted  a  priori,  and  generally  it  may  be  said  that  such  design  is 
still  an  art  rather  than  a  science. 

In  some  situations,  however,  there  are  certain  rational 

guidelines  to  this  design.  It  appears  for  some  configurations  that 

a  design  which  segments  the  chamber  into  smaller  sections  will  be 

successful  in  preventing  the  spinning  mode  in  the  chamber  although 

16 

the  standing  mode  may  still  be  present  .  This  is  true  even  if  the 
protrusion  distance  is  rather  small  compared  to  the  chamber  length 
provided  that  the  baffle  protrudes  part-way  into  a  region  of  active 
combustion.  If  this  is  true,  then  segmenting  a  chamber  by  baffles 
changes  the  acoustic  modes  of  oscillation  of  the  full  chamber  to 
those  permitted  by  the  baffle  cavities.  Even  standing  modes  can  be 
altered  in  type,  and  in  frequency  of  oscillation,  from  those  permit¬ 
ted  for  the  full  chamber.  This  is  important  as  larger  and  larger 
rocket  engines  are  considered  since  as  a  rocket  engine  size  increases 
the  natural  frequencies  of  the  chamber  decrease.  However,  it  is  well 
known  that  the  higher  the  frequency,  the  more  stable  will  be  the 
engine  ^  .  This  is  caused  by  the  frequency-dependent  energy  feed¬ 
back  to  an  acoustic  wave  from  the  combustion  process.  It  appears 
that  above  a  certain  frequency  (dependent  upon  propellant  and  injector 
type)  this  feedback  is  insufficient  to  sustain  the  wave  by  overcoming 
the  damping  due  to  the  exit  nozzle.  Therefore,  if  baffles  really  do 
alter  the  acoustic  nature  of  the  chamber  in  the  manner  described  above 
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(note  that  there  are  some  observed  exceptions  to  this  postulate),  it 
is  clear  that  a  sufficiently  large  number  of  cavities  can  increase 
the  frequencies  of  the  allowable  modes  to  a  point  where  combustion 
Instability  is  impossible  in  the  transverse  modes. 

If  it  is  assumed  that  this  design  principle  is  adopted,  it 
is  reasonable  to  ask  if  there  is  an  optimum  baffle  pattern  or  shape. 

Of  course,  "optimum"  must  be  defined.  By  this  it  will  be  meant 
that 

1.  The  natural  fundamental  frequency  of  the  cavities  is 
as  high  as  possible  for  any  given  number  of  cavities. 

2.  The  baffles  must  not  have  a  complicated  shape  that 
would  entail  manufacturing  problems  nor  cause  exces¬ 
sive  engine  weight. 

3.  There  will  be  a  dispersive  device  in  the  baffle  shape 
such  that  given  an  oscillation  at  the  natural 
frequency  of  the  cavity,  the  amplitude  will  not  be 
uniform  over  the  cavity  but  will  be  less  than  that 

at  the  point  of  maximum  amplitude. 

It  will  be  shown  that,  given  the  above-mentioned  assumption  concerning 
the  effect  of  a  baffle,  a  shape  meeting  the  specifications  does  exist. 

Analysis 

The  effect  of  the  combustion  process  in  determining  the 
natural  frequency  of  the  transverse  oscillations  is  of  the  order  of 
the  mean  Mach  number  of  the  flow  and  are  considered  negligible 
compared  to  unity.  This  would  be  a  small  correction  even  for  large 
Mach  numbers^  .  For  transverse  oscillations,  variations  in  the 
longitudinal  direction  are  also  of  this  same  order  and  are  neglected 
here.  A  linearized  analysis  gives  a  satisfactory  first  approximation 
to  the  natural  frequency  of  oscillation  within  the  cavity.  The 
oscillation  is  described  by  the  wave  equation 
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where  <p  is  the  velocity  potential,  t  is  the  time  coordinate 
and  x  and  y  are  the  transverse  dimensions  in  a  cartesian  system. 
The  velocity  potential  has  been  nondimens ionalized  with  respect  to 
steady-state  speed  of  sound  times  chamber  diameter,  time  with  respect 
speed  of  sound  divided  by  chamber  diameter  and  space  dimensions  with 
respect  to  chamber  diameter. 

A  solution  is  assumed  to  be  of  the  form 

0-  §Uy)e.,U 

Substitution  into  Equation  (1)  yields 

+  (2) 

Now,  the  frequency  oj  would  be  determined  by  the  shape  of  the 
cavity  by  means  of  the  application  of  boundary  conditions  to  Equation 
(2).  The  general  class  of  cavities  considered  here  is  shown  in  Figure 
B-l.  One  cavity  dimension  is  much  larger  than  the  other.  The  cavity 
may  or  may  not  be  synmetric  about  the  maximum  width  line.  For  con¬ 
venience,  R(0)  is  taken  as  zero  and  G  is  defined  as  the  maximum 
width  such  that  Y(0)  *  1  .  £«•  1  is  the  important  character¬ 

istic  of  this  type  of  configuration. 

The  gradients  in  the  y-direction  are  much  larger  than  the 
gradients  in  the  x-direction.  It  is  convenient,  therefore,  to  scale 
the  y-direction  by  the  factor  G  such  that 

J-C>1  tRW) 

Note  that  this  technique  was  originally  applied  to  the  problem  of  a 

36 

long,  slender  vibrating  membrane J  . 

Under  this  scale  transformation  Equation  (2)  becomes* 

u) 

The  boundary  conditions  state  that  the  normal  velocity  at 
the  baffles  (which  are  the  cavity  walls)  be  zero.  This  means,  as 
can  be  seen  from  Figure  B-l,  that 

*Terms  of  the  order  of  and  have  been  assumed  negligible  in 

the  region  of  interest  (vicinity  of  maximum  width).  This  is  con¬ 
servative  in  that  the  inclusion  of  these  terns  shows  that  the  design 
is  still  more  favorable. 
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This  may  be  transformed  to  the  following: 

—  C  ^  along  h  ■  O 

J  *1  Jx  TT  6  1 


«- 


*  #  #  *  c 


along  (4) 


In  solving  Equation  (3)  with  the  boundary  conditions  (4),  the  following 
form  of  the  solution  is  assumed 

§(*.n)-  eyu,¥ 

Now,  Equation  (3)  becomes 


tf.(jf)'  tE'u'-t' 

.U*  L _ J ....  -  J  +  1  -  //.  \  .  •»  J 


and  the  boundary  conditions  (4)  are  now 


(5) 


j?  -  «• 


n*  ° 


-  c  if  ♦£* #  ir  alon8 

The  following  series  solutions  for  Y  and  CO  are  assumed: 


V  «  C~‘  V,  4  HJ  ♦  c'V; 


(6) 


CO  «•  C  CJ  ,  4  CO0  4  C  CO,  4 _ 

These  series  are  substituted  into  Equation  (5)  and  we  find  after 
separation  that 

=  o 


HL 


which  Cleans  that  the  lowest  order  coefficient,  V.  p  is  a  function 
of  X  only.  The  separated  equation  which  governs  Y  is  found 
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to  be 


4h  alon8  >,“o 

*long  >i-Y<0 

The  final  transformation  to  be  used  is  ty  ■  log  8  .  This 
transforms  Equation  (7)  as  follows: 

♦  [<  ♦(#•/]  e-o 

The  boundary  conditions  (8)  are  written 


(7) 


(8) 


4^-  =  (^f-  ^)e  alon*  h-o 

*-  -  along  Vc») 

Note  that  x  appears  only  as  a  parameter  above,  so  that  the 
equation  may  be  treated  as  an  ordinary  differential  equation.  The 
lowest  eigenfunction  is 


arc  tan 


e=  [^T  " 

with  the  governing  relation  for 


«Hf)] 

y as  follows: 


‘'-♦(ftT-flr 


(9) 


(10) 


In  order  to  prevent  the  solution  from  going  to  infinity,  the  positive 

JY 

value  of  J  "*  is  taken  for  negative  values  of  x  while  the  negative 
slope  is  taken  for  positive  values  of  x  .  This  means  (0)  «  o  • 
The  differential  equation  for  v.  ^  cannot  be  solved  with  a 
general  Y(x)  .  However,  by  using  a  Taylor  series  for  Y(x)  the 
equation  may  be  solved  for  a  small  region  around  the  maximum  width 
point.  Noting  that  o  ,  we  say 


Yu) 


Y(o)  +  -T  fro)  **  +  0  CxO 
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Since  (0)  =  O  »  Equation  (10)  yields  O).,  x  ik-  Therefore, 
Equation  (10)  may  be  rewritten  at  follows 


Jx  “  (Y^O* 


TTX 


2 

x  . 


Note  this  means  that  the  phase  in  Equation 
Integrating,  we  have 


F 


(O) 


(9) 


is  of  order 


The  constant  of  integration  is  neglected  since  it  would  only  be  a 
multiplicative  factor  on  the  final  solution.  Tracing  back  through 
the  series  of  transformations  and  noting  that  Y(0)  -  1  we  find 


_  c-'V^V.'OCO 

S*  e 


e*fOCe) 


w, 


oco-  e  oco 


G*C  TT^  ♦  OCC) 


(ID 

2 

Terms  of  order  x  were  neglected  in  writing  the  argument  of  the 
cosine  function  since  their  effect  is  of  order  larger  than  £,  only 
outside  of  the  region  of  large  amplitude  as  shown  by  the  exponential 
term. 


The  physical  interpretation  of  Equation  (11)  is  simple. 

The  amplitude  of  the  oscillation  is  a  maximum  at  the  maximum  width 
line  (x  -  0)  and  decays  rapidly  in  an  exponential  manner  in  either 
direction  from  this  line  (the  amplitude  is  only  significant  in  a 
region  of  size  equal  to  vr  ).  This  decay  occurs  due  to  the 
dispersive  mechanism  provided  by  the  narrowing  of  the  cavity  away 
from  the  maximum  width  line.  The  difference  in  the  curvatures  (at 
the  maximum  width  line)  of  the  two  baffles  surrounding  the  cavity 
is  most  critical  in  determining  this  decay  rate.  The  larger  the 
difference  (or  the  larger  the  rate  of  "narrowing"),  the  greater  the 


(0) 


appearing  in  the  ex* 


decay  rate  as  shown  by  the  factor) 
ponential  term. 

The  natural  frequency  was  found  to  be 


This  means  the  frequency  is  approximately  determined  by  the  maximum 
width  £  and,  therefore,  is  very  large  when  C<<  1* 

Possible  Configuration 

One  possible  shape  is  shown  in  Fig.  B-2  (assuming  that 
three  baffles  are  sufficient  to  raise  the  frequency  of  the  acoustic 
mode  to  acceptable  values).  Circular  shapes  are  one  special  case 
of  this  elliptical  shape.  (It  should  not  be  construed  that  elliptical 
chambers  are  preferred  or  even  practical.)  The  frequency  is  maxi* 
mized  by  making  four  cavities  of  minimum  width.  This  can  be  done 
by  locating  the  maximum-cavity-width  centerline  transverse  to  the 
minimum  width  axis  of  a  symmetric  chamber.  The  difference  in  radii 
of  curvature  at  the  maximum  width  points  can  be  adjusted  equally 
between  the  cavities  or  adjusted  at  will.  The  only  constraint  on 
the  rest  of  the  shape  is  that  there  is  never  another  point  at  which 
the  maximum  width  is  equalled  or  exceeded.  There  are,  therefore,  an 
infinity  of  shapes  which  are  theoretically  equally  acceptable.  Of 
course,  the  analysis  requires  £<<1  and  how  good  the  theory  is 
for  C  not  too  small  can  only  be  determined  by  experiment.  Clearly, 
it  will  fail  for  £  «  1/2  since  a  well  defined  acoustic  mode 
exists  for  the  half  chamber  if  the  full  chamber  is  nearly  circular. 
Discussion  and  Conclusions 

Under  the  assumption  that  the  purpose  of  a  combustion 
chamber  baffle  is  to  segment  a  chamber  into  cavities  with  higher 
natural  frequencies  of  gas  oscillation  than  those  of  the  full 
chamber,  an  interesting  result  from  wave  theory  and  its  application 
have  been  presented.  Baffle  design  criteria  have  been  developed 
that  are  related  to  the  specification  of  two  quantities  at  the 
maximum  width  point  of  each  cavity,  i.e.,  this  maximum  width  and 
the  difference  in  the  radii  of  curvature  of  the  two  cavity  walls  at 
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this  point.  With  this  concept  it  is  possible  to  maximize  the  acoustic 

frequencies  allowable  in  the  chamber,  do  this  with  a  smooth  somewhat 

arbitrary  shape,  and  to  provide  a  degree  of  dispersion,  i.e.,  the 

amplitude  decays  with  distance  from  the  maximum  width  point. 

The  maximization  of  frequency  is  known  to  have  beneficial 

effects  concerning  combustion  instability.  The  arbitrariness  of 

0 

shape  away  from  the  maximum  width  point  allows  freedom  concerning 
manufacturing  and  weight  problems.  The  dispersive  mechanism  is 
Important  if  the  amplitude  level  is  important  in  determining  stability 
characteristics  of  an  engine.  However,  even  given  the  assumption 
concerning  the  baffle  effect  it  should  be  remembered  that  C  «•  1  is 
required  and  the  meaning  of  "very  much  less  than"  can  only  be 
determined  experimentally.  The  theory  is  at  best  an  asymptotic 
representation  of  the  true  state  of  affairs  and  should  not  be  blindly 
applied.  Since  the  allowable  number  of  baffles  in  an  engine  is 
necessarily  limited  (the  magnitude  of  C  is  limited)  care  must  be 
taken  in  Interpreting  experimental  results. 

Nomenclature 

profile  of  baffle 
time 

transverse  dimensions 
width  of  cavity 
maximum  width  of  cavity 
translated  y-dimension 
defined  in  Appendix 
velocity  potential 
amplitude  of  velocity  potential 
defined  in  Appendix 
frequency 
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Mixture  ratio  distribution  at  3"  from  the  injector 
face  for  the  (6x2)  .2  L  on  L  I.4T-  0  injector ;  9"  dio. 
chamber,  7"  dio.  injector,  1.4  mixture  ratio  flow. 
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Mixture  ratio  and  relative  mass  distribution  at  3 
inches  from  the  injector  face  along  the 
circumference  of  the  propellant^  inject  ion  circle 
for  a  ( 6  x 2 ).  2  L  on  L  I.  4  T  -  0  injector  t  1.4 
mixture  ratio . 
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Figure  IS -16 


Mixture  ratio  distribution  at  3"  from  the  injector 
face  for  the  (6x2)  .2  L  on  L  I.4T-  0  injector,  9"  dio. 
chamber,  7"  dia .  injector,  2.4  mixture  ratio  flow. 
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Figure  ET- 1  7 
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Mixture  ratio  and  relative  mass  distribution  at  3 
inches  from  the  injector  face  along  the 
circumference  of  the  propellant  injection  circle 
for  a  ( 6  X  2 )  .  2  L  on  L  1.4  T-0  injector  ,  2.4 
mixture  ratio . 
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Figure  IS -18 


Variable  -  length  ,  square  -  motor  ,  4x4  LonL  Type  m  injector, 

8001b  thrust,  300  psia  ,  15-IK  pulsed  0~0 

LOX,  Alcohol 
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Performance  vs  chamber  length  for  the  Type  UT  -  4  spud 
injector  with  LOX  orifices  of 
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Figure  IE~26 


Pseudo "  rocket  motor  for  vapor  displacement  studies 


Figure  IS- 27 
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Schematic  of  resonating  chamber 


Figure  UT-34 


Figure  JJT-35 


Resonating  chamber ,  light  receiving ,  Q  light  source  systems 
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Figure  IY-36 
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Figure  ET-39 


Typical  photomicrographs  of  spray  in  oscillating  flow  field 


Number-size  distributions  based  on  photographs  shown  on  Figure  IF  -39 
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Figure  IF-40 
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Schematic  of  the  burning  droplet  wake 
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Neutral  stability  line. 
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Figure  I U-2 
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Geometrical  configuration 
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